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ABSTRACT 

This  paper  discusses  the  electrical  and  magnetic  design  and  analysis  of  a  permanent 
magnet  generation  module  for  naval  applications.  Numerous  design  issues  are  addressed  and 
several  issues  are  raised  about  the  potential  improvements  a  PM  generation  system  can  offer.  A 
proposed  16  MW  PM  generation  module  design  is  presented  along  with  a  detailed  design 
methodology. 

Eighty  different  machines  and  power  conversion  modules  are  sized,  designed,  and 
analyzed  with  a  final  design  selected.  Specifically,  sizing  and  detailed  machine  design  and 
analysis  is  performed  examining  the  effects  of  numerous  parameters  including  number  of  phases, 
number  of  poles,  magnetic  geometry,  machine  dimensions,  and  material  types.  Analytical 
models  are  developed  to  study  rotor  losses  caused  by  stator  winding  time  and  space  harmonics 
and  slot  space  harmonics.  Power  electronics  and  conversion  modules  to  connect  the  high-speed 
generator  to  a  DC  distribution  system  are  designed  and  analyzed.  In  depth  simulation  of  the 
eighty  complete  systems  is  performed  using  the  software  programs  MATLAB  (Version  12.0, 
Mathworks)  and  PSIM  (Version  6.0,  Powersim,  Inc.). 

The  16  MW  permanent  magnet  generation  module,  consisting  of  the  generator  and 
associated  power  electronics,  provides  an  excellent  alternative  to  traditional  wound  rotor 
synchronous  machines.  The  final  design  offers  significant  reductions  in  both  weight  and 
volume.  Specifically,  it  is  estimated  that  the  PM  generation  module  has  a  7x  reduction  in 
volume  and  a  lOx  reduction  in  weight  compared  to  similarly  rated  wound  rotor  systems.  These 
reductions  can  provide  flexibility  to  naval  architects  since  power,  weight,  and  volume  are 
integral  parts  of  the  design  and  construction  processes.  However,  further  study  is  necessary  to 
verify  the  PM  generation  modules  thermal,  structural,  and  mechanical  performance. 
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Title:  Associate  Professor  of  Naval  Construction  and  Engineering 
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Chapter  1  Introduction 


1.1  Purpose 

The  purpose  of  this  thesis  is  to  design  and  analyze  a  permanent  magnet  generator  and 
power  module  for  naval  applications.  When  deciding  whether  to  implement  an  electrical 
technology  or  component  onto  a  naval  vessel,  the  size,  weight,  and  cost  are  the  major  factors  for 
successful  integration.  Significant  performance  improvements  and  cost  reduction  of  power 
electronics,  coupled  with  the  availability  and  decreasing  cost  of  high  energy  permanent  magnet 
(PM)  materials  makes  PM  generators  attractive  for  naval  usage.  These  machines  offer  numerous 
desirable  features,  including  light  weight,  small  size,  simple  mechanical  construction,  easy 
maintenance,  good  reliability,  and  high  efficiency  [1]. 

Before  analysis  of  a  generator  can  begin,  it  must  be  properly  designed  for  typical  naval 
power  requirements.  This  involves  sizing  the  generator  along  with  designing  the  associated 
power  electronics  for  connecting  the  machine  to  the  distribution  system.  A  specific  concern 
associated  with  PM  generators  is  possible  inefficiencies  and  excessive  heating;  in  particular  rotor 
losses  caused  by  space  and  time  harmonics  during  the  energy  conversion  processes.  The 
optimum  machine  design  is  one  that  delivers  the  required  power  through  a  matching  process 
between  the  generator  and  the  power  electronic  converter  [2], 

1.2  Problem 

The  Navy’s  commitment  to  develop  an  integrated  electric  power  system  for  the  next 
generation  warships  offers  the  expectation  of  using  the  installed  generation  capacity  to  power 
ship  propulsion,  advanced  weapons  components,  and  high  power  combat  control  systems  [3].  As 
these  electrical  loads  increase,  it  becomes  increasingly  important  to  efficiently  utilize  installed 
power  as  well  as  develop  smaller,  effective  power  generation  systems.  Navy  ships  are  extremely 
high  performance  systems  and  therefore  power  and  weight  considerations  are  integral  parts  of 
the  design  process. 

The  life  cycle  of  a  navy  ship  is  on  average  2-3  times  longer  than  a  commercial  ship  and 
therefore  navy  ships  undergo  excessive  modernization  and  upgrades  throughout  their  service  life. 
Many  of  the  newer  components  have  significantly  higher  power  requirements  than  the  originals 
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putting  a  burden  on  the  power  generation  system.  Therefore,  the  Navy  is  moving  toward 
designing  all  electric  ships  with  integrated  power  systems  (IPS)  and  increased  power  generation 
that  can  be  efficiently  managed  to  meet  future  demands. 

The  effective  integration  of  electrical  power  in  future  naval  ships  requires  the 
development  of  technologies  that  ensure  volume  and  mass  reduction  in  critical  mechanisms. 
Military  ships  require  higher  power  density  components,  impose  more  stringent  acoustic  and  EM 
signature  requirements,  and  subject  systems  to  harsher  environments  than  commercial 
applications  [4].  Rotating  generators,  coupled  with  prime  movers,  need  to  be  lighter  in  weight 
and  higher  in  power  density.  High-speed  PM  generators  provide  a  substantial  reduction  in  size 
and  weight  making  them  a  logical  choice  for  naval  applications.  Currently,  the  Navy  has  not 
designed  or  built  a  high  power  (megawatt)  PM  generator  and  therefore  the  need  exists. 

In  conjunction  with  constructing  a  PM  generator,  a  DC  bus  architecture  is  one  of  the 
preferred  schemes  for  the  future  [5].  DC  power  distribution  systems  can  offer  a  size  and  weight 
reduction  compared  to  high-power  AC  systems  [3].  With  a  DC  bus  distribution,  the  PM 
generator  can  be  optimized  independent  of  producing  60  Hz  frequency  as  required  in  the  past. 
However,  conversion  of  the  high-frequency  generator  AC  output  to  DC  requires  power 
electronics.  Rectification  of  AC  to  DC  presents  the  problem  of  creating  harmonics  in  the  input 
current  which  are  then  reflected  back  onto  the  generator  causing  rotor  losses.  In  addition,  the 
generator  produces  space  harmonics  which  also  produce  losses  in  the  rotor.  Therefore,  the  PM 
generator  and  power  electronics  module  (PEM)  need  to  de  designed  and  optimized  to  deliver 
constant  DC  power  while  minimizing  machine  losses. 

1.3  Background 

1.3.1  History 

The  Navy  has  designed  and  built  electric  ships  since  the  early  part  of  the  20th  century. 

The  original  advantages  perceived  for  electric  ships,  superior  performance,  reduced  manning, 
arrangement  flexibility,  and  fuel  efficiency,  are  still  relevant  today  [6],  In  the  early  part  of  the 
20th  century,  diesel-electric  submarines,  small  surface  ships,  and  some  battleships  and  carriers 
had  electric  propulsion.  By  the  late  1940s,  mechanical  drive  systems  became  popular  because  of 
improvements  in  metallurgy  and  manufacturing.  However,  the  capability  of  mechanical 
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transmissions  reached  its  limits  in  the  late  20  century  spawning  a  renewed  interest  in  electric 
drive  and  integrated  power  systems  for  military  applications. 

For  the  last  thirty  years,  the  commercial  industry  has  designed  and  operated  integrated 
power  systems  (IPS).  An  IPS  is  a  ship  architectural  paradigm  in  which  the  ship’s  power  and 
propulsion  are  provided  by  a  common  electrical  distribution  system  instead  of  having  a  separate 
mechanical  drive  for  propulsion.  For  military  applications,  an  IPS  provides  numerous  benefits: 

•  Decreased  life  cycle  costs  because  of  increased  fuel  economy  and  efficiency  (a  Navy  ship 
with  IPS  may  consume  10-25%  less  fuel  than  a  similar  ship  with  mechanical  drive  [7]) 

•  Increased  ship  design  and  arrangement  flexibility  since  the  ship  is  not  limited  to  having  a 
long  mechanical  shaft  line 

•  Reduced  system  complexity 

•  Higher  degree  of  modular  design  using  power  components 

•  Broad  industrial  base  for  implementing  IPS  design 

•  High  levels  of  automation  and  control 

•  Increased  power  available  for  non-propulsion  uses  since  a  Navy  ship  spends  a  large 
portion  of  its  time  operating  at  low  propulsion  levels  (approximately  95%  of  the  time) 

•  Increased  stealth,  survivability,  and  payload 

Since  an  IPS  provides  power  for  both  the  ship  loads  and  propulsion,  larger  generation  capability 
is  required.  As  the  Navy  builds  its  new  ships  with  IPS  architectures,  compact,  high-power 
generation  systems  must  be  examined  to  help  facilitate  implementation  of  the  new  designs. 

1.3.2  Power  Generation  &  Distribution 

Almost  all  naval  core  power  generators  to  date  are  air-cooled  50/60  Hz  machines  that  are 
military  derivatives  of  commercial  generators  and  are  therefore  size  and  weight  excessive  [8].  A 
typical  turbine  generator  system  is  shown  in  Figure  1  and  Table  1  contains  nominal 
characteristics  for  several  current  commercial  and  naval  generators. 
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Table  1:  Examples  of  Current  Generator  Characteristics 
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With  the  advent  of  high-power,  cost-effective  power  electronics,  it  is  no  longer  necessary  to 
generate  power  at  50/60  Hz  so  generators  can  be  optimized  independent  of  frequency.  High¬ 
speed,  power-dense  generators  become  the  logical  choice  for  naval  purposes.  Permanent  magnet 
machines  are  ideal  for  this  high-speed  application  due  to  their  simple  structure  and  high  power 
density  [9]. 

Since  the  generator  can  now  be  designed  to  produce  higher  frequencies,  distribution 
architectures  are  not  limited  to  being  60  Hz.  Either  a  high  frequency  AC  system  or  a  DC  system 
can  be  designed  through  the  use  of  power  electronics  modules,  with  the  DC  distribution  being 
preferred  because  of  its  advantages  in  size  and  weight.  In  this  thesis,  a  high-power  DC  zonal 
architecture  is  assumed  using  solid  state  converters  to  generate  AC  where  needed.  Each  zone  is 
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electronically  isolated  from  the  other  zones  with  automatic  fault  detection  and  reconfiguration  to 
provide  continuous  power  during  damaged  conditions.  The  PM  generator  sets  and  power 
electronic  conversion  modules  serve  as  the  backbone  thereby  providing  a  reliable  power  system 
for  navy  ships. 

1.4  Scope 

The  scope  of  this  thesis  is  limited  to  the  PM  generator  and  associated  power  electronics 
AC-DC  conversion  module.  The  following  is  accomplished: 

•  Determine  the  electrical  power  requirements  for  a  Navy  IPS  ship  in  order  to  properly  size 
the  generator 

•  Compare  typical  wound  rotor  machine  design  to  a  permanent  magnet  design  to  determine 
applicability  for  IPS  applications 

•  Conduct  material  analysis  and  selection  for  the  generator  design 

•  Perform  initial  PM  generator  detailed  design 

•  Design  the  power  electronics  conversion  module  to  perform  high-power  AC-DC 
conversion 

•  Conduct  detailed  analysis  of  rotor  losses  of  the  PM  generator,  in  particular  those  caused 
by  time  and  space  harmonics 

•  Perform  numerous  iterations  of  machine  and  power  electronics  designs  to  develop 
optimized  generation  scheme 
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Chapter  2  Power  Requirements  and  Machine  Selection 


2.1  Machine  &  Module  Requirements 

On  board  Navy  ships,  electricity  is  used  to  provide  power  to  virtually  all  components, 
including  mission  systems,  support  systems,  combat  systems,  and  communications  systems.  In 
addition,  as  ships  continue  to  be  upgraded  and  modernized,  more  power  is  needed  for  newer 
combat  systems  and  weapons  components.  Most  current  naval  platforms  have  some  form  of 
mechanical  propulsion  system  with  separate  ship  service  electrical  generators  supplying  the 
ship’s  power.  With  an  IPS  ship,  the  ship’s  generators  provide  power  for  propulsion  and  the 
ship’s  service  loads,  and  through  proper  utilization,  power  is  efficiently  managed. 

To  properly  size  the  PM  generator,  the  power  requirements  must  be  identified  and 
therefore  a  typical  load  list  is  developed  for  an  IPS  naval  ship  and  is  included  in  Appendix  A. 
The  overall  power  requirement  for  the  generator  is  16  MW.  Since  size  and  weight  are  important 
factors  and  the  generator  can  be  optimized  independent  of  frequency,  high-speed  operation  and 
maximum  power-density  are  desired.  Therefore,  the  highest  possible  speed  is  selected  while 
ensuring  the  PM  generator  is  compatible  with  both  gas  turbines  and  steam  turbines. 

Traditionally,  gas  turbines  run  at  much  higher  speeds  than  steam  turbines  causing  the 
steam  turbines  to  be  more  limiting.  From  information  collected  from  the  Elliot  Turbomachinery 
Company,  Inc.,  16,000  RPM  is  approximately  the  highest  speed  steam  turbine  that  can 
reasonably  be  constructed  at  the  megawatt  power  level  [13].  Therefore,  to  provide  a  degree  of 
conservatism,  13,000  RPM  is  selected  for  the  nominal  design  speed  for  the  PM  generator. 

The  power  electronics  module  (PEM)  converts  the  AC  voltage  from  the  generator  to 
700  VDC.  Overall,  the  PEM  and  generator  must  be  designed  so  that  losses  suffered  by  the 
permanent  magnets  on  the  generator  rotor  are  minimal.  Table  2  lists  the  general  requirements  for 
the  entire  system. 

2.2  Machine  Selection 

Military  ships  require  high  power  density  components  and  improved  acoustic  and 
electromagnetic  signature  requirements  while  subjecting  systems  to  harsh  environments  [14].  It 
is  therefore  important  to  ensure  the  power  generation  system  is  capable  and  efficient. 
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Table  2:  General  Machine  &  Module  Requirements 


1 

1  .  .  . . .  ■  ■  '  : 

. .  . 

SI  IS  11  1: 

Generator  Power 

16  MW 

Generator  Speed 

13,000  RPM 

PEM  Output  Voltage 

700±5  VDC 

PEM  Output  Ripple 

0.7  VDC  (@  16  MW) 

Generator  Rotor  Losses 

Minimal 

2.2. 1  Permanent  Magnet  versus  Wound  Rotor 

Reducing  the  size  and  weight  of  ship’s  turbine  generator  sets  offers  significant 
advantages  to  naval  architects.  Replacing  older  generators  with  lightweight  ones  could  make  it 
possible  to  decrease  the  size  of  some  generator  sets  by  as  much  as  50%  [15].  PM  generators 
therefore  become  an  attractive  alternative  compared  to  wound  rotor  machines  because  of  the 
availability  and  decreasing  cost  of  high  energy  PM  materials  along  with  improved  power 
electronics. 

A  wound  rotor  generator  normally  consists  of  armature  windings  on  a  stationary  stator 
frame  with  field  windings  on  an  inner  rotor.  The  rotor  is  turned  by  a  prime  mover,  usually  a  gas 
or  steam  turbine,  and  current  is  supplied  to  the  field  windings  through  brushes  or  a  brushless 
exciter.  As  the  current-carrying  field  windings  rotate  past  the  stator  windings,  current  is 
produced  in  the  stator  windings  through  Faraday’s  Law.  An  example  of  a  wound  rotor  machine 
is  shown  in  Figure  2  [8]  and  Figure  3  [16]. 


Figure  2:  Example  of  Wound  Rotor  Generator 
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Figure  3:  Cross  Section  of  Wound  Rotor  Generator 


Wound  rotor  generators  have  been  the  backbone  of  power  generation  for  the  U.  S.  Navy  because 
they  are  a  proven  technology  that  is  battle  tested.  They  offer  steady  voltage  regulation  using  the 
field  windings,  a  large  air  gap  for  producing  the  rotor  magnetic  flux,  low  fault  currents,  and  high 
power  capabilities. 

However,  the  machines  tend  to  be  complex,  weight  excessive,  and  they  require  field 
windings  which  limit  design  alternatives.  There  are  several  drivers  which  cause  these  problems 
to  occur.  First,  to  generate  the  necessary  magnetic  flux  levels,  wound  rotor  generators  have  large 
pole  pitches  to  support  the  required  field  windings.  These  pole  pitch  windings  in  turn  require 
larger  end  turns  and  thick  back  iron  to  support  the  magnetic  flux,  both  of  which  contribute  to 
increased  size.  Second,  because  of  the  winding  losses  in  the  rotor,  large  cooling  systems  can  be 
required  thus  increasing  the  number  of  support  components. 

High  speed  generators  offer  a  reduction  in  machine  size  and  weight  because  as  a 
machine’s  speed  increases,  its  size  decreases  for  a  given  output  power.  The  PM  generator  is 
ideal  for  high-speed  applications  because  of  its  simple  structure  and  high  power  density  [9].  In  a 
PM  generator,  the  rotor  field  windings  are  replaced  by  permanent  magnets  which  do  not  require 
additional  excitation.  As  the  permanent  magnets  are  rotated  by  the  prime  mover,  current  is 
produced  in  the  stator  windings.  An  example  of  a  PM  generator  is  shown  in  Figure  4  [8]. 
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Figure  4:  Example  of  PM  Generator 

PM  generators  offer  several  advantages:  they  have  no  rotor  windings  so  they  are  less 
complicated;  they  have  high  efficiencies;  the  gap  field  flux  is  not  dependent  on  large  pole  pitches 
so  the  machine  requires  less  back  iron  and  can  have  a  greater  number  of  smaller  poles;  and  they 
usually  require  smaller  and  fewer  support  systems.  Assuming  the  same  flux  density  and 
circumferential  arc,  doubling  the  number  of  poles  produces  the  same  radial  flux  but  requires  half 
the  stator  core  thickness,  as  shown  in  Figure  5. 


Figure  5:  Flux  vs.  Number  of  Poles 
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However,  PM  generators  also  have  some  disadvantages.  They  do  not  possess  field 
excitation  control  and  therefore  voltage  regulation  can  be  problematic.  This  can  be  corrected  by 
using  external  voltage  control  such  as  large  capacitor  banks  or  power  electronics,  as  well  as 
choosing  the  turns  on  the  stator  winding  properly  to  produce  the  anticipated  required  nominal 
voltage.  Additionally,  since  the  permanent  magnet  fields  cannot  be  turned  off,  there  exists  the 
risk  of  excessive  currents  in  the  event  of  an  internal  fault.  This  problem  can  also  be  solved 
through  the  design  of  the  turbine  governor  and  controller  or  dynamic  braking.  Overall,  the 
advantages  of  the  PM  generator  over  the  traditional  wound  rotor  generator  make  it  a  better 
alternative  for  high-speed  navy  applications.  A  summary  comparing  the  different  designs  is 
given  in  Table  3. 


Table  3:  Comparison  of  Wound  Rotor  and  PM  Generators 


wmslmm 

jjjj|  f| j ttiijMfe j 1 1 : 

Wound  Rotor 

Steady  voltage  regulation  with 
field  windings 

Weight  excessive 

High  power  capabilities 

Large  size 

Large  air  gap  for  flux 

Rotor  windings  &  associated 
losses 

Low  fault  currents 

Large  support  systems 

Proven,  robust  design 

Permanent  Magnet 

Less  complicated 

Lack  of  inherent  voltage 
regulation 

Reduced  size  and  weight 

Potential  fault  currents 

High  efficiency 

Magnet  losses 

No  excitation  supply  or  field 
windings 

High  speed  applicability 

2.2.2  Type  of  Permanent  Magnet  Machine 

There  are  numerous  layout  possibilities  for  permanent  magnet  machines  and  only  the 
most  common  are  discussed  here.  These  include  radial  flux  inner  rotor,  radial  flux  outer  rotor, 
and  axial  flux  designs.  In  most  PM  machines,  flux  crosses  from  the  rotor  to  the  stator  in  the 
radial  direction  [17].  The  first  type,  the  radial  flux  inner  rotor  design,  is  the  closest  configuration 
to  the  classical  AC  synchronous  generator.  An  example  of  this  design  is  shown  in  Figure  6  [17]. 
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Figure  6:  Example  of  Inner  Rotor  PM  Machine 

In  this  type  of  machine,  the  windings  are  placed  on  the  stator,  either  in  slots  or  in  a 
slotless  ring,  and  the  magnets  are  surface  mounted  on  the  rotor  or  buried  in  the  rotor.  Buried 
magnet  designs  often  result  in  rotors  that  are  larger  than  equivalent  surface-magnet  machines 
with  high-energy  magnets  [18].  Buried  magnet  machines  can  also  have  significant  structural 
issues  in  high-power  applications  [19].  When  the  magnets  are  surface  mounted  and  the  machine 
is  operated  at  high  speed,  the  magnets  are  often  secured  with  a  retaining  device  made  of  either 
alloy  steel  or  carbon-fiber.  Overall,  the  inner  rotor  machine  possesses  high  torque/power 
capability  and  good  heat  conduction  and  cooling  properties  making  it  ideal  for  high-speed, 
higher-power  applications  [18]. 

The  radial  flux  outer  rotor  machines  are  commonly  used  in  hard  disk  drives,  small 
computer  ventilation  fans,  and  some  blowers.  This  type  of  design  is  very  efficient,  low-cost, 
easy  to  manufacture,  and  applicable  for  low-power  applications  [18].  It  is  the  opposite  of  the 
inner  rotor  because  the  stator  in  on  the  inside  with  the  rotor  and  magnets  on  the  outside.  A  cross 
section  of  an  outer  rotor  machine  is  shown  in  Figure  7  [17]. 
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Figure  7:  Example  of  Outer  Rotor  PM  Machine 


The  axial  flux  machine  is  significantly  different  than  the  previous  two  because  flux  flows 
in  the  axial  direction  vice  radial  direction  and  the  windings  are  oriented  radially  vice  axially  (see 
Figure  8  for  an  example  diagram  [20]). 


The  main  advantages  of  this  design  are  their  low  cost,  flat  shape,  and  smooth  rotation.  However, 
if  axial-flux  machines  are  operated  at  high  speeds  (above  1000  RPM),  eddy-current  losses  and 
heating  can  become  excessive  [18].  Also,  stator  construction  is  difficult  because  it  must  be 
laminated  circumferentially.  An  example  of  this  design  is  the  turntable  for  a  record  player. 
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Overall,  because  of  its  inherent  advantages  in  heat  removal  and  cooling,  the  abundance  of 
manufacturing  capabilities,  and  its  high-power,  high-speed  applicability,  the  radial  flux  inner 
rotor  with  surface  mounted  magnets  is  selected  for  the  16  MW  PM  generator  design. 
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Chapter  3  Material  Selection  and  Machine  Initial  Design 


3.1  Material  Selection 

One  of  the  key  considerations  during  the  electromagnetic,  structural,  and  thermal  design  of 
a  permanent-magnet  machine  is  the  selection  of  the  magnet,  stator,  and  rotor  materials  [21]. 
Machine  output,  heat  rise,  weight,  and  cost  are  a  few  of  the  characteristics  which  are  directly 
influenced  by  selection  of  the  machine  materials  [22]. 

3.1.1  Permanent  Magnets 

The  size  and  performance  of  high-speed  PM  generators  depend  on  the  permanent  magnet 
material  properties  [9].  The  magnets  must  be  selected  to  provide  the  necessary  air  gap  magnetic 
field  and  ample  coercive  force  to  compensate  for  possible  damaging  effects  while  minimizing 
the  volume  of  material  because  of  cost  and  weight  considerations  [23]. 

Ferromagnetic  materials  are  the  most  common  substances  used  in  the  construction  of 
machines  and  their  properties  are  normally  described  using  B-H  curves  and  hysteresis  loops. 
These  curves  represent  an  average  material  characteristic  that  reflects  the  non-linear  property  of 
the  permeability  of  the  material  but  ignores  the  multi-valued  properties  [17],  An  example  of  a  B- 
H  curve  is  shown  in  Figure  9. 


Figure  9:  Example  of  B-H  Curve 


27 


Several  basic  magnetic  properties  are  of  critical  importance  for  the  permanent  magnets  in  a  PM 
machine: 

•  Remnant  Flux  Density  (Br):  It  is  the  value  of  the  flux  density  remaining  after 
magnetization  and  it  directly  influences  the  air  gap  flux  and  magnet  sizes. 

•  Coercivity  (He):  It  is  the  value  of  magnetizing  field  needed  to  reduce  the  flux  density  in 
the  magnet  to  zero  and  it  gives  a  first  order  estimate  of  a  magnet’s  resistance  to 
demagnetization. 

•  Energy  Product  (BHmax):  It  is  the  maximum  energy  product  of  the  magnet  and  it  is 
inversely  proportional  to  the  total  magnet  volume  required. 

•  Recoil  Permeability  (prec):  It  is  the  gradient  of  the  B-H  curve  and  it  gives  the  magnet’s 
ability  to  return  to  its  initial  magnetization  after  subjected  to  damaging  forces.  If  the 
magnet  goes  below  Hk,  then  it  will  recoil  along  a  lower  line  resulting  in  a  lower  magnet 
flux  density. 

•  Load  Line:  It  is  a  line  drawn  from  the  origin  to  the  magnet  operating  point  on  the 
hysteresis  curve  (Bm).  The  magnitude  of  the  slope  of  the  load  line  is  the  permeance 
coefficient. 

Permanent  magnet  materials  come  in  many  varieties  and  the  four  most  common  types  for 
machine  applications  are  Alnico,  Ferrites,  SmCo  material,  and  NdFeB  material.  Table  4  and 
Figure  10  show  the  characteristics  and  typical  B-H  curves  for  these  materials  [18]. 


Table  4:  Magnet  Material  Properties 
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mmm 

SUL.. 

BlIllBIl 

Coercivity  (He) 

kA/m 

40-130 

180-400 

800-1500 

800-1900 

Recoil  Permeability 

(Prec) 

1.9-7 

1.05-1.15 

Energy  Product  (BHmax) 

kJ/m3 

20-  100 

24-36 

140-220 

180-320 

Maximum  Temperature 

500-550 

250 

250-350 

100-200 

Br  Temperature 
Coefficient 

%/°c 

-0.01  to  -0.02 

-0.2 

-0.05 

-0.08  to -0.15 

The  rare-earth  magnets,  SmCo  and  NdFeB,  have  become  more  popular  for  high  performance 
applications  because  of  their  greater  power  density,  high  coercivity,  high  flux  densities,  and 
linearity  of  the  demagnetization  curves  [24], 
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Alnico  “—^-Ferrite  SmCo  . NdFeB 


Figure  10:  Typical  Magnet  B-H  Curves 

Between  the  two  rare-earth  permanent  magnets,  NdFeB  is  preferred  because  it  is  cheaper  and 
more  readily  available.  It  does  possess  some  adverse  characteristics  such  as  moderate  corrosion 
and  lower  resistance  to  temperature  effects,  but  these  can  be  controlled  using  surface  treatments 
and  proper  cooling  [25].  Therefore,  NdFeB  magnets  are  selected  for  use  in  the  PM  generator 
with  the  conservatively  assumed  values  listed  in  Table  5. 


Table  5:  Selected  Magnet  Properties 


Units  ,1;: 

flu  f:  Value  ■ 

Remanence  (Br) 

T 

1.2 

Coercivity  (Hc) 

kA/m 

900 

Recoil  Permeability  (prec) 

1.05 

Energy  Product  (BHmax) 

260 

Maximum  Temperature 

°c 

180 

Resistivity 

p£2/m 

1.43 
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3.1.2  Stator  and  Rotor  Material 

The  type  of  material  chosen  for  the  stator  and  rotor  is  important  because  it  impacts  the 
machine  losses  and  efficiency.  The  rotor  is  usually  built  from  the  same  material  as  the  stator  for 
ease  of  construction  but  is  can  be  made  of  any  economical  steel  provided  it  is  strong  enough  for 
the  given  function  [18].  No  one  material  is  optimum  for  every  application  and  the  normal 
criteria  for  selection  are  cost,  permeability,  core  losses,  and  saturation  flux.  It  is  important  that 
the  material  act  as  a  flux  guide  and  absorb  the  minimum  amount  of  magnetomotive  force  (MMF) 
so  that  the  flux  is  concentrated  in  the  air  gap.  In  addition,  the  material  should  minimize  core 
losses  including  hysteresis  and  eddy  current  losses. 

High-quality,  non-oriented,  electrical  grade  lamination  steels  are  typically  used  in  most 
machines  because  the  laminations  help  minimize  losses.  The  four  main  materials  are  low  carbon 
steels,  silicon  (Si)  steels,  nickel  (Ni)  alloy  steels,  and  cobalt  (Co)  alloy  steels.  Low  carbon  steels 
are  the  lowest  cost  and  are  used  in  high  volume  applications  where  high  core  losses  are 
acceptable.  Silicon  steels  usually  have  3%  silicon  which  increases  the  resistivity  to  reduce  eddy 
current  losses.  They  are  selected  and  specified  based  on  core  loss,  with  each  grade  (Ml 9,  M27, 
M36,  and  M43)  having  higher  core  losses  and  lower  cost  [22].  The  lamination  thickness  is  a 
tradeoff  between  cost  and  performance  and  the  most  common  sizes  are  0.014  in,  0.0185  in,  and 
0.025  in  (29  gauge,  26  gauge,  and  24  gauge). 

Nickel  alloys  are  either  49%  or  80%  nickel  and  they  have  lower  losses  than  the  silicon 
steel  but  are  much  more  expensive.  In  addition,  they  require  careful  handling  and  not  suited  for 
high  flux  density  environments  (above  0.8  T)  because  of  saturation.  The  cobalt  alloys  are  only 
used  in  extremely  high-performance  situations  such  as  military  aircraft  and  space  applications 
because  of  the  high  cost.  Table  6  summarizes  the  different  stator  materials  and  the  Ml 9,  29- 
gauge  electrical  silicon  steel  is  selected  for  the  PM  generator  because  it  is  economical,  its  thin 
laminations  minimize  losses,  and  it  has  a  saturation  flux  density  of  about  1.8  T  [2],  [18],  [22]. 
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Table  6:  Laminated  Steel  Properties 


ffgg5  ■ 

■  illl! 

If!  |§§§ 

a 

Ease  of 
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Low  Carbon 
Steel 

Fair 

Good 

Good 

Best 

0.5 

Si  Steel 

Good 

Good 

Fair 

Good 
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Thin  Si  Steel 

Better 

Good 

Fair 

Fair 
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Good 

Fair 
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3.2  Machine  Design  Parameters 
3.2. 1  Stator  Mechanical  Design 

The  stator  is  an  important  part  of  the  machine  because  it  serves  as  the  main  structural 
component,  it  provides  the  housing  for  the  armature  windings,  and  it  completes  the  flux  path  for 
the  magnetic  circuit.  The  main  consideration  in  the  mechanical  design  of  the  stator  is  whether  to 
make  it  slotted  or  slotless.  A  slotless  stator  has  the  armature  windings  located  in  the  air  gap  of 
the  machine  as  shown  in  Figure  11  [19]. 


Figure  11:  Slotless  Stator  Design 

One  of  the  advantages  of  the  slotless  construction  is  unique  winding  layouts  are  possible  to  meet 
specific  performance  goals.  Another  advantage  is  that  the  space  available  for  the  armature 
windings  increases  by  a  factor  of  about  two  since  there  are  no  stator  teeth.  This  produces  lower 
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conductor  losses  since  less  current  flows  in  each  winding.  The  flux  density  is  reduced,  however, 
because  the  effective  air  gap  is  much  larger  since  the  windings  are  in  the  air  gap.  Overall,  there 
exists  a  higher  electrical  loading  and  a  lower  magnetic  loading. 

One  disadvantage  of  the  slotless  design  is  there  are  no  good  conduction  paths  to  remove 
the  heat  generated  from  the  windings.  This  reduces  the  allowable  current  density  in  the  windings 
and  lowers  the  power  output.  Another  disadvantage  is  that  the  windings  are  directly  exposed  to 
the  rotating  flux  which  raises  the  possibility  of  additional  eddy-current  loss  in  the  conductors  and 
further  losses  due  to  circulating  currents  in  the  windings  [18].  Overall,  the  performance  of  a 
slotless  stator  is  almost  always  lower  than  that  of  an  equivalent  slotted  stator  design  and 
therefore  slotless  stators  do  not  appear  often  in  high-power  applications  [17]. 

Slotted  stators  are  the  traditional  stator  design  and  consist  of  openings  around  the  stator 
for  the  armature  windings  as  shown  in  Figure  12  [26].  The  openings  provide  rigid  housings  for 
the  conductors  and  associated  insulation. 


Figure  12:  Slotted  Stator  Design 

Stator  slots  vary  in  size  and  shape  with  the  most  common  configurations  being  rectangular  or 
trapezoidal.  In  this  paper,  the  slots  are  assumed  to  be  approximately  rectangular  as  shown  in 
Figure  13  and  contain  form-wound  windings  so  that  the  depression  width  is  the  same  as  the  slot 
top  width. 
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Figure  13:  Stator  Slot  Geometry 

Slotting  is  used  because  is  provides  a  good  way  to  achieve  a  narrow  air  gap  length  while  keeping 
the  winding  conductors  close  to  the  magnets  to  maximize  the  flux  linkage.  The  slots  also  greatly 
increase  the  surface  contact  area  between  the  windings  and  stator  steel  providing  a  path  of  low 
thermal  resistance  for  good  heat  conduction  which  is  important  for  keeping  the  windings  and 
magnets  cool. 

The  resulting  narrow  air  gap  from  the  slots  makes  the  permeance  greater  and  therefore 
the  air  gap  flux  density  greater  producing  a  more  powerful  machine.  In  addition,  the  depression 
in  the  slot  tops  help  control  parasitic  losses  in  the  rotor  by  improving  the  uniformity  of  the  air 
gap  field.  The  limits  of  the  size  of  the  slots  are  twofold:  the  magnetic  teeth  must  be  able  to  carry 
the  air-gap  flux  without  saturating  and  the  slots  must  be  large  enough  to  support  the  necessary 
current  density  in  the  windings.  Typical  limits  for  stator  current  density  are  shown  in  Table  7 
and  in  this  paper  it  is  assumed  that  the  limit  on  current  density  (J)  is  3000  A/cm  [18]. 


Table  7:  Stator  Current  Densities 


1  Cooling  Method  ’ 

3'  (A/ciii2) 

Natural  Convection 

Fan  Cooled 

800-  1200 

Liquid  Cooled 

2300  -  3200 

The  disadvantages  of  the  slots  are  that  cogging  torque  may  be  a  problem  and  it  can  be 
costly  to  insert  the  windings  if  proper  construction  techniques  are  not  used.  Overall,  however, 
slotted  designs  are  preferred  in  high-power  applications  and  therefore  a  slotted  stator  is  selected 
for  the  16  MW  PM  generator. 
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The  number  of  slots  in  the  machine  is  usually  a  function  of  the  number  of  phases  and 
windings  and  can  vary  based  on  the  application.  The  initial  design  of  the  generator  assumes  a 
three-phase  machine  but  the  number  of  phases  will  be  examined  to  determine  an  optimum 
design.  In  order  to  allow  for  high-power  operation  and  the  possibility  of  a  high  number  of 
phases,  36  slots  is  chosen  for  the  initial  generator  design. 


3.2.2  Rotor  Mechanical  Design 

For  high-speed  applications,  the  rotor  aspect  ratio,  defined  as  length-to-diameter  (L/D),  is 
a  critical  parameter.  If  it  is  relatively  low,  then  the  rotor  has  high  stiffness  and  good  dynamics 
but  a  large  diameter  which  increases  the  weight  and  makes  magnet  retention  extremely  difficult. 
Additionally,  the  centrifugal  force  on  the  surface-mounted  magnets  is  directly  proportional  to  the 
rotor  diameter  so  the  rotor  radial  size  must  not  be  excessive. 

Permanent  magnet  machines  offer  flexibility  in  selecting  pole  sizes  which  allows  for 
smaller  diameters.  They  are  therefore  ideal  for  high-speed  applications  because  they  can  have 
higher  L/D  ratios.  This  is  because  they  do  not  have  rotor  field  windings  which  have  end  turns 
necessitating  big  pole  pitches  and  large  diameters.  A  normal  L/D  ratio  for  a  wound  rotor 
machine  is  0.5  -  1.0  compared  to  1  -  3  for  a  PM  machine  [27].  Staying  close  to  these  ranges 
usually  provides  a  first  order  estimate  of  satisfactory  machine  dynamic  performance  and 
acceptably  low  vibrations  or  oscillations. 

The  rotor  radius  and  the  rotational  speed  also  determine  the  tip  speed  of  the  machine 
which  is  the  surface  velocity  of  the  rotor  (as  defined  by  Eqn  3-1). 

vtip  =  Rcom 

where  a>m  =  angular  speed  (rad/sec) 

R  =  rotor  radius  (m) 


Eqn  3-1 

For  most  rotating  machines,  the  upper  limit  on  tip  speed  is  between  100  -  250  m/s  depending  on 
the  design.  For  surface  magnet  PM  machines,  retaining  sleeves  are  sometimes  used  to  help  keep 
the  magnets  in  place  and  allow  for  higher  speeds.  These  sleeves  can  be  constructed  from  alloy 
steel,  carbon  fiber,  or  other  materials.  The  metal  sleeves  usually  provide  increased  mechanical 
performance  but  have  eddy  current  losses. 
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The  carbon  fiber  and  graphite  composite  sleeves  have  high  strength-to-weight  ratios 
which  produce  a  thin  sleeve  and  the  sleeve’s  lower  conductivities  yield  reduced  eddy  current 
losses.  However,  the  carbon  fiber  and  graphite  composite  sleeves  have  lower  temperature 
ratings  and  lower  thermal  conductivities  making  heat  removal  and  increased  cooling  for  the 
magnets  and  sleeve  important  issues  [28].  Overall,  the  use  of  a  retaining  sleeve  is  necessary  for 
the  16  MW  generator  since  it  is  operating  at  high-speed  and  this  allows  the  maximum  tip  speed 
limit  at  the  rotor  surface  to  be  200  m/s.  The  actual  material  for  the  retaining  sleeve  is  examined 
later  when  detailed  rotor  loss  analysis  is  performed. 


3.2.3  Number  of  Poles  and  Magnet  Pole  Design 

The  optimum  number  of  poles  is  a  complex  function  depending  on  a  number  of  factors 
including  the  magnet  material,  the  speed  of  rotation,  the  desired  output  frequency,  and  the 
mechanical  assembly  of  the  rotor.  An  even  number  of  poles  is  always  used  because  this  provides 
a  balanced  rotational  design.  As  the  number  of  poles  increases,  the  individual  pole  pitch  goes 
down  which  reduces  the  amount  of  stator  back  iron  needed  to  support  the  magnetic  flux.  In 
addition,  for  a  given  power/torque,  as  the  pole  number  rises,  the  required  magnet  volume 
decreases. 

Assuming  a  constant  mechanical  rotation  speed,  the  generated  electrical  frequency  is 
proportional  to  the  number  of  poles  as  shown  in  Eqn  3-2. 

N-(2p)  =  120*  f 

where  N  =  speed  (RPM) 

p  =  number  of  pole  pairs 
f  =  electrical  frequency  (Hz) 


Eqn  3-2 

If  a  PM  generator  is  going  to  be  the  source  for  a  DC  bus  through  a  rectifier  system,  a  high  pole 
number  is  desirable  because  as  the  electrical  frequency  increases,  support  components  such  as 
filter  capacitors  and  inductors  can  be  much  smaller.  Therefore,  for  a  given  rotational  speed,  one 
cheap  and  efficient  solution  is  to  have  a  higher  number  of  pole  pairs  and  frequency  [27]. 
However,  as  the  frequency  increases,  higher  stator  losses  result  because  core  losses  are 
proportional  to  frequency  squared.  In  addition,  as  the  pole  number  gets  larger,  the  number  of 
slots  per  pole  per  phase  decreases  and  can  cause  the  voltage  waveforms  to  become  less 
sinusoidal  so  all  factors  must  be  considered. 
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The  pole  arc  of  the  magnets  can  also  be  varied.  Magnets  seldom  span  the  full  pole  pitch 
because  the  flux  at  the  transition  between  north  and  south  poles  leaks  between  poles  without 
linking  the  coils  in  the  stator.  The  gaps  between  the  poles  usually  contain  non-magnet  pieces, 
such  as  soft-iron,  so  that  no  flux  crosses  over  the  air  gap  between  magnets.  A  full  pole  arc  is  0me 
=  180°E  and  produces  a  full  voltage  waveform  but  has  increased  harmonic  content.  As  the  pole 
arc  is  reduced  (up  to  20  -  30  %)  and  those  areas  are  filled  in  with  soft-iron  pieces,  the  resulting 
flux  waveform  is  more  sinusoidal  and  has  fewer  harmonics  and  therefore  lower  rotor  losses  [29]. 

The  magnet  poles  are  sometimes  skewed  to  reduce  cogging  torque  and  smooth  out 
variations  in  air  gap  reluctance,  flux,  and  voltage  waveforms.  Skewing  of  the  magnets  occurs 
axially  along  the  length  of  the  rotor  to  provide  a  constant  rotational  torque  and  prevent  pole 
pieces  from  exactly  lining  up  with  stator  teeth.  A  skew  factor  is  used  to  account  for  this  effect 
and  is  shown  in  Eqn  3-3. 


hn  = 


where  es  =  skew  angle,  radE 
n  =  harmonic  number 


Eqn  3-3 

As  the  pole  number  is  increased,  the  stator  conductors-per-pole  decreases  so  that  the  per- 
unit  inductance  and  synchronous  reactance  decreases  with  higher  pole  number.  This  can 
sometimes  result  in  improved  performance  of  the  machine  since  the  reactance  is  lower.  Overall, 
the  initial  16MW  generator  has  6  poles  but  this  is  examined  later  to  determine  an  optimal  design. 

3.2.4  Magnetic  Dimensions 

The  primary  magnetic  dimensions  that  affect  a  PM  machine  are  the  air  gap  and  the 
magnet  height.  These  two  parameters  play  a  major  role  in  determining  the  air  gap  magnetic 
field,  the  air  gap  flux  density,  and  the  induced  voltage  in  the  machine.  To  a  first  order 
approximation,  the  air-gap  flux  density  (Bg)  can  be  represented  by  Eqn  3-4  [30], 

The  radial  air  gap  is  usually  made  as  small  as  possible  to  maximize  the  air  gap  flux 
density,  minimize  the  flux  leakage,  and  produce  a  lower  reluctance  value  since  the  air  gap 
constitutes  the  largest  part  of  the  machine  permeance/reluctance.  However,  the  use  of  rare-earth 
permanent  magnets  (NdFeB  or  SmCo)  with  their  higher  flux  density  and  coercive  force  permit 
some  flexibility  in  the  size  of  the  air  gap. 
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Bg  =  - - Bf 

8  hm  +  8 

where  hm  =  magnet  height  (mm) 
g  =  air  gap  (mm) 

Br  =  magnet  remnant  flux  density  (T) 

I 

Eqn  3-4 

Once  the  permanent  magnet  material  is  selected,  the  desired  air  gap  flux  density  and 
induced  voltage  help  determine  the  magnet  height  needed.  If  the  magnet  height  is  too  large,  the 
air  gap  flux  density  might  be  significant  enough  to  cause  the  stator  core  material  to  saturate 
which  reduces  machine  performance.  The  goal  is  to  use  the  minimal  amount  of  magnet  material 
to  achieve  the  desired  effect  because  this  reduces  the  size  and  weight  of  the  machine  and 
decreases  the  magnet  material  cost.  Also,  losses  in  the  magnets  can  be  reduced  by  using  smaller 
magnets.  In  order  to  provide  uniform  magnetic  fields,  the  magnet  height  is  usually  larger  than 
the  air  gap  by  a  factor  of  5  -  10. 

3.2.5  Number  of  Phases 

In  general,  the  number  of  phases  affects  a  machine’s  power,  current,  and  voltage  ratings 
as  shown  in  Eqn  3-5.  If  the  power  is  fixed,  then  as  the  number  of  phases  increases,  the  phase 
voltage  and/or  current  decreases,  assuming  the  total  number  of  turns  is  constant. 

|P+  jQl  =qVI 

where  P  =  real  power  (W) 

Q  =  reactive  power  (VAR) 
q  =  number  of  phases 
V  =  RMS  phase  voltage  (V) 

I  =  RMS  current  (A) 

Eqn  3-5 

Most  motors  and  generators  are  three-phase  machines  because  it  is  the  industry  standard,  it  is  the 
most  common  form  of  power,  and  it  is  the  lowest  number  of  phases  that  produces  balanced 
torque  with  out  pulsations  in  rotating  machines.  However,  higher  utilizations  in  generators  can 
be  achieved  with  higher  phase  numbers  especially  if  the  generator  is  connected  through  power 
electronics  to  a  DC  bus  distribution.  This  is  because  the  higher  number  of  phases  produces 
lower  ripple  in  the  DC  bus  voltage. 
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However,  the  AC  line  current  harmonics  are  more  substantial  in  increased  phase 
machines  because  the  triple-n  harmonics  are  higher  order  as  the  phase  number  increases.  For 
example,  a  3-phase  machine  suppresses  harmonics  of  order  3n,  a  5-phase  machine  eliminates 
order  5n,  and  a  7-phase  machine  removes  order  7n.  Therefore,  in  higher  phase  machines,  a 
greater  number  of  large  harmonics  result  in  the  AC  line  current.  Also,  as  the  number  of  phases 
increases,  the  phase  inductances  and  reactances  change  since  there  are  a  greater  number  of 
windings  influencing  each  other. 

Most  machines  are  usually  designed  with  the  phases  balanced  meaning  that  they  have  an 
evenly-spaced  phase  distribution  around  the  stator  of  the  machine.  This  produces  voltage 
waveforms  that  are  identical  in  shape  from  phase  to  phase  but  differ  by  a  phase  offset  angle.  In 
order  to  initially  size  the  PM  generator,  it  is  assumed  to  have  three  phases  but  this  will  be 
optimized  later  in  conjunction  with  the  power  electronics  module.  It  is  also  assumed  that  the 
phases  are  always  balanced. 

3.2.6  Slots  per  Pole  per  Phase 

The  number  of  slots  per  pole  per  phase  (m)  is  an  extremely  important  design  parameter 
when  considering  generator  design  and  it  is  calculated  using  Eqn  3-6.  It  is  used  to  help 
determine  the  relationship  and  interactions  between  the  rotor  poles  and  the  stator  windings  as 
well  as  shape  the  generated  back  voltage  of  the  machine.  When  m  is  an  integer,  the  machine  is 
an  integral  slot  machine  and  when  m  has  a  fractional  part,  it  is  a  fractional  slot  machine. 


2pq 

where  Ns  =  number  of  slots 

p  =  pole  pairs 
q  =  number  of  phases 

Eqn  3-6 

In  an  integral  slot  machine,  the  back  EMFs  of  all  of  the  coils  making  up  a  phase  winding 
are  in  phase  with  each  other  and  add  up  so  that  the  final  voltage  amplitude  is  the  direct  sum  of 
the  individual  coil  voltages.  In  a  fractional  slot  machine,  the  back  EMF  of  all  of  the  coils  are  not 
in  phase  so  the  net  voltage  has  a  different  shape  than  the  individual  winding  voltages.  Varying 
the  number  of  slots/pole/phase  is  one  method  used  to  produce  a  more  sinusoidal  voltage 
waveform  and  reduce  the  harmonics  generated  by  the  machine. 
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3.2.7  Stator  Windings 


The  stator  windings  are  the  location  where  the  generator  voltage  is  induced  due  to  the 
time  varying  magnetic  flux  caused  by  the  permanent  magnets  on  the  rotor.  In  a  slotted  machine, 
the  winding  arrangement  is  used  to  help  shape  the  back  voltage  to  produce  a  more  sinusoidal 
waveform.  The  windings  can  be  distributed  by  three  methods:  pitch,  skew,  or 
breadth/distribution. 

The  pitch  of  a  winding  (a)  refers  to  the  angular  displacement  between  the  sides  of  a  coil, 
usually  expressed  in  electrical  degrees  or  radians.  When  the  individual  coil  pitch  differs  from 
180°  E,  the  winding  is  said  to  be  short-pitched  or  fractional-pitched.  This  causes  angular 
segments  where  the  back  voltage  is  zero  because  the  flux  linkage  is  constant  and  can  help 
produce  a  sinusoidal  waveform  when  multiple  coils  are  connected.  It  also  has  the  advantage  of 
lowering  the  coil  resistance  and  making  the  stator  end  windings  more  manageable. 

Windings  in  the  stator  can  also  be  skewed  axially  along  the  length  of  the  machine.  This 
requires  the  stator  slots  to  be  more  intricately  designed  which  complicates  the  mechanical 
construction  of  large  machines.  Therefore,  since  the  generator  being  designed  is  a  large,  high- 
power  machine,  skewing  of  the  stator  windings  is  not  used  but  skewing  of  the  rotor  is  employed. 

The  breadth  of  a  stator  winding  results  from  the  coils  occupying  a  distribution  or  range  of 
slots  within  a  phase  belt.  A  stator  winding  normally  consists  of  several  coils  each  separated  by 
an  electrical  angle  y.  The  distribution  of  the  coils  causes  each  to  link  the  rotor  flux  slightly  out  of 
phase  with  each  other  so  when  they  are  added  together,  they  produce  a  more  sinusoidal 
waveform. 

Within  each  stator  slot,  there  are  geometric  size  constraints  which  determine  how  many 
conductors  can  be  placed  in  a  slot.  In  smaller  machines,  coils  are  composed  of  round  insulated 
wires  that  are  placed  in  the  stator  slot  along  with  insulation  material.  A  slot  fill  factor  (A*)  is  used 
to  determine  how  much  of  the  slot  cross-sectional  area  is  occupied  by  winding  material  as  shown 
in  Eqn  3-7. 

^  _  WindingArea 
s  TotalSlotArea 


Eqn  3-7 

In  larger  machines,  form-wound  windings  are  used  for  ease  of  construction  and  for  better 
performance.  A  sketch  of  what  a  form-wound  winding  looks  like  is  shown  in  Figure  14. 
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Typically,  machines  contain  two  coil  sides  per  slot  making  the  winding  a  double-layer  design 
[17],  [18],  [19].  Overall,  slot  fill  factors  vary  in  value  from  0.30  -  0.70,  depending  on  the 
number  and  size  of  the  conductors  in  the  slots  as  well  as  the  amount  of  labor  utilized.  In  this 
paper,  a  slot  fill  factor  of  0.50  is  assumed. 
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Figure  14:  Example  of  Form-Wound  Winding 

In  conductors  that  carry  high-frequency  currents,  skin  effect  can  become  an  issue  and 
affect  the  operation  of  the  machine.  Skin  effect  is  caused  by  eddy  currents  in  the  windings 
themselves  due  to  the  changing  magnetic  field.  These  eddy  currents  force  the  current  flowing  in 
the  conductor  to  crowd  to  the  outer  edges  of  the  conductor.  This  in  turn  causes  the  current  to 
flow  through  a  smaller  cross-sectional  area  and  increase  the  resistance  of  the  conductor. 
However,  the  generator  under  design  is  expected  to  operate  at  less  than  2  kHz  and  for 
frequencies  below  12  kHz,  Rac/Rdc  <  1.01  so  skin  effect  can  be  neglected  [32]. 

Within  a  phase,  stator  windings  can  be  connected  in  wye  or  delta  patterns  as  well  as 
series  or  parallel.  Almost  all  machines  use  series,  wye-connected  windings  because  they  provide 
the  safest  alternative.  This  is  because  in  a  delta  or  parallel  connection,  the  back  EMFs  can 
produce  circulating  currents  which  can  result  in  addition  losses,  heating,  or  damage.  Therefore, 
wye  series  connected  windings  are  selected  for  use  in  the  designs  in  this  paper. 
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3.3  Machine  Calculated  Parameters 

3.3.1  Basic  Model 

Since  the  machine  is  assumed  to  balanced,  parameters  can  be  determined  on  a  per-phase 
basis  and  then  applied  to  all  of  the  phases.  Each  phase  of  the  machine  can  therefore  be  modeled 
as  shown  in  Figure  15. 


Figure  15:  Per  Phase  Model 

The  armature  resistance  (Ra)  is  the  resistance  of  the  windings  of  the  machine  and  it  is  usually 
relatively  small.  The  synchronous  inductance  (Ls)  of  the  machine  comes  from  the  inductance  of 
the  windings  and  is  composed  of  the  air  gap  inductance,  the  slot  leakage  inductance,  and  the  end- 
turn  inductance.  The  back  voltage  (Ea)  is  produced  through  the  flux  linkage  in  the  windings 
from  the  rotating  magnetic  field  in  the  machine.  Lastly,  Va  is  the  terminal  voltage  and  is  found 
using  basic  circuit  analysis  once  the  other  parameters  are  known. 

3.3.2  Winding  Resistances 

The  stator  coils  in  the  machine  are  made  of  copper  and  therefore  have  some  resistance  to 
the  current  flow.  This  resistance  of  the  copper  phase  windings  is  calculated  using  Eqn  3-8. 


where  I  =  length  of  conductor 
a  =  winding  conductivity 
A  =  winding  cross-sectional  area 

Eqn  3-8 

The  length  of  the  conductor  comes  from  the  windings  traveling  twice  the  length  of  the  machine 
and  twice  around  the  end  turns  of  the  machine.  It  is  assumed  that  the  end  turns  follow  roughly  a 
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circular  path  from  one  side  of  the  machine  to  the  other  where  the  radius  of  the  circle  is  the 
distance  to  one  half  the  stator  slot  height.  The  cross-sectional  area  of  the  conductor  is  obtained 
from  the  slot  area  and  slot  fill  factor  as  shown  in  Eqn  3-9,  assuming  form-wound  windings. 

A.'L 
A  =  ■■■ 
ac  2-Nc 

where  As  =  slot  area 

Nc  =  turns  per  coil 


Eqn  3-9 

3.3.3  Winding  &  Magnet  Factors 

As  discussed  in  section  3.2.7,  windings  are  normally  not  full-pitched  or  concentrated  but 
rather  are  short-pitched  and  have  breadth  associated  with  them.  To  account  for  these  effects,  a 
winding  factor  (kw)  is  utilized  which  is  the  ratio  of  flux  linked  by  an  actual  winding  to  the  flux 
linked  by  a  full-pitch,  concentrated  winding  having  the  same  number  of  turns.  The  winding 
factor  is  the  product  of  a  pitch  factor  (kp)  and  a  breadth/distribution  factor  (kb)  as  shown  in  Eqn 
3-10. 

1c  —  V  •  1c  i 

Rwn"  K  pn  ^bn 

Eqn  3-10 

The  pitch  factor  accounts  for  the  windings  spanning  a  electrical  degrees  vice  spanning  a 
full  180°  E  as  shown  in  Figure  16  [26]. 


Figure  16:  Short-Pitch  Coil 
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The  pitch  factor  is  the  ratio  of  the  flux  produced  by  a  short-pitch  coil  to  the  flux  produced  by  a 
full-pitch  coil.  Short-pitching  is  an  important  means  for  eliminating  harmonics  and  improving 
the  power  quality  of  the  machine.  The  pitch  factor  can  be  derived  with  the  final  result  shown  in 
Eqn  3-11. 


f  na  ^ 

^  n-n  ' 

=  sin 

•sin 

— 

l  2  J 

1  2  J 

where  n  =  harmonic  number 


Eqn  3-11 

The  breadth  factor  explains  the  effect  of  the  windings  occupying  a  distribution  or  range 
of  slots  within  a  phase  belt.  A  phase  winding  normally  consists  of  numerous  coils  connected 
together  linking  flux  slightly  out  of  phase  with  each  other  as  shown  in  Figure  17  [26], 


The  breadth  factor  can  be  derived  either  magnetically  or  geometrically  to  obtain  Eqn  3-12. 


sin 


n-m-y 


^bn  = 


V  2  J 


m-sin 


^ny 


where  n  =  harmonic  number 

m  =  slots  per  pole  per  phase 
y  =  coil  electrical  angle 


Eqn  3-12 

In  addition  to  estimating  different  winding  effects,  the  geometry  of  the  magnetic  air  gap 
must  be  represented.  Field  methods  are  utilized  along  with  vector  potential  analysis  to  develop 
expressions  that  account  for  different  magnetic  gap  geometries.  Reference  [19]  contains  detailed 
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derivations  of  the  magnetic  gap  factor  (kgn)  for  several  magnet  and  slot  configurations.  The 
equation  for  the  slotted  stator,  surface  magnet  configuration  is  shown  in  Eqn  3-13. 


where  Rs  =  outer  magnetic  boundary  R2  =  outer  boundary  of  magnet 

R|  =  inner  magnetic  boundary  R1  =  inner  boundary  of  magnet 

Eqn  3-13 

3.3.4  Flux  and  V oltage 

The  primary  significance  of  the  magnetic  flux  linkage  in  a  machine  is  that  it  induces 
voltage  across  a  winding  whenever  the  flux  varies  with  time  as  explained  through  Faraday’s 
Law.  The  first  step  in  the  process  is  to  determine  the  air  gap  flux  density.  The  flux  from  the 
magnet  poles  crosses  the  air  gap  to  the  stator  windings  but  some  flux  leaks  along  the  way  and 
this  is  accounted  for  using  a  leakage  factor  (Ki  ~  0.95  for  surface  magnets).  In  addition,  the  flux 
path  is  normally  dominated  by  the  air  gap  reluctance  since  the  reluctance  of  the  stator  steel  is 
much  less  than  that  in  the  air  gap.  However,  a  reluctance  factor  (Kr  ~  1 .05  for  surface  magnets) 
is  used  to  compensate  for  the  small  effects  of  the  steel  reluctance  on  the  air  gap  flux. 

The  presence  of  the  slots  in  the  stator  also  affects  the  air  gap  flux  density  because  of  the 
difference  in  permeance  caused  by  the  slots.  The  flux  crossing  the  air  gap  in  a  slot  region  travels 
farther  before  reaching  the  highly  permeable  stator  back  iron.  Carter’s  coefficient  (Kc)  is  used  to 
account  for  this  effect  [17].  The  air  gap  flux  density  is  also  affected  by  the  magnet  geometry  in 
the  air  gap  as  previously  described  by  Eqn  3-13.  Since  the  magnet  poles  rotate  north/south,  the 
air  gap  flux  density  shape  can  be  approximated  as  shown  in  Figure  18.  This  can  be  represented 
as  a  Fourier  series  using  only  odd  components  because  of  half-wave  symmetry  as  shown  in  Eqn 
3-18.  Overall,  the  air  gap  flux  density  is  calculated  using  Eqn  3-14  through  Eqn  3-18. 
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Figure  18:  Air  Gap  Flux  Density 


where  ws  =  average  slot  width 
wt  =  tooth  width 

*s  =  ws+wt 

Eqn  3-14 
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Eqn  3-15 
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where  PC  =  permeance  coefficient 

=  flux  concentration  factor  (A  m/Ag) 


Eqn  3-16 
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where  (irec  =  recoil  permeability 
Br  =  remnant  flux  density 


Eqn  3-17 
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CO 


B(e)  =  Y,  Bn-sin(np9) 


n  =  1 

n  odd 


where 


np0 
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Bn  =  -'BgVsi\  .  , 

9m  =  magnet  physical  angle 
n  =  harmonic  number 


sin 


'  nn^ 
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Eqn  3-18 

Once  the  flux  density  is  known,  the  flux  must  be  calculated.  Given  a  machine  uniform  in 
the  axial  direction,  the  flux  linked  by  a  single,  full-pitched  coil  which  spans  an  angle  from  0  to 
7i/p  is  represented  by  Eqn  3-19.  Assuming  Bfiux  is  sinusoidally  distributed,  the  peak  flux  for  this 
ideal  coil  is  given  by  Eqn  3-20. 


O  = 
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P 

Bflux  Rs  "Lst  de 

0 


where  B  flux  =  radial  flux  through  coil 


flux 

P 


Eqn  3-19 


Eqn  3-20 

Given  that  there  are  Na  coils  in  a  stator  phase  winding  and  including  all  the  real  winding  effects, 
the  total  flux  linkage  is  shown  in  Eqn  3-21.  Through  Faraday’s  Law,  the  back  EMF  for  the 
machine  is  given  by  Eqn  3-22. 


oo 
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Eqn  3-21 


46 


oo 


Ea=  ^  Vn  sin(np0) 


n  =  1 

n  odd 


where  Vn  =  ^n  =  co0^n 


Eqn  3-22 

With  permanent  magnet  excitation,  the  field  cannot  be  controlled  as  in  a  wound-rotor 
machine  so  the  number  of  turns  in  the  stator  phase  windings  must  be  chosen  so  that  the  machine 
EMF  is  close  to  the  nominal  system  voltage.  Another  option  is  to  use  power  electronics  to 
convert  the  machine  EMF  to  the  system  voltage  thereby  providing  steady  regulation. 

Another  effect  that  is  sometimes  an  issue  is  armature  reaction.  It  is  caused  by  current 
flowing  in  the  stator  windings  which  creates  a  magnetic  field  that  tends  to  distort  the  magnetic 
field  established  by  the  permanent  magnets.  For  surface-magnet  machines,  because  the  magnet 
recoil  permeability  is  approximately  one  and  the  magnet  height  is  large  compared  to  the  air  gap, 
the  armature  reaction  flux  density  is  small.  As  long  as  the  stator  teeth  are  not  highly  saturated 
due  to  the  permanent  magnets,  armature  reaction  is  negligible  [17],  [18]. 


3.3.5  Machine  Inductances 

In  a  slotted  permanent  magnet  machine,  there  are  three  distinct  components  of 
inductance:  air  gap  inductance,  slot  leakage  inductance,  and  end-tum  inductance.  The  most 
accurate  means  for  calculating  these  parameters  is  finite  element  analysis  but  analytical  methods 
provide  almost  as  good  results  (within  a  couple  of  percent)  and  are  used  in  this  paper  [33], 

The  air  gap  inductance  is  usually  the  largest  portion  of  the  total  inductance  and  it  is  due 
to  the  interaction  of  the  stator  windings  with  the  flux  crossing  the  air  gap.  To  calculate  the  air 
gap  inductance,  a  full-pitch,  concentrated  winding  carrying  a  current  I  is  initially  examined 
which  leads  to  an  air  gap  flux  density  shown  in  Eqn  3-23. 


oo 

Bflux=  X  Bn'sin(nP0) 


n  =  1 

n  odd 


where 


B„ 


4  FO  jV 

Ml  (g  +  hm)  2p 


Eqn  3-23 
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When  this  concept  is  expanded  to  polyphase  windings  with  balanced  operation,  the  air  gap  flux 
density  becomes  Eqn  3-24. 


oo 

Bflux=  2  Vsin(nP0) 


n  =  1 
n  odd 


where 
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Eqn  3-24 

The  flux  can  be  found  using  equation  Eqn  3-19  and  the  total  flux  linkage  is  X  =  NaO.  With  all 
real  winding  effects  included,  the  air  gap  inductance  is  then  given  by  Eqn  3-25. 


q  4  P0'Bs'Bst'^a  '^wn 

2"”  nV-(«  +  hJ 


Eqn  3-25 

In  addition  to  the  air  gap,  the  coil  currents  generate  a  magnetic  field  that  crosses  from  one 
side  of  the  slot  to  the  other  producing  a  slot  leakage  inductance.  For  calculating  the  slot  leakage 
inductance,  it  is  assumed  that  the  slot  is  rectangular  with  slot  depressions  (Figure  13)  which 
results  in  a  slot  permeance  per  unit  length  shown  in  Eqn  3-26  [17],  [18],  [19]. 
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Eqn  3-26 

Assuming  m  slots  per  pole  per  phase  and  a  standard  double  layer  winding,  it  can  be  shown  that 
the  slot  leakage  inductance  is  given  by  Eqn  3-27  through  Eqn  3-29  [19]. 


Las  =  2-p-Lst-Perm 
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Eqn  3-27 


Lam  =  2-p- Lgt  •  Perm  Ngp  •  Nc  (mutual) 


Eqn  3-28 
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Eqn  3-29 

The  end  turn  inductance  is  the  smallest  of  the  three  components.  It  is  created  by  the 
magnetic  field  that  surrounds  a  coil  after  it  leaves  one  slot  and  before  it  enters  another  slot. 

Since  it  is  extremely  difficult  to  accurately  determine  because  of  complex  winding  patterns,  a 
rough  approximation  is  used.  It  is  assumed  that  the  end  turns  are  semi-circular  with  a  radius 
equal  to  one-half  the  mean  coil  pitch.  Using  reference  [17],  the  total  end  turn  inductance  per 
phase  is  shown  in  Eqn  3-30. 
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Eqn  3-30 

The  total  inductance  for  the  phase  is  the  sum  of  the  three  inductances,  ignoring  other  small 
factors.. 


^ag  +  ^slot  +  ^e  Xs  w0'^s 


Eqn  3-31 

\ 

3.3.6  Basic  Losses 

Losses  in  a  machine  consist  of  core  losses,  conductor  losses,  friction  and  windage  losses, 
and  rotor  losses.  Rotor  losses  include  magnet  losses  and  retaining  can  losses,  they  require 
detailed  waveform  and  harmonic  analysis,  and  therefore  are  discussed  later  in  Chapter  5. 

3.3.6. 1  Core  Losses 

High-speed  generator  stator  core  losses  (per  weight)  can  be  greater  than  normal  machines 
because  of  the  higher  frequencies.  These  losses  are  minimized  by  using  laminated  steels  in  the 
stator  construction  as  discussed  in  3.1.2  and  by  not  generating  frequencies  that  are  too  high. 

Core  losses  consist  of  hysteresis  and  eddy  current  losses.  Hysteresis  loss  results  from  the 
steel  not  wanting  to  change  magnetic  state.  As  the  flux  density  varies,  the  material  traverses  the 
B-H  curve  and  energy  is  lost.  Eddy  current  loss  is  also  caused  by  the  variation  in  flux  density. 
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Electrical  currents  are  induced  within  the  ferromagnetic  material  and  circulate  dissipating  power 
because  of  the  resistivity  of  the  material. 

Because  there  are  usually  various  imperfections  in  materials,  the  best  way  to  approximate 
core  losses  is  to  use  empirical  loss  data.  If  the  flux  density  is  estimated  for  each  part  of  a 
machine  and  the  mass  of  the  steel  calculated,  empirical  core  loss  data  can  be  used  to  estimate  the 
total  losses.  Empirical  data  for  M-19,  29  gauge  material  is  obtained  as  shown  in  Figure  19. 
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Figure  19:  Core  Loss  Data 

An  exponential  curve  fit  is  then  applied  to  the  data  to  obtain  an  equation  for  estimating  the  core 
losses  (Eqn  3-32)  [9],  [35],  [36], 


pc  =  p0- 


The  values  for  the  bases  and  exponents  in  Eqn  3-32  are  listed  in  Table  8. 


Eqn  3-32 
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Table  8:  Core  Loss  Parameters 


llili  Parameter 

Value 

Base  Power  (Po) 

36.79  W/lb 

Base  Flux  Density  (Bo) 

LOT 

Flux  Density  Exponent  (sb) 

2.12 

Base  Frequency  (fo) 

1000  Hz 

Frequency  Exponent  (sf) 

1.68 

3. 3. 6.2  Conductor  Losses 


Conductor  losses  arise  from  the  current  flowing  in  the  stator  windings.  The  resistance  of 
the  windings  is  calculated  using  Eqn  3-8  from  section  3.3.2.  The  conductor  losses  are  then  found 
using  the  traditional  power  equation  for  a  resistance  (Eqn  3-33). 


Pa  =  q*I„ 


Eqn  3-33 

3.3.6. 3  Friction  &  Windage  Losses 

For  rotors  operating  at  high-speed,  friction  and  windage  in  air  can  cause  losses  which 
result  in  inefficiency  and  heat  production.  These  losses  are  calculated  using  the  power  necessary 
to  overcome  the  drag  resistance  of  a  rotating  cylinder  as  given  by  Eqn  3-34  [32]. 


3  4 

*\vind  -  Cf7t-p  air®  ^  '^st 

where  Cf  =  friction  coefficient 
pair  =  density  of  air 


Eqn  3-34 

The  friction  coefficient  depends  on  numerous  factors  such  as  surface  roughness  and  flow  region. 
Since  the  air  gap  is  a  small  annulus  and  the  rotor  is  spinning  at  high  speed,  it  is  assumed  that  the 
air  in  the  gap  is  in  the  turbulent  region.  Therefore,  the  coefficient  of  friction  can  be 
approximated  by  Eqn  3-35  [34]. 

Cf  =  0.0725  -Rey^20 
where  Rey  =  Reynold's  Number 


Eqn  3-35 
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3.4  Machine  Sizing  Methods 


3.4.1  Basic  Sizing  Method 

Whenever  a  machine  is  being  designed,  it  is  important  to  perform  some  back-of-the- 
envelope  calculations  to  gain  insight  into  initial  sizing  estimates.  Most  generators  are 
constrained  by  two  competing  design  parameters,  air  gap  magnetic  shear  stress  and  rotor  tip 
speed.  Air  gap  magnetic  shear  stress  (t)  is  the  magnetic  shear  force  developed  per  unit  gap  area 
and  is  constrained  by  magnetic  design  and  thermal  management  [8],  It  is  proportional  to  the 
product  of  the  surface  current  density  and  magnetic  flux  density  as  shown  in  Eqn  3-36. 

X  Kz-  Bg 

where  x  =  shear  stress  (psi) 

Kz  =  surface  current  density 
Bg  =  air  gap  flux  density 


Eqn  3-36 

Typical  values  for  air  gap  shear  stress  for  different  types  of  generators  are  shown  in  Table  9  [8], 

[1 1],  [16],  [18],  [31].  For  the  basic  sizing  calculations,  15  psi  is  assumed  since  the  generator  is  a 
large  liquid-cooled  machine. 


Table  9:  Air  Gap  Shear  Stress  Values 


Small  Air-Cooled 

1  -  5  psi 

Large  Air-Cooled 

HHHEHBfJSSHHI 

Large  Liquid-Cooled 

High-Temperature  Superconducting 

Rotor  tip  speed  is  discussed  previously  in  section  3.2.2  and  given  by  Eqn  3-1  with  the 
assumed  limit  of  200  m/s.  The  fundamental  machine  power  equation  is  utilized  to  derive  the 
rotor  radius  and  stack  length  of  the  machine  (Eqn  3-37). 

P=2-7tRLst-Tvtip 

where  R  =  rotor  radius 
Lst  =  stack  length 


Eqn  3-37 
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In  order  to  simplify  the  equation,  the  L/D  ratio  is  substituted  in  for  L*.  Using  the  air  gap  shear 
stress,  the  rotor  tip  speed  limit,  and  the  power  rating  of  the  machine  (16  MW),  the  basic  power 
equation  is  iterated  to  obtain  a  final  L/D  ratio,  rotor  radius,  and  stack  length  while  matching  the 
desired  rotational  speed  of  the  machine  (13,000  RPM).  Using  an  air  gap  flux  density  of  0.8  T,  a 
pole  pair  value  of  3,  a  slot  height  of  15  mm,  and  a  slot  fill  fraction  of  0.5,  the  frequency  and 
current  density  of  the  machine  are  found.  The  detailed  MATLAB  code  is  contained  in  Appendix 
B  with  the  results  shown  below. 


Basic  Machine  Design 

Input  Parameters: 

Power  =  16000.0  kW 

Shear  Stress  = 

15.0  psi 

L/D  Ratio  = 

2.85 

Tip  Speed  = 

200.0  m/s 

Pole  Pairs  = 

3.0 

Air  Gap  Bg  = 

0.80  T 

Output: 

Rotor  Radius  = 

0.147  m 

Stack  Length  = 

0.838  m 

Speed  = 

13000  RPM 

Frequency  = 

650.0  Hz 

Ja 

1757.67  A/cm2 

3.4.2  Detailed  Sizing  Method  One 

Once  basic  sizing  of  the  machine  is  complete,  in-depth  analysis  is  conducted  to  ascertain 
the  overall  performance  and  scale  of  the  16  MW  generator.  Two  detailed  sizing  methods  are 
developed  using  MATLAB  code.  To  help  gain  an  understanding  for  current  PM  machine 
designs,  numerous  articles,  transactions,  and  proceedings  are  examined  and  a  database  is 
developed  containing  comprehensive  specifications  for  over  20  different  PM  machines  (see 
Appendix  C).  This  database  is  used  to  help  develop  input  parameters  for  the  first  method  (Table 
10),  although  many  of  these  values  are  obtained  after  numerous  iterations  through  the  code. 


Table  10:  Input  Parameters  for  Sizing  Method  1 


1  Parameter 

SliSllalue 

lililtaliliiiir 

Required  Power 

16  MW 

Number  of  Phases 

3 

Rotational  Speed 

13,000  RPM 

Number  of  Slots 

36 

Power  Factor  Angle 

0  deg 

Slots  Short-pitched 

1 

Rotor  Radius 

0.147  m 

Peripheral  Tooth  Fraction 

6.5 

Magnet  Height 

25  mm 

Slot  Depth 

25  mm 

Stack  Length 

0.838  m 

Slot  Depression  Depth 

0.5  mm 

Pole  Pairs  (p) 

3 

Slot  Depression  Width 

N/A 

Magnet  Br 

1.2  T 

Stator  Back  Iron  Ratio 

0.7 

Magnet  Angle 

Turns  per  Coil 

1 

Magnet  Skew  Angle 

Slot  Fill  Fraction 

0.5 

53 


Air  Gap _ 

Steel  Density 
Magnet  Density 


4  mm 
7700  kg/m 
7400  kg/m 


Winding  Conductivity 

Conductor  Density 


6.0  x  10v  S/m 
8900  kg/m3 


Once  the  input  parameters  are  entered,  the  first  step  in  sizing  is  to  generate  the  geometry 


of  the  machine.  This  includes  determining  the  number  of  slots  per  pole  per  phase  (Eqn  3-5), 
number  of  armature  turns  (Eqn  3-38),  tooth  width,  slot  dimensions,  stator  back  iron  dimensions, 
coil  pitch,  and  winding  end  turn  geometry. 


Na  =  2  pm- Nc 

where  Nc  =  Turns  per  coil 

Na  assumes  each  slot  has  2  half  coils 


Eqn  3-38 

Next,  the  electrical  frequency  and  rotor  surface  speed  are  determined  using  Eqn  3- land  Eqn 
3-39. 


-  P'N 
60 

where 


co  =  2-Tif  co  =  pcom 

co  =  electrical  frequency  (rad/sec) 
com  =  mechanical  frequency  (rad/sec) 


Eqn  3-39 

Winding,  skew,  and  magnetic  gap  factors  are  then  estimated  as  discussed  in  sections  3.2.3  and 
3.3.3.  The  magnitude  of  the  air  gap  magnetic  flux  density  (Bg)  is  determined  accounting  for 
slots,  varying  reluctances,  and  flux  leakage  per  Eqn  3-14  through  Eqn  3-17. 

The  magnetic  flux  and  back  voltage  magnitudes  can  then  be  calculated.  The  fundamental 
component  of  the  magnetic  flux  links  the  stator  windings  to  create  useful  voltage.  Therefore, 
only  the  fundamental  components  of  Eqn  3-21  and  Eqn  3-22  are  used  to  determine  the  internal 
voltage  of  the  generator  as  shown  in  Eqn  3-40  and  Eqn  3-41. 


B  l  = 
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Pe  m 

V  2  j 


Eqn  3-40 
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Ea  -  o)Q  X 


p 


Eqn  3-41 

Utilizing  equations  from  section  3.3.5,  the  machine  inductances  and  reactances  are 
established.  The  lengths,  volumes,  and  masses  of  components  and  the  overall  generator  are 
calculated  using  basic  geometric  equations  and  relationships.  A  15%  service  mass  fraction  is 
added  to  the  total  mass  estimate  to  account  for  the  additional  services  associated  with  large 
liquid-cooled  machines  [31].  Once  the  mass  of  each  of  the  stator  parts  is  known,  the  core  losses 
are  estimated  in  accordance  with  section  3.3.6. 1. 

The  terminal  voltage  and  current  of  the  machine  must  then  be  calculated  accounting  for 
conductor  losses  and  windage  losses.  The  vector  relationship  (Figure  20)  between  terminal 
voltage  (Va),  internal  voltage  (Ea),  and  the  synchronous  reactance  voltage  drop  is  utilized  to 
obtain  Eqn  3-42.  The  armature  resistance  is  usually  ignored  because  it  is  much  smaller  than  the 
synchronous  reactance. 


Figure  20:  Voltage  Vector  Relationship 


va  =  JEa2  “  (x s  1  a-cosv)2  -  X g-I a-siny 

Eqn  3-42 
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The  machine  efficiency  (r|)  is  then  easily  obtained.  Appendix  D  contains  the  detailed  MATLAB 

code  and  the  results  are  shown  below. 

PM  Machine  Design,  Version  1:  Surface  Magnet,  Slotted  Stator 
Machine  Size: 


Machine  Diameter  = 

0.472  m 

Machine  Length  = 

1.003  m 

Rotor  radius  = 

0.147  m 

Active  length  = 

0.838  m 

Slot  Avg  Width  = 

16.493  mm 

Slot  Height  = 

25.000  mm 

Back  Iron  Thick  = 
Machine  Ratings: 

34.300  mm 

Tooth  Width  = 

15.403  mm 

Power  Rating  = 

16000.0  kW 

Speed  = 

13000  RPM 

Va  (RMS)  = 

2341  V 

Current  = 

2267.5  A 

Ea  (RMS)  = 

2925  V 

Arm  Resistance  = 

0.00526  ohm 

Synch  Reactance  = 

0.768  ohm 

Synch  Induct  = 

0.188  mH 

Stator  Cur  Den  = 

2199.7  A/cm2  Tip  Speed  = 

200  m/s 

Efficiency  = 

0.992 

Power  Factor  = 

1.000 

Phases  = 

Stator  Parameters: 

3 

Frequency  = 

650.0  Hz 

Number  of  Slots  = 

36 

Num  Arm  Turns  = 

12 

Breadth  Factor  = 

0.966 

Pitch  Factor  = 

0.966 

Tooth  Flux  Den  = 
Slots/pole/phase  = 

Rotor  Parameters: 

1.59  T 

2.00 

Back  Iron  = 

1.14  T 

Magnet  Height  = 

25.00  mm 

Magnet  Angle  = 

50.0  degm 

Air  gap  = 

4.00  mm 

Pole  Pairs  = 

3 

Magnet  Remanence  = 

1.20  T 

Aig  Gap  Bg  = 

0.80  T 

Magnet  Factor  = 

Machine  Losses: 

0.949 

Skew  Factor  = 

0.989 

Core  Loss  = 

11.6  kW 

Armature  Loss  = 

81.2  kW 

Windage  Loss  = 

Machine  Weights: 

30.7  kW 

Rotor  Loss  = 

TBDkW 

Core  = 

396.49  kg 

Shaft  = 

438.05  kg 

Magnet  = 

64.74  kg 

Armature  = 

89.58  kg 

Services  = 

148.33  kg 

Total  = 

1137.18  kg 

3.4.3  Detailed  Sizing  Method  Two 

In  order  to  provide  a  check  on  the  methodology  of  the  previous  sizing  procedure,  a 
second  MATLAB  code  is  constructed.  The  second  method  is  developed  using  a  combination  of 
processes  from  references  [2],  [17],  [18],  [19],  [26],  and  [37]. 
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Table  11:  Input  Parameters  for  Sizing  Method  2 


I  Parameter 

Value 

16  MW 

Current  Density  (Ja) 

Rotational  Speed 

13,000  RPM 

Number  of  Phases 

3 

Power  Factor  Angle 

0  deg 

Slots/Pole/Phase  (m) 

2 

Maximum  VtiP 

200  m/s 

Slots  Short-pitched 

1 

Bsat 

1.65  T 

Avg  Slot  Width  (ws) 

16.0  mm 

Pole  Pairs  (p) 

3 

Slot  Depth 

25  mm 

Magnet  Br 

1.2  T 

Slot  Depression  Depth 

0.5  mm 

Permeance  Coefficient  (PC) 

5.74 

Slot  Depression  Width 

N/A 

Magnet  Skew  Angle 

Turns  per  Coil 

1 

Air  Gap 

4  mm 

Slot  Fill  Fraction 

0.5 

■ 

WSmMSM, 

Winding  Conductivity 

■BENESH 

Conductor  Density 

wr-'rrrw 

Many  of  the  same  equations  used  in  the  first  method  are  utilized  but  the  process  differs  in  some 
of  the  input  parameters  and  calculations.  The  second  sizing  method  inputs  are  shown  in  Table 
1 1  with  different  entries  from  method  1  being  vtiP,  Bsat,  ws,  m,  Ja,  and  PC. 

First,  given  the  maximum  vtjP  (200  m/s)  and  the  rotational  speed  (13,000  RPM),  the 
electrical  frequency  and  rotor  radius  are  computed.  Next,  the  winding  and  skew  factors  are 
determined  similar  to  the  first  sizing  method.  The  magnet  dimensions  (hm  and  0m)  are  then 
determined  along  with  the  air  gap  flux  density.  Using  Bsat,  the  PC,  and  Eqn  3-43  through  Eqn 
3-46,  the  magnet  arc  width  and  height  are  iterated  until  the  tooth  width  (wt)  equals  the  average 
slot  width  (ws).  This  differs  from  the  first  method  where  the  peripheral  tooth  fraction  is  used  to 
set  the  tooth  width  equal  to  the  slot  top  width  (wst). 


where  am  =  magnet  pitch  coverage  coefficient 

Eqn  3-43 

^  m=  &  e 

Eqn  3-44 
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Eqn  3-46 

The  permeance  coefficient  (PC)  normally  varies  between  5-15  with  higher  PCs  only 
used  for  high  performance  space  and  aircraft  applications  [18],  [37].  A  PC  value  of  4  -  6  is 
typical  for  large  PM  machines  where  the  air  gap  flux  density  is  approximately  60  -  80%  of  the 
remnant  flux  density  [17].  The  PC  is  input  as  5.74  to  coincide  with  the  PC  calculated  in  the  first 
sizing  method. 

Once  the  magnet  dimensions  are  known,  the  geometry  of  the  machine  is  generated  using 
similar  equations  as  the  first  method  except  for  the  stator  core  back  iron  depth.  In  the  first 
method,  the  input  variable  “stator  back  iron  ratio”  is  used  whereas  in  the  second  method,  the 
saturation  flux  density  is  applied  as  shown  in  Eqn  3-47. 

7i  D  g-6  m  B  g 

C=  4'P  ®  sat 


Eqn  3-47 

Utilizing  a  power  balancing  procedure,  the  stack  length,  terminal  voltage,  and  terminal 
current  are  iterated  to  obtain  a  complete  design.  Core  losses,  conductor  losses,  and  windage 
losses  are  calculated  in  accordance  with  sections  3.3.6. 1,  3. 3.6.2,  and  3. 3.6. 3.  The  power 
crossing  the  air  gap  of  the  machine  is  determined  using  Eqn  3-48  [2],  [37]. 

^gap  =  4-Jt-ke-k"kw-kg-ks-sin(0m^-— -(Kz-Bg)-Ds  ‘^st 

where  ke  =  electrical  power  waveform  factor 
kj  =  current  waveform  factor 


Eqn  3-48 
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The  electrical  power  waveform  factor  is  defined  in  Eqn  3-49  and  the  current  waveform  factor  is 
shown  in  Eqn  3-50.  These  factors  depend  on  the  back  EMF  and  current  waveforms  generated  in 
the  machine  and  a  detailed  discussion  of  these  factors  is  in  reference  [2]. 


k 


e 


_j_ 

T 


EpkIpk 


ki  = 


pk 


rms 


Eqn  3-49 


Eqn  3-50 

During  the  balancing  procedure,  the  machine  inductances  and  reactances  are  established 

per  section  3.3.5.  The  lengths,  volumes,  and  masses  of  components  and  the  overall  generator  are 

calculated  using  basic  geometric  equations  and  relationships  similar  to  the  first  sizing  method. 

Once  the  balancing  procedure  is  complete,  the  efficiency  of  the  machine  is  determined. 

Appendix  E  contains  the  detailed  MATLAB  code  and  the  results  are  shown  below. 

PM  Machine  Design,  Version  2:  Surface  Magnet,  Slotted  Stator 
Machine  Size: 


Machine  Diameter  = 

0.484  m 

Machine  Length  = 

0.979  m 

Rotor  radius  = 

0.147  m 

Active  length  = 

0.813  m 

Slot  Avg  Width  = 

16.000  mm 

Slot  Height  = 

25.000  mm 

Back  Iron  Thick  = 

39.587  mm 

Tooth  Width  = 

15.863  mm 

Machine  Ratings: 
Power  Rating  = 

16000.0  kW 

Speed  = 

13000  RPM 

Va  (RMS)  = 

2604  V 

Current  = 

2064.1  A 

Ea  (RMS)  = 

3007  V 

Arm  Resistance  = 

0.00534  ohm 

Synch  Reactance  = 

0.723  ohm 

Synch  Induct  = 

0.177  mH 

Stator  Cur  Den  = 

2064.1  A/cm2 

Tip  Speed  = 

200  m/s 

Efficiency  = 

0.993 

Power  Factor  = 

1.000 

Phases  = 

3 

Frequency  = 

650.0  Hz 

Stator  Parameters: 
Number  of  Slots  = 

36 

Num  Arm  Turns  = 

12 

Breadth  Factor  = 

0.966 

Pitch  Factor  = 

0.966 

Tooth  Flux  Den  = 

1.70  T 

Back  Iron  = 

1.05  T 

Slots/pole/phase  = 
Rotor  Parameters: 

2.00 

Magnet  Height  = 

26.17  mm 

Magnet  Angle  = 

47.7  degm 

Air  gap  = 

4.00  mm 

Pole  Pairs  = 

3 

Magnet  Remanence  = 

1.20  T 

Aig  Gap  Bg  = 

0.85  T 
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Magnet  Factor  = 
Machine  Losses: 

0.958 

Skew  Factor  = 

0.989 

Core  Loss  = 

11.9  kW 

Armature  Loss  = 

77.5  kW 

Windage  Loss  = 
Machine  Weights: 

29.7  kW 

Rotor  Loss  = 

TBDkW 

Core  = 

438.62  kg 

Shaft  = 

424.57  kg 

Magnet  = 

62.89  kg 

Armature  = 

85.50  kg 

Services  = 

151.74  kg 

Total  = 

1163.32  kg 

3.4.4  Comparison  of  Methods 

Both  sizing  methods  produce  similar  sized  generators  with  the  results  agreeing  within  1  - 
5%  on  most  parameters.  The  main  reasons  for  the  differences  are  twofold.  First,  the  calculations 
for  the  magnet  dimensions  are  different.  In  the  first  method  these  parameters  are  input  whereas 
in  the  second  they  are  determined  through  an  iterative  process  using  material  properties.  Second, 
the  procedure  for  determining  the  slot  and  tooth  widths  differs  slightly. 

Overall,  both  methods  underestimate  the  overall  dimensions  (length  and  diameter)  and 
weight.  This  is  because  the  sizing  programs  do  not  include  calculations  for  portions  such  as 
structure,  frames,  mounts,  and  maintenance  access.  These  pieces  are  added  in  section  6.1  to  get 
the  total  sizes  and  weights  of  the  PM  generators,  specifically  including: 

•  50%  weight  factor  added  to  total  weight 

•  0.8  m  added  to  overall  length  and  0.4  m  added  to  overall  diameter 

In  later  chapters  when  waveforms  are  examined  and  machine  optimization  occurs,  the  first  sizing 
method  is  mostly  employed  with  the  second  method  used  for  verification. 
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Chapter  4  Power  Electronics  and  Conversion 


4.1  Background 

In  conjunction  with  employing  a  PM  generator,  a  DC  bus  architecture  is  one  of  the 
preferred  schemes  for  the  future  [5].  DC  power  distribution  systems  can  offer  an  advantage  in 
size  and  a  weight  over  high-power  AC  systems  [3].  Therefore,  a  power  electronics  module  is 
required  to  rectify  the  AC  generator  output  and  then  convert  it  to  the  appropriate  DC  distribution 
voltage  of  700  VDC.  However,  conversion  of  the  high-frequency  generator  AC  output  to  DC 
generates  harmonic  voltages  and  currents  on  the  AC  side  and  ripple  on  the  DC  side.  These 
harmonics  can  be  reduced  through  proper  design  of  the  conversion  module. 

Since  the  generator  is  a  permanent  magnet  machine,  the  field  cannot  be  controlled  like  in 
a  wound-rotor  machine.  The  power  electronics  module  (PEM)  is  needed  not  only  to  convert  the 
machine  EMF  to  the  DC  system  voltage  but  also  provide  steady  regulation.  Also,  since  the 
frequency  of  the  generator  output  is  higher  than  the  normal  60  Hz  systems,  the  filter 
requirements  for  the  PEM  are  less  since  the  components  can  be  smaller. 

Many  options  exist  for  the  type  and  model  of  the  PEM,  and  a  combination  of  a  rectifier 
and  DC-DC  buck  converter  are  chosen.  This  design  provides  a  two-step  conversion  from  the 
generator  to  the  DC  bus  and  allows  for  good  DC  bus  voltage  regulation  while  reducing 
harmonics  reflected  back  to  the  generator.  A  block  diagram  of  the  system  is  shown  in  Figure  21. 
The  calculations  for  the  rectifier  and  buck  converter  are  performed  for  the  ideal  case  and  ignore 
parasitics  such  as  equivalent  series  resistance,  inductance,  and  capacitance. 
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4.2  Rectification 

The  rectifier  stage  converts  the  AC  voltage  from  the  PM  generator  to  a  DC  voltage. 

There  are  two  main  choices  when  designing  the  rectification  stage:  active  or  passive  rectification 
and  series  or  polyphase  rectification.  Active  rectification  uses  controllable  components,  such  as 
thyristors  or  MOSFETs,  to  actively  rectify  the  AC  voltage.  This  increases  the  cost  and 
complexity  of  the  rectifier  but  provides  voltage  regulation  and  control.  Passive  rectification 
consists  of  diode  bridges  to  convert  the  AC  to  DC  while  relying  on  some  other  circuit  to  perform 
the  voltage  regulation.  A  passive  rectification  scheme  is  used  in  this  paper  because  of  its 
simplicity,  potential  for  lower  input  harmonics,  and  lower  cost. 

Series  rectification  is  when  each  phase  is  rectified  independently  from  the  others  and  the 
output  DC  voltage  is  obtained  by  a  series  connection  of  each  single  phase  rectifier.  Each  phase 
conducts  during  the  full  period  resulting  in  high  power  losses  and  higher  stresses  on  components. 
Also,  the  output  DC  voltage  is  extremely  high  since  each  rectified  phase  voltage  is  added 
together.  Parallel  rectification  is  when  the  lines  with  the  most  positive  and  negative 
instantaneous  voltages  provide  the  forward  bias  to  turn  on  two  diodes.  The  output  voltage 
waveform  corresponds  to  the  instantaneous  difference  between  two  line  voltages.  For  sinusoidal 
voltage  supplies,  the  equation  for  the  output  DC  voltage  is  approximately  Eqn  4-1  [38]. 


where  Vs  =  peak  amplitude  of  phase  voltage 

Eqn  4-1 

The  losses  are  lower  with  parallel  rectification  and  its  output  voltage  provides  a  better 
match  with  the  buck  converter.  Therefore,  a  parallel  rectification  scheme  is  used  in  this  paper. 
For  the  initial  PM  generator  from  section  3.4.2,  the  rectifier  is  shown  in  Figure  22. 
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Figure  22:  Basic  3-phase  Rectifier 

The  DC-side  inductor  and  capacitor  are  necessary  to  obtain  a  good  DC  output  voltage 
with  low  ripple  (large  C)  and  an  acceptable  AC  line  current  waveform  (large  L)  [39].  The  goal  is 
to  produce  maximum  power  factor  with  minimal  line-current  harmonics.  A  power  factor  of  0.96 
with  low  line  harmonics  can  be  achieved  by  selecting  the  inductor  and  capacitor  values  using 
Eqn  4-2  and  Eqn  4-3  [40]. 


where  Lon=0.10 

Vref  =  nominal  l-n  RMS  source  voltage 

■ref  =  P  ref /V ref  (RMS  current  from  V  ref) 

Pref  =  nominal  output  power 
w  =  nominal  source  frequency 


_wj 

''Vcf  *ref  / 


where  Con  =  100.0 

Vref  =  nominal  l-n  RMS  source  voltage 

•ref  =  P  ref ^ref  (RMS  current  from  V  ref) 

Pref  =  nominal  output  power 
fref  =  nominal  source  frequency 


Eqn  4-2 


Eqn  4-3 


63 


For  the  initial  PM  generator  from  section  3.4.2  where  Vref  =  2925  V,  Pref  =  16  MW,  and  fref  =  650 
Hz,  the  values  for  the  DC-side  inductor  and  capacitor  are  Lr  =  82.3  pH  and  Cr  =  287.7  mF. 

These  values  will  change  as  the  system  is  optimized. 


4.3  DC-DC  Conversion 

After  the  rectifier  converts  the  AC  voltage  to  DC  voltage,  the  DC  voltage  must  be 
stepped  down  to  the  bus  voltage  of  700  VDC.  A  high-power  buck  converter  is  ideal  for  this 
application.  The  power  circuit  has  a  basic  topology  to  which  other  circuit  components  are  added 
to  perform  functions  such  as  filtering  and  over-voltage  protection  [41].  The  overall  goal  is  to 
alter  electrical  energy  provided  by  an  input  system  to  that  required  by  an  output  system. 


The  system  consists  of  a  two-stage  input  filter,  a  converter,  and  a  low-pass  output  filter  as 
shown  in  Figure  23.  Duty  ratio  control  using  a  feedback  controller  could  be  used  to  provide  DC 
bus  voltage  regulation  during  normal  operation. 


Figure  23:  Basic  Buck  Converter 

4.3.1  Buck  Converter 

The  buck  converter  is  considered  a  high-frequency  DC-DC  switching  converter.  The 
difference  between  the  switching  frequency  and  the  frequency  of  the  external  waveforms  allows 
the  use  of  low-pass  filters  in  the  input  and  output  to  help  improve  performance  [41].  The  buck 
converter  is  designed  to  take  a  high  DC  voltage  and  reduce  it  to  a  lower  DC  voltage  under 
varying  loads  and  over  a  range  of  input  voltages.  The  level  of  voltage  reduction  is  controlled  by 
the  high  frequency  switching.  The  switch  opens  and  closes  at  a  specific  frequency  with  the  ratio 
of  on-time  to  the  period  defined  as  D,  the  duty  ratio. 


For  the  buck  converter,  the  voltage  conversion  ratio  is  dictated  by  Eqn  4-4. 


Vout 


=  D 


*out  _  1 


Eqn  4-4 

It  is  assumed  that  the  required  output  power  from  the  generator  and  PEM  can  be  as  high  as  the 
rated  16  MW  and  as  low  as  0.1  MW.  Based  on  this,  output  currents  and  resistive  load  values  are 
determined  and  are  listed  in  Table  12. 


Table  12:  Buck  Converter  Load  Values 


PKaielResistance'' . 

0.1  MW 

700  VDC 

142.9  A 

4.90  Q 

16  MW 

700  VDC 

22,857.1  A 

0.031  Q 

The  switching  frequency  for  the  converter  is  selected  to  be  50  kHz  (T  =  20  psec).  This  is 
sufficiently  fast  to  ensure  passive  components  are  not  excessively  large  but  also  not  too  high 
speed  so  that  there  are  not  excessive  switching  losses.  Additionally,  current  switching  ratings  for 
advanced,  high-power  IGBT  devices  are  limited  to  the  hundred  kilohertz  range.  By  selecting  50 
kHz,  this  provides  a  margin  below  this  maximum  switching  frequency  and  ensures  devices  are 
presently  available  or  will  be  in  the  near  future. 


4.3.2  Output  Filter 

A  low-pass  output  filter  is  required  to  help  reduce  the  output  voltage  ripple  to  within  the 
required  specification  listed  in  Table  2.  The  output  voltage  from  the  converter  normally 
fluctuates  between  0  and  Vout,  so  the  filter,  consisting  of  an  inductor  and  capacitor,  is  used  to 
minimize  this  problem.  It  is  important  that  the  comer  frequency  of  the  filter  be  much  lower  that 
the  switching  frequency  (fc  «  fsw)  so  that  the  switching  frequency  ripple  is  eliminated  in  the 
output  voltage  [42].  Normally,  a  factor  of  at  least  lOOx  is  used  [43],  so  to  ensure  an  adequate 
margin  for  parasitics  and  fulfill  the  strict  DC  bus  voltage  ripple  requirement,  a  factor  of  400  is 
selected  so  that  fc  ~  125  Hz.  For  an  LC-filter,  Eqn  4-5  holds. 

tc--T= 

l-TlyfUC 


Eqn  4-5 
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Therefore,  given  the  value  selected  for  fc: 


LC=  1.6210  6 


Eqn  4-6 

With  this  equation  the  output  filter  inductor  and  capacitor  are  sized  using  a  buck 
converter  sizing  program  [44].  It  is  assumed  that  an  inductor  current  ripple  of  0.1%  at  maximum 
power  (16  MW)  is  satisfactory  for  sizing  the  inductor  and  Eqn  4-7  is  used  where  Vp  =  1.9  V,  the 
forward-bias  diode  voltage  [44]. 


1 

u 


-\v 

\  xnjnfci 


W  4  V^\ 
I  in  max  F  I 


/ 


Eqn  4-7 

For  the  initial  generator  of  section  3.4.2,  this  produces  an  inductor  size  of  551  pH.  Using  this 
inductance  along  with  Eqn  4-6,  a  capacitor  value  of  2.94  mF  is  calculated. 

In  order  to  validate  the  filter’s  performance,  a  rough  estimate  for  the  voltage  ripple  is 
calculated.  It  is  assumed  that  the  ripple  component  of  inductor  current  flows  through  the 
capacitor  and  the  average  component  of  the  output  current  flows  through  the  load  resistor  [42]. 
Therefore,  the  peak-to-peak  voltage  ripple  can  be  written  as  shown  in  Eqn  4-8. 

AV  AQ  1  1  aiL  Ts 

°"  c  "c  2  2  2 


Eqn  4-8 

The  output  voltage  ripple  is  estimated  to  be  0.019  V  which  is  within  the  specification  listed  in 
Table  2.  This  overall  process  is  used  later  to  design  the  converter  output  filter  as  the  PM 
generator  system  is  optimized. 


4.3.3  Input  Filter 

The  input  filter  is  needed  to  prevent  the  converter  switching  fluctuations  from  affecting 
the  voltage  and  current  from  the  rectifier.  The  filter  selected  is  a  two-stage  filter  with  an  L-C 
filter  and  an  R-C  damping  leg.  From  the  switching  frequency,  the  angular  frequency  is  co  = 
3.142xl05  rad/sec.  In  order  to  examine  worst  case,  the  16  MW,  22,857  A-output  is  used,  and 
with  the  corresponding  duty  ratio  of  0.1023  for  the  generator  of  section  3.4.2,  an  average  current 
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of  2,338.3  A  is  calculated.  Since  the  fundamental  of  the  ripple  is  attenuated  much  less  than  the 
higher  harmonics,  the  goal  is  to  suppress  the  fundamental  of  the  ripple  along  with  some  margin. 
The  fundamental  of  the  average  current  is  calculated  in  Eqn  4-9. 

4  4 

iconv  1  = <  iconv  !>—  —  *  (2338.3)  =  2911. 2  Amps 

K  K 

Eqn  4-9 

It  is  assumed  that  the  buck  converter  input  current  ripple  must  be  kept  very  small  (<  10  A  at 
maximum  power)  to  ensure  no  detrimental  effects  on  the  rectifier.  Therefore,  at  the  50  kHz 
switching  frequency,  the  maximum  gain  of  filter  is  G  as  calculated  in  Eqn  4-10. 

10.0  Amps 

-  •  ripple  ^ 

G=— = - -  = - 2 - =  0.00168 

iconv 1  2977.2 


Eqn  4-10 

In  order  to  account  for  the  contribution  of  higher  harmonics  and  non-idealities  associated  with 
actual  components,  a  50%  margin  is  used  for  G  resulting  in  G  =  2.52x10 '3. 

The  output  impedance  of  the  input  filter  is  also  a  concern  because  it  must  be  sufficiently 
small  so  that  it  does  not  adversely  affect  the  converter’s  performance.  It  is  assumed  that  the 
output  voltage  ripple  of  the  filter  must  be  less  than  10  V  at  maximum  power.  From  this,  a  value 
of  Zo  at  50  kHz  is  calculated  (Eqn  4-11). 

10.0  V 

Z  =  — - - =  0.00168 

0  2977.2 


Eqn  4-11 

Using  the  same  margin  of  50%  for  the  same  reasons,  the  output  impedance  is  1.12xl0‘3. 

Initially,  the  L-C  low  pass  filter  portion  is  examined  to  calculate  values  for  Lf  and  Cf. 
Then,  a  well-designed  filter  has  Cb  ~  10*Cf,  and  Rf  is  determined  using  reference  [39].  The 
calculations  are  listed  below. 
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Lf 

Iin 

Lf 

Zo  C 

- - o 

7)  Cf  = 

=  & 

> - 

Zo  = 


Zo  = 


ZLZC 
ZL  +  ZC 


s-L 


1  +  s-LC 


G  = 


•in  = 

k n_ 

•out 


ZC  +  ZL 


out 


ZC  +  ZL 


lout 


1.1210  3  = - — 


1  +  s-LC 


2.5210  3  = 


1  +  s-LC 


1.1210  3+  1.1210  3s2LC=s-L 


2.5210  3 +  2.5210  3-s2-L-C=l 


sL=  0.444 
Lf  =  1.415uH 


s  LC  =  395.83 

s-C=  891.51  Cj-  =  2.84  mF 


Eqn  4-12 


Cb  =  n  -Cf  =  10  Cf  =  28.4  mF 


Ro=  I- 

cf 


0.0223 


R 


f  _  (2  +  n)-(4  +  3n) 

Ro  J  2-n2-(4+n) 


Rf  =  0.215  R0  =  4.79  10 
Choose  1 0*Rf  to  be  safe 


R^  =  4.79  mQ 


Eqn  4-13 

The  performance  of  the  input  filter  is  determined  by  examining  the  filter’s  transfer 
function  and  verifying  it  has  sufficient  attenuation  at  the  switching  frequency  and  also  minimal 
peaking.  The  transfer  function  for  the  filter  is  given  in  Eqn  4-14  [45]. 
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( 

'Rfcbs) 

1+1 

(  i\ 

(LfCfCbRfs J 

|  + 

Lf(cf  +  Cc)-s2 

+  (r  fCb-s)  +  1 

Eqn  4-14 

Figure  24  shows  the  Bode  plot  for  the  initial  input  filter  and  it  operates  satisfactorily  (Appendix 
F  contains  the  detailed  MATLAB  code).  This  overall  procedure  for  designing  the  converter 
input  filter  is  utilized  later  to  as  the  PM  generator  system  is  optimized. 


Figure  24:  Bode  Plot  for  Converter  Input  Filter 

4.3.4  Converter  Control 

One  control  technique  that  can  be  used  is  duty  ratio  control  utilizing  a  feedback 
controller.  This  arrangement  is  presented  as  an  illustration  of  one  control  methodology  for  the 
power  electronics  module  but  implementation  is  beyond  the  scope  of  this  paper. 

The  control  scheme  involves  negative  feedback  control  using  linearized,  averaged  state- 
space  models.  Figure  25  shows  a  block  diagram  of  this  where  K  represents  the  controller  and  H 
represents  the  buck  converter.  It  is  designed  so  that  the  closed-loop  system  is  stable  and  well 
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damped  for  the  input  voltages  and  output  loads.  Voltage-mode  control  (duty-ratio  control)  is 
implemented  using  a  Proportional  Integral  Derivative  (PID)  controller. 


Figure  25:  Block  Diagram  of  Control  Scheme 


The  majority  of  control  systems  used  in  industry  applications  are  based  on  this  type  of 
controller.  In  general,  proportional  and  integral  control  can  be  applied  independently  or  together 
(P,  I,  or  PI),  while  derivative  control  can  only  be  applied  in  combination  with  one  or  both  of  the 
other  two  (PD,  PID).  Testing  of  a  PID  controller  normally  involves  plotting  the  system  response 
to  a  simple  input,  such  as  a  step  input,  and  adjusting  the  gains  based  on  the  output. 

The  first  step  in  developing  the  control  scheme  is  to  obtain  a  linearized,  averaged  model 
for  the  buck  converter  and  derive  the  transfer  function  from  perturbation  in  duty  ratio  to 
perturbation  in  output  voltage.  The  transfer  function  represents  H  in  Figure  25.  Next  the  PID 
controlled  is  developed,  consisting  of  three  gains,  Kp,  Ki,  and  Kd,  and  this  represents  K  in  Figure 
25.  The  H  and  K  are  then  combined  to  form  the  complete  transfer  function  for  the  system. 

Initially,  the  state  variables  are  assigned  (iL  and  vc)  and  the  state-space  equations  are 
determined  (neglecting  the  input  filter). 

L'jiL  =  (V1  “  -  VC'*1  ~  °7VC  =  I  ‘L"  “7  ‘q(t)  “  f‘L"  vl(1  "  q(t)) 

at  at  \  K  j  \  K  j 


(vrq(t)  -  vc) 


Eqn  4-15 

Then,  the  circuit  is  averaged  and  it  is  assumed  that  the  state  variables  have  slow  variation  and 
small  ripple,  and  d  =  avg(q(t)).  All  of  the  variables  are  now  averaged  quantities  and  Eqn  4-15 
becomes  Eqn  4-16. 
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Eqn  4-16 


dt^L  =  l  (V1  d  -  VC> 


d  1 
— vr  =  — 

dt  C  C 


(  VC^ 

vlL  R  J 


The  variables  are  linearized  resulting  in  the  equations  shown  in  Eqn  4-17. 


vC=  Vc+  v'c 


‘L  =  lL  +  *L 
Linearized  Equations 


d  =  D  +  d’ 


V1  =  Vi  +  v’j 


dt  L  L  1  L  L  '■ 


d  ,  1  1  , 

— v  r  =  — It - v  r 

dt  C  C  L  R  C  c 


Eqn  4-17 

Since  vc  is  the  output  voltage,  utilizing  the  Laplace  operator  and  solving  for  the  transfer  function 
yields  Eqn  4-18. 


2  1  1  > 

s  + - s  + - 

RC  LC  J 


out 


H  = 


d  ,  vi 

- v  i  + - d 

LC  1  LC 


V, 


out  LC 


d’  2  1  1 

s  + - s  + - 

RC  LC 


For  the  controller  the  transfer  function  is  Eqn  4-19. 

K I  where  e'  is  Vref  -  Vout 


£ 

e' 


K  =  —  =  Kp  +  Kds  +  — 


Eqn  4-18 


K  = 


KjyS  +  Kp-  s  +  Kj 


Eqn  4-19 


Overall,  the  total  transfer  function  is  then  Eqn  4-20. 
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=  H-K 


v  out 
e' 

v'out  I^(KDs2+KPs  +  Ki) 


e'  3  1  2  1 

s  + - s  + - s 

RC  L-C 


Eqn  4-20 

Using  the  total  system  transfer  function,  the  gains  on  the  controller  are  adjusted  to  obtain 
the  desired  steady-state  and  transient  operation.  This  type  of  system  can  readily  be  designed  and 
optimized  to  provide  the  required  overall  performance,  but  this  is  beyond  the  scope  of  this  paper. 


4.4  Conversion  Losses 

An  integral  part  of  a  DC-DC  converter’s  performance  is  how  efficiently  it  operates 
during  the  energy  conversion  process.  Losses  that  occur  include  winding  losses  due  to 
conduction,  switching  losses,  and  inductor  core  losses.  A  detailed  analysis  of  all  of  these  losses 
is  beyond  the  scope  of  this  paper  and  they  are  only  described  qualitatively  and  estimated  using 
analytical  equations. 

Practical  magnetic  components,  such  as  inductors,  exhibit  losses  generated  by  winding 
conduction,  hysteresis,  and  eddy  currents.  Resistive  components  in  the  converter  as  well  as 
parasitics  (equivalent  series  resistance)  produce  I2R  losses.  Ideal  inductors  and  capacitors  are 
employed  in  this  paper  and  parasitic  losses  are  not  calculated. 

The  transistors  and  diodes  have  conduction  losses  based  on  their  on-state  resistances. 
The  equations  to  estimate  these  losses  are  shown  in  Eqn  4-2 land  Eqn  4-22.  Lastly,  the 
transistors  have  switching  losses  associated  with  the  transistors  changing  states.  If  the  voltage 
and  current  waveforms  are  assumed  to  rise  and  fall  approximately  linearly,  the  switching  losses 
are  determined  as  shown  in  Figure  26. 


where  lL  =  average  current 
D  =  duty  cycle 

Rds  =  IGBT  on-state  resistance 

Eqn  4-21 
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pdiode  =  0  -  D>  >L  Rd 

where  lL  =  average  current 
D  =  duty  cycle 

Rd  =  diode  on-state  resistance 

Eqn  4-22 


4.5  Component  Sizes  and  Weights 

The  PEM  detailed  in  the  previous  sections  consists  of  IGBTs,  diodes,  capacitors, 
inductors,  and  resistors.  The  sizes  and  weights  of  the  resistors  are  neglected  because  they  are 
significantly  smaller  than  the  other  components.  To  ascertain  reasonable  estimates  for  the  other 
devices,  numerous  references  are  utilized  [15],  [46],  [47],  [48],  [49],  [50],  [51],  [52],  [53],  [54]. 
In  addition,  a  30%  service  fraction  is  added  to  the  total  component  mass  and  volume  calculations 
to  account  for  heat  sinks,  cooling,  etc. 

The  power  transistor  is  the  main  switching  device  in  the  converter  circuit.  In  many 
applications,  IGBTs  are  used  because  they  have  the  current  density  and  low  loss  of  bipolar 
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transistors  with  the  high-speed  and  high  input  impedance  of  MOSFETs  [54].  Power  diodes  are 
needed  to  complete  the  buck  converter  circuit  and  provide  proper  operation.  Table  13  lists  the 
assumed  device  parameters  for  the  IGBTs  and  diodes  used.  The  IGBT  and  diode  values  are  not 
from  exact  devices  from  particular  manufacturers  but  rather  represent  data  assembled  from 
numerous  sources.  In  order  to  achieve  the  high  power  levels,  the  semiconductor  devices  must  be 
connected  in  parallel.  Specifically,  since  the  IGBTs  and  diodes  are  rated  for  1,200A  each,  19 
modules  (each  rated  for  845  kW)  are  connected  in  parallel  to  achieve  the  required  22.5  kA.  This 
can  be  readily  accomplished  using  IGBTs  provided  the  switching  is  properly  controlled  [55]. 

The  inductors  and  capacitors  are  needed  to  perform  energy  storage  and  filtering 
functions.  In  high-power  applications  they  are  usually  the  dominant  components  in  terms  of 
mass  and  volume.  For  estimating  their  contribution  to  the  converter,  component  weight  and 
volume  energy  densities  are  utilized.  Table  13  lists  the  assumed  characteristics  for  the  inductor 
and  capacitor  sizing  estimations. 


Table  13:  Power  Electronics  Module  Component  Characteristics 


Maximum  Collector  Emitter  Voltage 

6500  V 

Maximum  Collector  Current 

1200  A 

Mass  (including  heat  sink) 

1200  g 

Volume  (including  heat  sink) 

Turn  On  /  Rise  Time 

450  ns 

Turn  Off /Fall  Time 

400  ns 

On-State  Resistance 

0.06  Q 

msm mmsssmmm^am 

1.9  V 

Maximum  Switching  Frequency 

-100  kHz 

Peak  Blocking  Voltage 

6500  V 

Average  Forward  Current 

1200  A 

Mass  (including  heat  sink) 

450  g 

Volume  (including  heat  sink) 

IBBSIBflH 

On-State  Resistance  (estimated) 

0.0073  Q 

be mmssmmmmm 

1.0  V 

.  ! '  _  . .. _  ..  . 1 

Mass  Density 

■IIH1I39HHU 

8  J/cc 

Mass  Density 

3  kg/J 

Volume  Density 

0.08  J/cc 
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In  addition  to  using  the  component  values  to  size  the  power  conversion  module,  a  weight  factor 
(l.Ox)  and  a  volume  factor  (0.20x)  are  employed  to  account  for  portions  such  as  structure, 
frames,  mounts,  and  maintenance  access. 
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Chapter  5  Waveforms,  Models,  and  Machine/Module  Optimization 


5.1  Initial  Generator  Waveforms 

The  PM  generator  produces  back  EMF  waveforms  that  are  dependent  on  a  number  of 
factors  as  discussed  in  Chapter  3.  The  goal  is  to  produce  a  voltage  waveform  that  closely 
resembles  a  sinusoidal  waveform  with  a  low  total  harmonic  distortion  (THD)  because  this  results 
in  minimal  harmonic  content  which  reduces  losses  in  the  machine.  THD  is  a  measure  of  the 
distortion  in  a  waveform  caused  by  undesirable  frequency  components.  It  is  calculated  as  shown 
in  Eqn  5-1. 

THD  = 

Eqn  5-1 

The  back  EMF  waveforms  are  generated  using  Eqn  3-18,  Eqn  3-19,  and  Eqn  3-22  which 
are  shown  again  below. 


oo 

B(e)=  Yj  Bn  •sin  (npe) 


n  =  1 
n  odd 


where 


4  ,  .  nP0m  . 

B  = - B  k  sin  — —  sin 

n  mi  g  ^  2  J 

6m  =  magnet  physical  angle 

n  =  harmonic  number 


Eqn  5-2 


oo 

x{q)=  y  xn  •sin(np0) 


n  =  1 
n  odd 


where 


^•^t-Na-VW^n 

P 


Eqn  5-3 
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n  =  1 
n  odd 


where  y  =  !Lx  =  an.x 
Yn  ,/n  0  n 

at 


Eqn  5-4 

MATLAB  code  is  developed  to  create  the  generator  EMF  waveforms  as  well  as 
determine  the  harmonic  content  and  THD  of  the  generator  output  waveform.  This  code  is 
contained  in  0.  For  the  initial  PM  generator  designed  in  section  3.4,  the  air  gap  flux  density, 
back  EMF  waveform,  and  harmonic  content  (line-neutral)  are  shown  in  Figure  27,  Figure  28,  and 
Figure  29.  As  a  frame  of  reference,  MIL-STD-1399  requires  Type  1 60  Flz  AC  systems  to  have  a 
total  THD  of  less  than  5  %  line-line  [56].  This  limit  does  not  directly  apply  to  the  PM  generator 
AC  bus  since  it  is  being  rectified  to  DC  before  supplying  loads.  However,  it  is  used  as  a 
guideline  because  it  provides  a  reasonable  limit  to  ensure  proper  rectification  and  performance 


Figure  29:  Initial  Generator  Harmonic  Content 


5.2  Rotational  Stress  and  Retaining  Sleeve 

Since  the  PM  generator  is  spinning  at  high  speed,  the  rotor  and  permanent  magnets  are 
subjected  to  extremely  high  centrifugal  forces.  These  forces  can  cause  significant  amounts  of 
damage  if  the  magnets  and  rotor  are  not  properly  restrained.  The  rotational  components  can  be 
strengthened  by  enclosing  them  in  a  retaining  sleeve/can  which  also  increases  the  air  gap  length, 
The  centrifugal  force  on  the  magnets  due  to  the  rotor  spinning  is  calculated  in  Eqn  5-5. 

M  -v  2 
m  vmag 

C6n=  R  +  hm 

where  Mm  =  mass  of  magnets 

vmag  =  velocity  of  magnets  (m/s) 


Eqn  5-5 

Using  the  inner  surface  area  of  the  retaining  sleeve,  this  force  is  converted  to  an  outward 
pressure.  Treating  the  retaining  sleeve  as  a  thin-walled  vessel,  the  hoop  stress  felt  by  the  sleeve 
is  determined  as  shown  in  Figure  30  and  Eqn  5-6. 


Figure  30:  Retaining  Sleeve  Hoop  Stress 

rn 

V  Fv  =  -2  cc  h  L+  P  r  L  sin(e)  Ld0  =  0 

'  J  r\ 


vert 


Pr 


°c  = 


80 


Eqn  5-6 


To  provide  a  margin  for  mechanical  tolerances  and  imperfections,  a  safety  factor  (SF)  of  1.2  is 
applied  to  the  hoop  stress  calculated  in  Eqn  5-6  to  get  a  final  SF  stress  for  the  retaining  sleeve. 

The  retaining  sleeve  can  be  made  from  many  different  types  of  materials  including  metal 
alloys  and  composites.  A  disadvantage  of  a  metallic  sleeve  is  eddy  currents  are  induced  in  the 
sleeve  by  variations  in  the  flux  density  caused  by  the  stator  slots.  The  advantages  of  a  metallic 
sleeve  are  that  it  shields  the  magnets  from  most  of  the  flux  density  variations  and  it  has  a  high 
thermal  conductivity  for  heat  removal. 

A  composite,  such  as  carbon-fiber,  provides  reasonable  strength  while  having  lower 
losses  since  the  lower  conductivity  reduces  the  eddy  currents.  However,  the  composite  sleeve 
has  low  thermal  conductivity  and  does  not  shield  the  magnets  from  the  flux  variations.  This 
results  in  increased  losses  in  the  magnets  themselves.  In  general,  when  the  sleeve  conductivity  is 
low,  the  rotor  losses  are  in  the  permanent  magnets,  while  when  the  sleeve  conductivity  is  high, 
the  PM  losses  decrease  and  the  sleeve  losses  increase.  The  materials  considered  for  use  as  in  the 
retaining  sleeve  are  listed  in  Table  14  [10],  [57],  [58],  [59],  [60],  [61],  [62],  [63],  [64].  The 
numbers  for  the  carbon  fiber  composite  are  in  the  middle  of  a  range  of  typical  values  (0.5  -  18.0 
pQ-m)  since  wide  variations  exist  depending  on  how  the  composite  is  made. 


Table  14:  Retaining  Sleeve  Materials 


ij  |  HU  Material';- 

Stainless  Steel 

90 

0.72 

Aluminum  Alloy 

75 

0.05 

Titanium  Alloy 

110 

0.78 

Inconel 

132 

0.98 

Carbon  Fiber 

100 

9.25 

0  contains  detailed  MATLAB  code  to  perform  the  retaining  sleeve  stress  calculations. 


The  initial  PM  generator  from  section  3.4.2  is  analyzed  with  the  results  below. 


Retaining  Sleeve  Stress: 

Stress  Limits: 

Stainless  Steel  =  90.0  ksi 

Titanium  Alloy  =  1 10.0  ksi 

Inconel  =  132.0  ksi 

Actual  Sleeve  Stress: 

Sleeve  Thickness  Actual  Stress 
0.50  mm  1136.9  ksi 

1.00  mm  568.4  ksi 

1.50  mm  379.0  ksi 


Aluminum  Alloy  =  75.0  ksi 

Carbon  Fiber  =  100.0  ksi 


SF  Stress 
1364.2  ksi 
682.1  ksi 
454.7  ksi 


81 


2.00  mm 

284.2  ksi 

341.1  ksi 

2.50  mm 

227.4  ksi 

272.8  ksi 

3.00  mm 

189.5  ksi 

227.4  ksi 

3.50  mm 

162.4  ksi 

194.9  ksi 

4.00  mm 

142.1  ksi 

170.5  ksi 

4.50  mm 

126.3  ksi 

151.6  ksi 

5.00  mm 

113.7  ksi 

136.4  ksi 

5.50  mm 

103.4  ksi 

124.0  ksi 

6.00  mm 

94.7  ksi 

113.7  ksi 

6.50  mm 

87.5  ksi 

104.9  ksi 

7.00  mm 

81.2  ksi 

97.4  ksi 

7.50  mm 

75.8  ksi 

90.9  ksi 

8.00  mm 

71.1  ksi 

85.3  ksi 

It  is  evident  that  the  initial  4  mm  gap  length  is  not  sufficient  once  the  retaining  sleeve  is 
considered  (since  the  air  gap  dimension  includes  the  retaining  sleeve  thickness).  The  machine 
must  therefore  be  revised  to  allow  for  the  retaining  can.  In  this  modification  process,  the 
aluminum  alloy  is  removed  as  a  sleeve  material  possibility  since  it  has  lower  strength  and  much 
higher  conductivity  than  the  other  materials  which  together  would  produce  a  larger  air  gap  and 
higher  eddy  current  losses. 

The  PM  generator  will  be  redesigned  taking  the  retaining  sleeve  and  hoop  stress  limits 
into  consideration.  The  updated  PM  generator  will  be  a  bigger  machine  with  a  much  larger  air 
gap,  greater  magnet  height,  lower  Bg,  lower  voltage,  and  higher  current  density. 

5.3  Rotor  Losses 

5.3. 1  Model  for  Time  Harmonics  &  Winding  Space  Harmonics 

The  permanent  magnets  used  in  high-speed  generators  are  electrically  conductive  and 
therefore  support  eddy  currents.  The  retaining  sleeves  are  sometimes  made  from  electrically 
conductive  material  that  also  can  carry  eddy  currents.  These  eddy  currents  are  primarily  caused 
by  fluctuations  in  the  magnetic  flux  density  produced  by  time  and  space  harmonics  of  the 
winding  currents.  The  currents  produce  losses  which  can  potentially  cause  excessive  heating  or 
demagnetization  of  the  permanent  magnets.  An  analytical  model  is  developed  using  the  winding 
and  current  harmonics  and  the  surface  impedance  to  estimate  the  rotor  losses. 
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Figure  31  shows  the  PM  generator  geometry  “flattened  out”  into  rectilinear  coordinates. 
This  is  an  accurate  representation  provided  the  dimensions  are  on  a  radial  scale  that  is  much 
smaller  than  the  radius  of  the  machine  so  that  curvature  is  not  important  [65]. 


Magnets 

Rotor  Hub 


Figure  31:  General  Magnet  Loss  Model 


The  direction  of  rotation  is  in  the  positive  x-direction,  the  radial  direction  is  y,  and  the  armature 
current  flows  in  the  axial  dimension,  z.  In  addition  to  the  arrangement  of  Figure  31,  the 
geometry  shown  in  Figure  32  is  also  utilized. 


Figure  32:  Layer  of  Material 


The  following  assumptions  are  made  in  developing  the  rotor  analytical  loss  model: 

•  Layers  of  material  extend  to  ±  infinity  in  the  ±  x  direction 

•  Layers  effectively  extend  to  negative  infinity  in  the  negative  y  direction 
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•  Motion/rotation  is  in  the  +  x  direction 

•  The  physical  constants  of  the  layers  are  homogeneous,  isentropic,  and  linear 

•  The  ferromagnetic  material  does  not  saturate 

•  The  machine  is  long  axially  so  magnetic  variations  in  the  z  direction  are  ignored  (H  and 
B  only  vary  in  x,  y  directions) 

•  All  currents  flow  in  the  z  direction 

•  The  rotor  and  stator  are  constructed  of  laminated  steel  so  their  conductivity  in  the  z 
direction  is  negligible 

•  The  time  and  space  variations  are  approximately  sinusoidal 

•  Flux  density  at  y  =  infinity  is  zero 

•  A  traveling  flux  wave  harmonic  can  be  represented  by  an  equivalent  traveling  current 
sheet  (Kz)  on  the  surface  of  the  stator 

•  The  normal  component  of  the  flux  density  is  continuous  at  all  interfaces 

•  The  tangential  component  of  the  flux  density  is  continuous  at  all  interfaces  except  at  the 
stator/air  gap  where  it  is  increased  by  the  current  sheet  density 

•  The  magnetic  flux  density  crossing  the  air  gap  and  magnets  is  perpendicular 

•  The  effect  of  magnet  eddy  currents  on  the  magnetic  flux  density  is  negligible  -  this  is 
accurate  below  10  kHz  [66] 

•  The  magnet  flux  density  is  constant  over  the  magnet  breadth 

One  of  Maxwell’s  Law  states  that  there  is  no  magnetic  charge  and  therefore  flux  lines 
must  close  on  themselves,  represented  by  Eqn  5-7. 

-> 

VB  =  0 

Eqn  5-7 

Knowing  this,  the  flux  density  B  can  be  represented  as  the  divergence  of  the  magnetic  scalar 
potential  A  (Eqn  5-8). 

->  -> 

B  =  V  x  A 

Eqn  5-8 

Substituting  in  poH  for  B,  taking  the  gradient  of  Eqn  5-8,  and  using  Ampere’s  Law  for  current 
density  yields  Eqn  5-9. 
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p0-H  =  V  x  A 


|1qV  xH=VxVxA 


2 

|^0-V  x  H  =  V\V  A/  -  V  A 

»  2~> 

(Jq-J  =  -V  A 

■*  2~* 

Hq-g-E=  -V  A 

Eqn  5-9 

Faraday’s  Law  is  then  employed  to  determine  a  relationship  between  E  and  A  as  shown  below. 

j  ■) 

V  x  E  =  — B  (Faraday's  Law) 

dt 

*)  j  “ ^ 

V  x  E  =  -V  x  —  A 

dt 

■)  A 

E  =  —  A 


Combining  Eqn  5-9  and  Eqn  5-10  gives  Eqn  5-11. 


Eqn  5-10 


2  “7  j  "7 

V  -A  =  p  0  G  —  A 
u  dt 

Eqn  5-11 

Since  H  and  B  vary  only  in  the  x  and  y  directions,  A  and  E  are  in  the  z  direction  and  A  can  be 
represented  by  Eqn  5-12. 


Az=Re[(Az(y)).eJ<“‘-kx)] 

where  k  =  2  %/X 


Eqn  5-12 


Performing  the  mathematical  operations  and  solving,  Eqn  5-12  becomes  Eqn  5-13. 
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2  d 

V  ’Ay  —  (JLn*0> - A  7 

z  u  dt  z 


^  1  _J 

Since  b  =  V  x  A  and  A  is  in  z-dir,  then  by  =  —  A7  and  bv  =  — 


x  A  z 
dy 


dx 


v2  - 


( d2  a2  d2  A 

+  —i 


dx  dy  dz  ) 


.2  d* 

=  -k  +  — 

dy 


-k2-Az(y)  +  -^Az(y)  =  j(OH0oAz(y) 
dy 


"~2 Az(y)  -  (k2  +  j-co-|x0-o J-Az(y)  =  0 
dy 


Let  y  =  _^k~  +  j-co-^Q-a 
Solution  of  the  form: 

jy  ,  a 


Az(y)  =  Ap-e  +An-e 
{( 


A,  =  Re||A„.e™  +  Ane-W}=i<“,-b<^ 


Solving  for  the  magnetic  flux  densities  and  magnetic  fields  produces  Eqn  5-15. 

yy^  ej-(cot-kx) 


Bx  =  “Az  =  Re[(Y  Ap  eYy  “Y  An  e 


By  =  ~ Az  =  Rel 


dx 


jy 


Hx=Re| 


y-A 


p.e^y  _ 


j'kApe  +  j'kAne  J'e 

rAn  -  A 


-  yy  i  j-(cot-kx) 


yy 


L  v 


Po  Po 


j-(cot-kx) 


J 


Hy  =  Re| 


jk  AP  yy  jk  An  - 

- L.gi-'  + - e 


yy 


L  v 


PO  PO 


i(cot-kx) 


Eqn  5-13 


Eqn  5-14 
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Let 


Y'Ap  .  "YAn 

Hp  =  —  and  H„ - 


R0 


t>0 


Hx=R{(Hp' 

Hy  *  R{(jfHP 


eW  +  Hne 


■  kx) 


-  H  |  -  Hr  e 

Y 


-  yy  I  i(cot-kx) 


Hv  =  Re|  Hx  e 


j(cot-kx) 


Hy  =  Re| 


Hye 


j(cot-kx) 


Eqn  5-15 


The  layer  from  Figure  32  can  now  be  examined  utilizing  a  “surface  coefficient.”  A  surface 
coefficient  is  defined  as  the  ratio  of  the  y-directed  to  x-directed  magnetic  field  amplitude  (a  = 
Hy/Hx).  At  the  bottom  of  the  layer  where  y  =  0,  the  surface  coefficient  is  given  by  Eqn  5-16. 

,  Hy  k  Hp  -  Hjj  (“p 


H 


X  Y  Hp  +  Hn  ocb  =  j — 

Y 


«P 

Hn 


+  1 


Eqn  5-16 


At  the  top  of  the  layer  (y  =  h),  the  surface  coefficient  is  shown  in  Eqn  5-17. 


f 


at(y  =  h)  =  j— 
Y 


V1’’  -  Hn' 


-  yh 


,Yh 


Hpe  +V 


-  yh 


.  k 

«t  =  J— 
Y 


f  i_r  \ 

_p  eYh_e-Yh 
H„ 


”p  yh  -  yh 
— -e'  +  e  ' 


VHn 


J 


Eqn  5-17 

Using  Eqn  5-16  to  solve  for  the  ratio  Hp/Hn,  a  final  expression  for  the  top  surface  coefficient  is 
determined  (Eqn  5-18),  and  it  is  applicable  to  any  uniform  region. 
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«  i+I  =j± 

b  H„  T  «„ 


Hpf.k  A  .k 

—  J--«b  =J-  +  «b 
Hn  V  Y  J  Y 


(  jk 
—  +  ab 

_1 _ eYh  _  e-  Yh 

jk 

. - ab 

.  k  y 

«t  =  J-  — - - 

Y  i±  +  „ 

+  ab 

A - eYh  +  e"Yh 


H  J”  +  ab 
p  _  y 

Hn  ”  .  k 

n  J—  "«b 
Y 


[  j  k  )  yh  f  jk  )  -yh 

a  =  ;.l  _W _ )_ _ L 1 _ l _ 

1  y  f  j  k  ^  yh  f  j  k  ^  -  Xh 

—  +  ab  e'  +  - - ab  e 

_W  J  Y  J 


j — sinh  (yli)  +  abcosh  (yh) 

.  k  y 

«t  =  J—  — - 

Y  j — cosh  (yli)  +  ab-sinh  (yli) 

Y  J 

Eqn  5-18 

If  the  region  being  examined  is  positioned  on  top  of  a  ferromagnetic  surface,  such  as  the  magnets 
on  the  steel  rotor  shaft,  the  boundary  condition  at  the  bottom  of  the  layer  (ab  —*■  infinity  as  Hx  — ► 
0)  produces  Eqn  5-19. 


a  t  =  j— coth(yh) 

Y 


Eqn  5-19 


In  the  case  of  the  air  gap  where  the  conductivity  is  zero,  Eqn  5-18  reduces  to  Eqn  5-20. 
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«t  =  > 


f  j-sinh(Th)  +  a^-cosh^)^ 
jcosh(')ti)  +  a  ij-sinh(')h) 


Eqn  5-20 

The  electric  field  (E)  is  z-directed  and  the  magnetic  flux  density  and  field  (B  and  H)  are 
y-directed  so  a  relationship  is  determined  between  them  using  Faraday’s  Law.  The  surface 
impedance  is  then  determined.  Surface  impedance  (Zs)  is  the  ratio  of  the  z-directed  electric  field 
to  the  z-directed  current.  Eqn  5-21  and  Eqn  5-22  show  the  detailed  calculations. 


V  x  E  =  — B 
dt 


(Faraday's  Law) 


V  x  E  =  -Un  — H 
u  dt 

7ex  -  xEz  ■  -t*0~:Hy 

dz  dx  dt  J 


jkEz=  -jcop0-Hy 
co„ 


k. 


-HO'Hx 


Eqn  5-21 


Z.  = 


Esurf  Ez  kjj 


-P()Hy 


^surf  Hx 
C0n 

Zs  =  —^0a 

% 


-H„ 


Eqn  5-22 

Since  electromagnetic  power  flow  into  the  rotor  is  the  desired  quantity,  Poynting’s  vector 
is  employed  (Eqn  5-23). 

*  ■>  -> 

S  =  Ex  H 


Eqn  5-23 


Dissipation  in  the  rotor  is  in  the  negative  y  direction  producing  Eqn  5-24. 
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-syHEz'Hx-E*Hz) 

There  is  no  z-directed  magnetic  field  (Hz  =  0) 
-Sy  =  -%HX 

In  time  average  form,  Eqn  5-24  becomes  Eqn  5-25. 

-<  Sy>=  -"Re(Ez  •  Hx) 

-<  Sy>=  ~—\Hx\2  ■  Re(— ) 

v  r  t  J 


Eqn  5-24 


-<  Sy>=  •  Re(Zs) 

Eqn  5-25 

Eqn  5-25  yields  the  power  dissipated  at  the  stator  surface.  This  is  the  correct  result  for 
the  rotor  because  there  is  no  mechanism  for  dissipating  power  between  the  stator  and  rotor.  The 
power  estimated  by  Poyn ting’s  theorem  flows  directly  from  the  stator  to  the  rotor  [65]. 

In  order  to  calculate  the  rotor  losses,  the  above  analytical  model  is  applied  to  the 
geometry  in  Figure  31.  For  this  model,  the  stator  is  assumed  to  be  a  smooth  surface  without  slots 
because  the  slot  effects  are  considered  later  in  section  5.3.2.  The  first  step  is  calculating  the 
surface  coefficient  at  the  bottom  of  the  magnet  layer.  It  is  assumed  that  this  is  formed  by  the 
highly  permeable  rotor  shaft  below  the  magnets.  This  assumption  allows  the  surface  coefficient 
at  the  top  of  the  magnet  layer  to  be  calculated  using  Eqn  5-19. 

Next,  the  surface  coefficient  at  the  top  of  the  retaining  sleeve  is  determined  using  the 
sleeve  material’s  conductivity  and  Eqn  5-18.  Traveling  across  the  air  gap  to  the  surface  of  the 
stator,  the  surface  coefficient  is  estimated  using  Eqn  5-20.  The  surface  impedance  is  then 
computed  utilizing  Eqn  5-22. 

Once  the  surface  impedance  is  known,  the  rotor  losses  caused  by  the  stator  winding  time 
and  space  harmonics  are  calculated  using  Eqn  5-25.  Only  the  more  significant  harmonics  are 
considered  as  shown  in  Figure  33.  The  higher  order  harmonic  effects  are  ignored  because  the 
products  of  higher  order  time  harmonics  (small  number)  and  higher  order  space  harmonics 
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(small  number)  produce  negligible  contributions  (small  number  *  small  number  =  very  small 
number).  Appendix  I  contains  the  detailed  MATLAB  code  that  performs  all  of  the  above 
calculations. 


Figure  33:  Relevant  Harmonics 


5.3.2  Model  for  Stator  Slot  Effects 

The  stator  slots  cause  variations  in  the  magnetic  field  which  produce  losses  in  the 
retaining  sleeve  and  magnets  of  the  rotor.  Accurate  calculation  of  the  losses  in  the  retaining 
sleeve  is  extremely  difficult.  Several  different  methods  have  been  developed  and  in  this  paper, 
the  technique  from  reference  [18]  is  employed.  Figure  34  shows  the  geometry  of  the  retaining 
can  and  the  currents  that  are  induced. 


Figure  34:  Retaining  Sleeve  Induced  Currents 
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As  the  rotor  spins  past  the  slot  openings  of  the  stator,  the  air  gap  flux  density  undergoes 
modulation  due  to  the  change  in  reluctance.  The  dip  in  Bg  (shown  in  Figure  35)  travels  along  the 
B-waveform  which  is  otherwise  moving  synchronously  with  the  rotor. 


Figure  35:  Flux  Density  Variation 

The  rotation  of  the  rotor  generates  an  E-field  in  the  sleeve  and  a  subsequent  axial  current 
density  J  =  E/p.  When  this  current  density  is  integrated  over  the  volume  of  the  can,  the  average 
loss  per  unit  area  is  determined  (Eqn  5-26)  [18]. 

w  = 

Eqn  5-26 

It  is  evident  from  Eqn  5-26  that  as  the  slot  width  increases,  the  width  of  the  flux  density  dip  (P) 
gets  larger  causing  the  sleeve  losses  to  increase.  The  above  equation  only  considers  eddy 
currents  flowing  in  the  axial  direction  but  there  are  also  circumferential  components.  These 
portions  are  accounted  for  using  a  factor  Ks  as  shown  in  Eqn  5-27  where  the  total  can  losses  are 
determined  [18]. 


Z.[B  N  (R  +  M  »  where 

3600  p  ^2  > 
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Pcan  =  Kg-  w-A 


where  a  =  *2.(r  +  h  J-Lst 


tanh 

PLst 

P'Lst 

2(R  +  hm)_ 

Eqn  5-27 


One  way  to  reduce  the  retaining  can  losses  is  to  split  the  sleeve  cylinder  into  separate 
rings.  This  reduces  the  effective  flow  path  length  for  the  eddy  currents  thus  lowering  the  losses. 
There  is  a  limit  to  how  much  this  can  be  done  since  as  the  number  of  rings  increases,  the 
construction  becomes  more  difficult.  In  this  paper  it  is  assumed  that  the  maximum  feasible 
number  of  sleeve  sections  is  ten. 

When  metal  alloys  are  used  as  the  retaining  sleeve  material,  they  effectively  shield  the 
magnets  from  the  flux  density  variations,  whereas  there  is  little  or  no  shielding  with  the  carbon 
fiber  composite  because  of  its  low  conductivity.  Therefore,  with  a  metal  alloy  can,  the  rotor 
losses  occur  largely  in  the  can  compared  to  the  carbon  fiber  sleeve,  where  the  losses  mostly 
occur  in  the  magnets.  These  magnet  losses  are  calculated  using  methods  similar  to  Eqn  5-26  and 
Eqn  5-27  assuming  that  the  eddy  currents  flow  in  the  top  10%  of  the  magnet  volume.  Appendix 
J  contains  the  MATLAB  code  that  performs  the  detailed  rotor  loss  calculations  caused  by  the 
stator  slot  effects. 


5.4  Complete  System  Model  &  Design  Procedure 

The  complete  model  for  the  power  generation  module  consists  of  the  PM  generator,  the 
rectifier,  the  filtering/energy  storage  components,  the  buck  converter,  and  the  load.  The  program 
PS1M  (Version  6.0,  Powersim,  Inc.)  is  utilized  to  construct  the  circuit  model,  perform  detailed 
time  series  and  waveform  analysis,  and  determine  voltage  and  current  values.  An  example  of  a 
7-phase  version  of  this  model  system  is  shown  in  Figure  36.  In  each  of  the  models,  the  converter 
and  rectifier  are  simulated  as  one  IGBT/diode  module  for  simplicity  but  for  actual 
implementation/construction,  19  modules  are  required  as  described  in  section  4.5. 
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GENERATOR  RECTiriER  INPUT  FILTER  CONVERTER  OUTPUT  FILTER/ OUTPUT 

Figure  36:  System  Model 

The  procedure  for  designing  and  analyzing  the  power  generation  module  consists  of  the 
following: 

•  Design  a  PM  generator  using  the  methodology  from  sections  3.4.2  and  5.2,  iterating  to 
ensure  the  sizes,  current  density,  magnetic  saturation,  and  sleeve  thickness/stresses  are  all 
satisfactory 

•  Verify  the  PM  generator  has  proper  output  waveforms  and  THD  per  section  5. 1 

•  Calculate  the  rotor  losses  caused  by  stator  slot  harmonics  as  outlined  in  section  5.3.2 

•  Determine  the  rectifier  circuit  parameters  utilizing  the  procedures  from  section  4.2 

•  Design  the  buck  converter  and  associated  input  and  output  filters  employing  the  methods 
of  section  4.3 

•  Construct  the  PSIM  circuit  using  the  correct  generator  parameters  (number  of  phases, 
back  voltage  (Ea),  and  synchronous  inductance  (Ls)),  rectifier  values,  converter 
parameters,  and  filtering  components 

•  Simulate  and  analyze  the  PSIM  model  verifying  the  system  meets  the  required 
specifications  for  power  and  output  voltage 

•  Determine  the  generator  line  harmonic  current  magnitudes  from  the  PSIM  model 

•  Using  the  harmonic  current  values,  calculate  the  rotor  losses  caused  by  the  winding  time 
and  space  harmonics  as  detailed  in  section  5.3.1 

•  Calculate  the  switching  losses  and  conduction  losses  of  the  power  conversion  module  per 
section  4.4 

•  Determine  the  overall  losses  and  efficiency  of  the  system 


94 


•  Calculate  the  masses  and  volumes  of  the  generator  and  power  conversion  module  using 

the  methodology  of  sections  3.4.2  and  4.5 

5.5  Optimization 

Optimizing  the  overall  power  generation  module  involves  designing  numerous  PM 
generators  and  associated  power  conversion  electronics  to  achieve  an  ideal  result.  Using  the 
design  procedure  from  section  5.4,  machine  parameters  are  varied  producing  different  generator 
designs.  For  each  unique  generator,  a  power  conversion  module  is  devised  so  that  its  values  are 
matched  with  the  output  of  the  generator.  This  results  in  a  wide  range  of  power  generation 
modules. 

The  main  machine  parameters  that  are  varied  are  the  number  of  poles,  number  of  phases, 
and  the  retaining  sleeve  material  as  shown  in  Figure  37.  This  results  in  80  different  machines 
and  associated  power  conversion  modules.  For  each  of  these  generators,  the  magnet  thickness 
(hm),  air  gap  (g),  magnet  angle  (0m),  and  stack  length  (Lst)  are  iterated  to  ensure  the  sizes,  current 
density,  magnetic  saturation,  sleeve  thickness/stresses,  output  waveforms,  and  THD  are  all 
satisfactory. 

Rotational  speed  (13,000  RPM)  and  magnet  skew  angle  (10°)  are  held  constant.  The 
number  of  slots  (N)  is  set  to  36  for  the  3  and  5-phase  machines  and  to  72  for  the  7,  9,  and  1 1- 
phase  machines  to  properly  fit  the  windings  and  ensure  a  reasonably  sinusoidal  back  voltage 
waveform.  Other  parameters  including  the  number  of  slots/pole/phase  (m),  number  of  armature 
turns  (Na),  back  voltage  (Ea),  synchronous  inductance  (Ls),  electrical  frequency  (f),  breadth 
factor  (kb),  pitch  factor  (kp),  skew  factor  (ks),  magnet  factor  (kg),  and  air  gap  flux  density  (Bg)  all 
change  as  a  result  of  the  machine  variations. 
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Figure  37:  Machine  Optimization  Parameters 

Each  machine  is  optimized  to  have  minimal  back  voltage,  size,  and  mass  while  ensuring 
it  remains  within  the  current  density  and  magnetic  saturation  limits.  The  power  conversion 
modules  are  then  developed  matching  component  ratings  to  those  of  the  PM  generators.  The 
power  electronics  are  designed  to  produce  maximum  power  factor  and  minimal  line  harmonics 
as  discussed  in  section  4.3.  This  yields  the  lowest  possible  losses  in  the  PM  generator  thus 
maximizing  efficiency  and  minimizing  rotor  heat  generation. 
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Chapter  6  Results  and  Analysis 


6.1  General 

As  discussed  in  section  5.5,  the  main  machine  parameters  that  are  varied  are  the  number 
of  poles,  number  of  phases,  and  the  retaining  sleeve  material  and  resulting  in  80  different 
machines  and  associated  power  conversion  modules.  These  PM  machines  and  power 
conversions  modules  are  developed  with  varying  results.  Appendix  K  and  Appendix  L  contain 
the  detailed  specifications  for  each  of  the  generators  and  power  conversion  modules.  Appendix 
M  and  Appendix  N  have  the  in-depth  losses  and  weights  for  each  design.  Appendix  O  lists  the 
mass  and  volume  results  for  the  rectifiers  and  Appendix  P  contains  the  mass  and  volume  data  for 
the  converters.  A  summary  of  the  range  of  values  for  the  designs  is  contained  in  Table  15  and 
Table  16. 


Table  15:  General  Module  Specifications 


Back 

Voltage  Ea 
(RMS) 

Synch 
Inductance 
Ls  (mH) 

Gen  Volt 
THD 

Stack 

Length 

Lst  (m) 

Overall 

Machine 

Length 

L(m) 

Machine 
Diameter  D 
(m) 

Air  Gap  Flux 
Density 
Bg(T) 

PCM  Vol 
(m3) 

Minimum 

1,136.0 

0.020 

3.12% 

0.850 

1.732 

0.834 

0.61 

6.54 

Maximum 

3,380.0 

0.348 

23.05% 

1.080 

2.057 

0.900 

0.77 

7.96 

1,977.0 

0.110 

11.61% 

0.943 

1.861 

0.866 

0.68 

7.14 

Table  16:  General  Module  Results 


Machine 

Weight 

_  Jkg)_ 

BIS 

Total 

Weight 

(kg) 

Machine 

Losses 

(kW) 

WB2LM 

1 

■m 

mmm 

1 

1 Wm 

Line 

Current 

THD 

Minimum 

1,388.6 

2,302.7 

esebu 

241.7 

1,688.3 

2,126.7 

6.80% 

Maximum 

2,302.9 

3,954.7 

6,101.8 

1 ,259.8 

3,293.7 

3,535.4 

24.19% 

1,754.6 

3,025.5 

4,780.2 

463.8 

2,124.0 

2,587.8 

14.56%  [ 

In  general,  the  PCM  losses  are  too  high  due  to  the  IGBT  switching  losses.  This  is 
because  high-frequency  hard  switching  is  used  in  the  power  conversion  module  as  opposed  to 
soft  switching  or  a  lower  frequency.  Soft  switching  involves  using  snubber  circuits  or  resonant 
converters  to  minimize  device  switching  losses  by  constraining  the  switching  of  the  power 
devices  to  time  intervals  when  the  voltage  across  the  device  or  the  current  through  it  is  nearly 


97 


zero.  Soft  switching  can  reduce  switching  losses  as  much  as  50-80%  [67],  [68].  In  addition  to 
lowering  switching  losses,  soft  switching  reduces  EMI,  permits  higher  switching  frequencies, 
and  achieves  higher  efficiencies.  Analyzing  different  soft-switching  topologies  is  beyond  the 
scope  of  this  paper  but  the  effects  soft  switching  has  on  the  PCM  are  studied.  Using  a  lower 
switching  frequency  is  not  examined  because  this  would  produce  a  larger  converter  since  the 
passive  components  would  have  to  be  bigger. 

In  order  to  select  a  final  design,  the  different  variants  are  examined  to  ascertain  which 
one  is  optimal.  Parameters  such  as  weights,  machine  losses,  PCM  losses,  and  THD  are  used  to 
help  facilitate  the  selection  process. 

6.2  Number  of  Phases 

To  facilitate  studying  the  effects  the  number  of  phases  has  on  the  power  generation 
module,  several  parameters  are  held  constant.  The  retaining  sleeve  material  and  number  of  poles 
are  arbitrarily  set  to  inconel  and  12.  This  does  not  affect  the  analysis  because  the  trends  related 
to  varying  the  number  of  phases  are  consistent  across  all  80  power  modules  as  indicated  in  the 
detailed  results  contained  in  the  appendices. 

As  the  number  of  phases  varies  from  3  to  1 1 ,  the  weights  of  the  PM  generator  and  the 
power  conversion  module  change,  but  only  slightly.  Figure  38  shows  this  and  there  is  a 
maximum  variation  of  approximately  242  kg,  or  5.7%.  The  1 1-phase  machine  produces  the 
lowest  weight  power  module  but  only  by  a  small  amount. 


Weights  vs.  Phase  Number 
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Figure  38:  Weights  vs.  Number  of  Phases 
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The  machine  losses  fluctuate  significantly  as  the  phase  number  changes  with  the  3-phase 
machine  having  48.8%  of  the  losses  of  the  1 1  phase  machine  (see  Figure  39).  The  core, 
windage,  and  armature  losses  are  approximately  constant  for  the  different  phase  generators,  but 
the  rotor  losses  vary  considerably.  This  disparity  is  caused  predominately  by  the  AC  line  current 
harmonics  from  the  rectifier  being  much  greater  in  the  higher  phase  machines.  Their  triple-n 
harmonics  are  higher  order  thus  allowing  lower  order  harmonics  to  have  a  greater  effect, 
specifically  the  3rd  harmonic. 

The  higher  rotor  losses  are  also  caused  by  the  geometry  of  the  machine.  As  the  phase 
number  increases,  the  number  of  slots/pole/phase  (m)  decreases  which  causes  the  winding 
factors  to  have  less  of  an  effect  in  reducing  higher  harmonic  effects.  The  exception  is  the  jump 
from  5-phases  to  7-phases  because  the  number  of  slots  is  doubled  from  36  to  72  so  that  m 
increases  initially.  Overall,  the  3-phase  generator  has  the  lowest  rotor  and  total  machine  losses. 


Machine  Losses  vs.  Phase  Number 
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Figure  39:  Machine  Losses  vs.  Number  of  Phases 
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PCM  Losses  vs.  Phase  Number 
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Figure  40:  PCM  Losses  vs.  Number  of  Phases 
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The  PCM  losses  also  differ  greatly  depending  on  the  number  of  phases  (Figure  40).  The 
1 1 -phase  machine  has  the  lowest  total  losses  due  to  it  having  the  lowest  IGBT  switching  losses. 
This  is  because  it  has  a  smaller  rectified  voltage  which  reduces  the  stresses  and  switching  losses 
in  the  IGBTs  (see  section  4.4).  However,  if  soft  switching  is  implemented,  the  PCM  losses  are 
similar  for  the  different  phase  machines  varying  by  a  maximum  of  63.5  kW  instead  of  the  843 
kW  associated  with  hard  switching.  Therefore,  the  lower  phase  machines  become  more 
attractive  if  soft  switching  is  utilized. 

Lastly,  as  the  number  of  phases  increases,  Figure  41  shows  that  in  general,  the  generator 
voltage  and  line  current  THD  increase.  This  is  caused  by  the  same  effects  as  with  the  rotor 
losses,  the  slots/pole/phase  and  the  AC  line  harmonics  from  the  rectifier.  As  mentioned  in 
section  5.1,  the  specification  from  MIL-STD-1399  for  voltage  THD  is  normally  5%  so  using 
these  standards  as  guidelines,  only  the  3-phase  and  the  7-phase  power  modules  are  satisfactory. 
The  line  current  THD  should  be  made  as  low  as  possible  because  this  directly  affects  the  amount 
of  losses  and  heating  in  the  rotor.  The  typical  limit  for  line  current  THD  is  10%  making  the  3- 
phase  generator  the  only  one  with  acceptable  THD  performance  [69]. 
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THD  vs.  Phase  Number 
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Figure  41:  THD  vs.  Number  of  Phases 


6.3  Retaining  Material 

Similar  to  the  analysis  of  the  number  of  phases,  several  parameters  are  held  constant 
when  examining  the  various  retaining  sleeve  materials.  Specifically,  the  number  of  phases  is 
arbitrarily  set  to  3  and  the  number  of  poles  is  fixed  at  6.  This  does  not  affect  the  analysis 
because  the  trends  related  to  varying  the  retaining  material  are  consistent  across  the  continuum  of 
power  modules  as  shown  in  the  detailed  results  contained  in  the  appendices. 

With  respect  to  the  weights  of  the  PM  generator  and  the  PCM,  the  type  of  retaining 
sleeve  has  little  effect.  Using  a  stainless  steel  retaining  can  produces  the  largest  power  module 
but  it  is  not  appreciably  larger  than  the  others.  When  the  retaining  sleeve  is  constructed  from 
titanium,  inconel,  or  carbon  fiber,  the  difference  in  total  weight  is  less  than  0. 1 1%.  Figure  42 
shows  the  effects  the  retaining  sleeve  material  has  on  the  weights  of  the  power  generation 
module. 
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Figure  42:  Weights  vs.  Retaining  Sleeve  Material 


Figure  43:  Machine  Losses  vs.  Retaining  Sleeve  Material 

The  type  of  retaining  sleeve  material  does  affect  the  losses  in  the  PM  generator.  The 
core,  windage,  and  armature  losses  vary  only  slightly  but  the  rotor  losses  fluctuate  a  great  deal. 
Inconel  has  the  lowest  rotor  losses  (77.3  kW)  while  stainless  steel  has  the  highest  (107.1  kW). 
The  carbon  fiber  has  the  smallest  losses  in  the  actual  can  but  larger  losses  in  the  magnets.  This  is 
because  the  carbon  fiber  provides  no  shielding  for  the  magnets  from  flux  variations  and  winding 


current  harmonics.  Overall,  the  inconel  sleeve  produces  the  lowest  rotor  and  total  machine 
losses  as  shown  in  Figure  43. 

The  PCM  losses  are  not  influenced  considerably  by  the  kind  of  retaining  can  material 
(Figure  44).  Depending  on  the  type,  the  losses  vary  by  about  157.2  kW,  or  4.2%.  When  soft 
switching  is  utilized,  the  losses  fluctuate  by  only  23.7  kW  or  2.5%.  Therefore,  the  PCM  losses 
are  not  a  major  factor  in  selecting  the  retaining  can  material. 
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Figure  44:  PCM  Losses  vs.  Retaining  Sleeve  Material 

Lastly,  the  selection  of  a  retaining  sleeve  material  is  not  a  factor  in  the  generator  voltage 
THD  or  the  AC  line  current  THD.  The  voltage  THD  varies  by  only  0.5%  and  the  line  current  by 
less  than  0.15%  depending  on  the  type  of  material  as  indicated  in  Figure  45.  It  is  evident, 
however,  that  selected  a  3-phase,  6-pole,  36-slot  machine  would  not  be  wise  since  the  voltage 
THD  is  above  the  5%  limit  in  all  cases. 
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THD  vs.  Retaining  Sleeve  Material 
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Figure  45:  THD  vs.  Retaining  Sleeve  Material 
6.4  Number  of  Poles 

The  number  of  poles  in  a  machine  is  an  important  factor  because  it  affects  numerous 
parameters  including  the  electrical  frequency,  the  magnet  pole  pitch,  and  the  magnetic  air  gap 
flux  density.  To  facilitate  studying  the  effects  of  the  number  of  poles,  several  parameters  are 
held  constant.  The  retaining  sleeve  material  and  number  of  phases  are  arbitrarily  set  to  inconel 
and  3.  This  does  not  affect  the  analysis  because  the  trends  related  to  varying  the  number  of  poles 
are  consistent  across  the  range  of  machines  as  indicated  in  the  detailed  results  in  the  appendices. 

First,  the  weights  of  the  generator  and  PCM  are  examined  to  ascertain  the  effects  of  the 
number  of  poles.  As  Figure  46  indicates,  as  the  number  of  poles  increases  the  weights  of  both 
the  machine  and  the  PCM  decrease.  The  generator  weight  goes  down  because  assuming  the 
same  flux  density  and  circumferential  arc,  a  greater  number  of  poles  produces  the  same  radial 
flux  but  requires  less  stator  core  back  iron  (see  Figure  5).  The  PCM  weight  decreases  because 
the  higher  frequencies  and  lower  rectified  voltage  associated  with  the  increased  pole  numbers 
reduce  the  size  and  energy  requirements  of  the  passive  components,  especially  the  capacitors. 
Therefore,  the  total  power  generation  module  weight  decreases  as  the  number  of  poles  increases 
making  the  12-pole  machine  the  most  attractive  in  terms  of  weight. 


104 


Weights  vs.  Pole  Number 
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Figure  46:  Weights  vs.  Number  of  Poles 

Next,  looking  at  Figure  47,  as  the  number  of  poles  goes  up,  the  losses  in  the  machine 
increase,  some  more  than  others.  The  core  losses  increase  because  the  higher  frequencies  cause 
the  eddy  currents  and  hysteresis  to  have  a  greater  effect.  However,  the  range  in  core  losses  is 
only  7  kW  so  this  increase  is  not  a  major  contributor.  The  windage  losses  are  not  significantly 
affected  by  the  number  of  poles  differing  by  only  3  kW  across  the  pole  variation. 

Machine  Losses  vs.  Pole  Number 
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Figure  47:  Machine  Losses  vs.  Number  of  Poles 

The  armature  losses  do  fluctuate  some  as  the  pole  number  changes  with  the  6-pole 
machine  having  the  lowest  (1 17.6  kW)  and  the  12-pole  machine  having  the  highest  (137.9  kW). 
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This  is  due  to  the  higher  pole  machines  having  higher  stator  current  densities  as  the  machine 
performance  is  optimized.  The  rotor  losses  also  increase  as  the  number  of  poles  gets  higher 
mainly  because  of  the  winding  time  and  space  harmonics.  As  the  pole  number  increases,  the 
number  of  slots/pole/phase  (m)  decreases  which  causes  the  winding  factors  to  have  less  of  an 
influence  in  reducing  higher  harmonic  effects. 

Overall,  all  of  the  different  pole  machines  have  less  than  160  kW  of  rotor  losses  with  the 
6-pole  machine  having  the  lowest.  However,  as  mentioned  in  section  6.3,  the  generator  voltage 
THD  is  more  significant  in  this  machine  causing  the  higher  pole  machines  to  become  viable 
options. 

The  number  of  poles  also  produces  sizeable  effects  on  the  PCM  losses.  As  the  pole 
number  increases,  the  PCM  losses  decrease  due  to  much  lower  IGBT  switching  losses.  This  is 
because  the  higher  pole  machines  have  lower  rectified  voltages  which  reduce  the  stresses  and 
switching  losses  in  the  IGBTs  (see  section  4.4).  The  difference  between  having  6-poles  and  12- 
poles  is  754.3  kW  in  PCM  losses,  a  22.9%  reduction.  If  soft  switching  is  implemented,  the  PCM 
losses  are  still  lower  for  the  higher  pole  machines,  but  the  difference  is  now  98.1  kW,  a  10.3% 
reduction.  Overall,  as  shown  in  Figure  48,  the  12-pole  machine  is  the  best  option  when 
considering  PCM  losses. 
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Figure  48:  PCM  Losses  vs.  Number  of  Poles 

The  generator  voltage  THD  and  the  AC  line  current  THD  are  both  affected  by  the 
number  of  poles.  Examining  Figure  49,  the  10-pole  and  12-pole  generators  are  the  only  variants 
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that  produce  voltage  waveforms  that  have  acceptable  THD  levels,  3.74%  and  4.95% 
respectively.  The  AC  line  current  THD  values  are  satisfactory  for  all  of  the  pole  numbers  with 
the  12-pole  machine  having  the  lowest  value  (6.82%).  If  the  THD  is  examined  in  conjunction 
with  the  machine  losses,  the  12-pole  machine  is  the  best  option.  It  has  the  top  overall  THD 
performance,  has  25.8  kW  lower  rotor  losses  compared  to  the  10-pole  generator,  and  has  14.6 
kW  lower  total  machine  losses  compared  to  the  10-pole  design. 
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Figure  49:  THD  vs.  Number  of  Poles 
6.5  Final  Power  Module 

6.5.1  PM  Generator 

Based  on  the  previous  analysis  and  the  goal  of  producing  an  optimal  power  generation 
module,  a  final  design  is  selected.  The  sizing  codes  and  analysis  from  sections  3.4.1,  3.4.2, 

3.4.3,  and  3.4.4  are  also  utilized  in  this  process.  The  final  PM  generator  is  a  3-phase,  12-pole 
machine  with  an  inconel  retaining  sleeve.  The  final  design  is  not  a  fully  optimized  in  terms  of  all 
parameters  but  represents  the  best  selection  using  the  methodology  of  this  paper.  Detailed 
characteristics  for  the  PM  generator  are  listed  in  Table  17. 
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Table  17:  PM  Generator  Final  Design  Parameters 


B  . . 

1 

Total  Machine  Diameter 

0.848  m 

Total  Machine  Length 

1.868  m 

Rotor  Radius 

0.147  m 

Active  Length 

1.000  m 

Slot  Average  Width 

16.973  mm 

Slot  Height 

25.000  mm 

Back  Iron  Thickness 

17.150  mm 

Tooth  Width 

15.882  mm 

1  'V ;; v ;  'i' 

pv: 

Power  Rating 

16,000  kW 

Speed 

Va(RMS) 

1,739  V 

Phase  Current 

3,169.3  A 

Ea  (RMS) 

2,380  V 

Armature  Resistance 

0.00457  n 

Synchronous  Reactance 

0.377  Q 

Synchronous  Inductance 

0.046  mH 

lEMEsssmsmm 

■fMSMLfTSRU 

Tip  Speed 

200  m/s 

Efficiency 

0.979 

Power  Factor 

0.960 

Phases 

3 

Frequency 

1,300  Hz 

1 

f||f 

■  ■ 

Number  of  Slots 

36 

Number  of  Armature  Turns 

12 

Fund  Breadth  Factor 

1.000 

Fund  Pitch  Factor 

0.866 

Tooth  Flux  Density 

1.30  T 

Back  Iron  Flux  Density 

0.93  T 

1.00 

29.00  mm 

Magnet  Angle 

20  deg 

Air  Gap 

5.00  mm 

Pole  Paris 

6 

Magnet  Remanence 

1.20  T 

Air  Gap  Bg 

0.69  T 

Magnet  Factor 

0.961 

Skew  Factor 

0.955 

Ml/ 

Core  Loss 

19.9  kW 

Armature  Loss 

137.9  kW 

Windage  Loss 

34.4  kW 

Rotor  Loss 

158.9  kW 

g|||  *  '  ' 

'  I.  _ : ..  ..  '  -. 

■-  — 

■ 

’.ill 

Sfi 

Core 

292.49  kg 

Shaft 

522.73  kg 

Magnet 

72.59  kg 

Armature 

82.48  kg 

Services 

145.54  kg 

Structure 

557.92  kg 

Total 

1,673.75  kg 

A  diagram  of  the  PM  generator  is  shown  in  Figure  50  with  the  stator  winding  layout 
given  in  Table  18.  The  numbers  in  the  winding  layout  table  refer  to  the  slot  number  locations  of 
Figure  50. 
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Figure  50:  Diagram  of  PM  Generator  Final  Design 


Table  18:  Winding  Layout 
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6.5.2  Power  Electronics  Module 

In  conjunction  with  selecting  the  final  PM  generator  design,  the  associated  power 
electronics  module  is  determined.  Figure  5 1  shows  the  complete  power  module,  including 
individual  component  values. 
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Figure  51:  Power  Module  Final  Design  Diagram 

Along  with  Figure  51,  additional  characteristics  for  the  final  power  module  are  given  in  Table  19 
with  detailed  data  contained  in  the  appendices. 


Table  19:  PCM  Final  Design  Parameters 


I  i(  i ii]|i ! fi i >' pi-  'Hii ■  ■  ii s  i:  1 tf  i  f  m  i  tf i piiiiii "  , 

IGBT  Switching  Losses 

2,103.6  kW 

Diode  Conduction  Losses 

168.2  kW 

IGBT  Conduction  Losses 

267.6  kW 

Total  PCM  Losses 

2,539.4  kW 
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IGBT  Switching  Losses 

420.7  kW 

Diode  Conduction  Losses 

168.2  kW 

IGBT  Conduction  Losses 

267.6  kW 

Total  PCM  Losses 
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6.5.3  Performance/Waveforms 


As  discussed  in  section  5.4,  all  80  power  modules  are  constructed  and  simulated  using 
PSIM.  In  order  to  conserve  space,  only  the  simulation  results  for  the  final  design  are  presented 
here.  The  module  is  simulated  operating  at  the  maximum  load  of  16  MW  and  at  the  minimum 
load  of  100  kW. 

For  the  16  MW  loading,  Figure  52  shows  the  generator  voltage  waveforms.  Since  the 
THD  is  low,  the  curves  fairly  closely  resemble  sine  waves. 


Figure  52:  PM  Generator  Voltage  Waveforms 


Figure  53  through  Figure  55  show  plots  for  various  parameters  for  the  complete  module  loaded 
at  16  MW.  The  output  voltage  fluctuates  approximately  0.09  V  (0.0013%)  which  is  within 
specification,  the  output  current  varies  by  less  than  22  A  (0.0096%),  and  the  AC  line  currents 
have  low  harmonic  content. 


Ill 


Figure  54:  Output  Current  (16  MW) 
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Figure  55:  AC  Line  Current  (16  MW) 


Figure  56  through  Figure  58  display  the  module  performance  at  the  minimum  load  of  100 
kW.  The  output  voltage  ripple  is  1 .4  V,  higher  than  at  full  load  but  still  not  significant.  The 
output  current  also  has  some  variation  (about  25  A)  but  this  is  acceptable  at  the  lower  power 
levels.  The  AC  line  current  has  much  higher  harmonic  content  because  of  notches  that  occur  in 
the  waveforms  caused  by  the  DC  link  capacitor  charging  and  discharging  at  the  low  load. 
However,  the  magnitudes  of  the  harmonics  are  insignificant  since  the  magnitude  of  the 
fundamental  is  low  and  therefore  the  rotor  losses  and  heating  are  not  adversely  affected. 
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Figure  56:  Output  Voltage  (100  kW) 
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Figure  57 :  Output  Current  (100  kW) 
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Figure  58:  AC  Line  Current  (100  kW) 


6.6  Comparison 

The  premise  for  designing  and  employing  a  high-speed  PM  generator  and  power 
conversion  module  is  to  offer  a  light  weight,  power  dense  alternative  to  existing  wound  rotor 
machines.  In  Table  1  in  section  1.3.2,  parameters  for  some  typical  wound  rotor  machines  are 
presented.  In  order  to  compare  them  with  the  new  design  power  module,  these  machines  are 
scaled  to  16  MW  using  the  general  relationship  of  Eqn  6-1. 

4 

Pwr  Vol 3 


Eqn  6-1 


Table  20  then  provides  an  evaluation  of  the  different  machines  and  power  modules. 


Table  20:  Comparison  of  Machines/Modules 


[  Machine 

Phag|| 

fMW) 

^ijg 

tm 

H(m)ht 

Volume 

Tmt)1 

i 

3 

16 

3,600 

4.4 

3.7 

3.2 

52.0 

40.8 

2 

3 

16 

3,600 

5.0 

3.6 

3.1 

54.3 

34.3 

3 

3 

16 

3,600 

5.5 

3.2 

4.0 

69.0 

47.3 

4 

3 

16 

3,600 

4.6 

2.8 

3.7 

47.6 

41.0 

PM  Gen 

3 

16 

13,000 

1.87 

0.85 

0.85 

1.4 

1.67 

PCM 

N/A 

16 

N/A 

2.0 

1.66 

2.0 

6.6 

2.32 

Total  PM 

3 

16 

13,000 

3.87 

Varies 

Varies 

8.0 

3.99 
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Examining  Table  20,  the  new  design  offers  significant  reductions  in  both  weight  and 
volume.  Specifically,  the  new  total  power  module  is  1 1.6  -  16.8%  the  volume  and  8.5  -  1 1.7% 
the  weight  of  the  other  machines.  This  translates  to  about  a  7x  reduction  in  volume  and  a  lOx 
reduction  in  weight.  These  reductions  can  provide  flexibility  to  naval  architects  since  power, 
weight,  and  volume  are  integral  parts  of  the  design  and  construction  processes. 
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Chapter  7  Conclusions  and  Recommendations 


7.1  Design  Lessons  Learned 

Numerous  lessons  are  learned  from  both  the  generator  and  power  electronics  design 
processes.  The  entire  16  MW  power  module  is  a  complex  system  with  many  interrelated 
components  that  affect  the  overall  performance. 

7.1.1  PM  Generator 

The  final  generator  is  the  product  of  examining  multiple  parameters  to  determine  an 
optimal  design.  Various  lessons  are  learned  relating  to  different  portions  of  the  machine 
including  the  number  of  slots,  the  number  of  poles,  the  number  of  phases,  the  magnetic 
geometry,  the  windings  and  output  voltage,  and  optimization  methods. 

7. 1.1.1  Number  of  Slots 

Initially,  it  was  thought  that  the  number  of  slots  could  be  held  constant  over  the  range  of 
machine  designs.  However,  this  is  not  possible  because  as  the  number  of  phases  increases, 
additional  slots  are  necessary  to  fit  the  greater  number  of  windings.  One  of  the  most  important 
parameters  that  can  be  used  to  help  reduce  the  harmonics  in  a  rotating  machine  is  the  number  of 
slots/pole/phase  (m).  Therefore,  the  number  of  slots  can  be  used  to  help  optimize  this  parameter 
and  improve  the  shape  of  the  output  voltage  waveform  and  THD. 

The  number  of  slots  also  influences  the  electrical  loading  of  the  machine.  The  number  of 
slots  directly  affects  the  slot  size  which  helps  determine  the  current  density  in  the  machine.  This 
is  usually  one  of  the  limits  on  the  machine’s  performance  because  of  cooling.  Overall,  the 
number  of  slots  is  an  important  parameter  that  must  be  carefully  selecting  when  designing  a 
machine. 

7. 1.1. 2  Number  of  Poles 

Traditionally,  the  number  of  poles  is  driven  by  the  system  frequency  and  prime  mover 
speed.  With  a  high-speed  generator  rectified  to  a  DC  distribution  system,  the  number  of  poles 
becomes  an  important  design  variable.  Therefore,  throughout  the  design  progression,  careful 
consideration  should  be  given  to  the  number  of  poles.  This  proved  to  be  true  throughout  this 
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paper  since  the  number  of  poles  affects  the  electrical,  magnetic,  and  structural  performance, 
including  the  electrical  frequency,  the  voltage  waveforms,  the  magnetic  flux,  the  magnet  volume, 
the  magnet  pole  pitch,  the  air  gap  size,  and  the  stator  back  iron  thickness. 

It  is  recommended  that  initially  the  designer  select  a  range  of  acceptable  output  electrical 
frequencies  which  puts  bounds  on  the  allowable  number  of  poles.  From  there,  the  number  of 
poles  can  be  varied  to  help  determine  the  size,  waveforms,  and  performance  of  the  generator.  In 
general,  if  the  rotational  speed  is  held  constant  as  the  number  of  poles  increases: 

•  The  number  of  slots/pole/phase  (m)  changes  affecting  the  output  waveforms  and  THD 

•  The  weight  of  the  machine  decreases 

•  The  electrical  frequency  increases 

7 . 1 . 1 . 3  Number  of  Phases 

An  odd  number  of  phases  is  usually  the  best  option  for  rotating  machines  because  it 
produces  balanced  operation  with  limited  mechanical  fluctuations.  Initially,  it  was  assumed  that 
a  higher  number  of  phases  (greater  than  3)  would  produce  a  more  optimal  power  module.  It  is 
discovered  that  this  is  not  necessarily  the  case.  As  the  number  of  phases  increases,  m  decreases 
which  can  adversely  affect  the  voltage  THD.  Also,  the  AC  line  current  harmonics  are  more 
substantial  in  increased  phase  machines  because  the  triple-n  harmonics  are  higher  order  as  the 
phase  number  increases.  However,  a  big  advantage  of  having  a  higher  number  of  phases  is  that 
it  reduces  the  requirements  for  the  PCM  because  of  lower  rectifier  voltage  ripple  and  lower 
phase  voltages  and  currents. 

It  is  concluded  that  the  benefits  of  3-phases  and  a  higher  number  of  phases  can  be 
simultaneously  realized  through  a  modified  design.  Multiple  3-phase  sets  of  windings  can  be 
located  equidistant  around  the  stator  and  then  connected  to  the  PCM  to  achieve  the  desired 
results.  An  example  of  this  would  be  three  3-phase  sets  of  windings,  one  set  at  0°,  120°,  and 
240°;  one  at  45°,  165°,  and  285°;  and  one  at  90°,  210°,  and  330°. 

7 . 1 . 1 .4  Magnetic  Geometry 

Designing  the  magnetic  geometry  in  a  PM  generator  is  an  intricate  process.  Careful 
consideration  should  be  given  to  the  type  of  retaining  material,  the  magnet  material,  and  the 
desired  voltage.  It  is  determined  that  the  magnet  height  and  air  gap  length  greatly  affect  the 
overall  machine  performance,  and  if  they  are  not  properly  selected,  the  generator  does  not 
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operated  properly.  A  typical  starting  point  in  a  design  is  to  have  the  magnet  height  be 
approximately  5  -  lOx  the  length  of  the  air  gap  [17],  [18],  [30].  This  ensures  that  the  magnetic 
flux  is  adequate  while  not  saturating  the  core  material. 

7. 1 . 1 .5  Windings  and  Output  Voltage 

As  discussed  in  section  3.2.7,  wye  series  connected  windings  are  used  in  this  paper.  This 
produces  a  generator  output  voltage  that  is  too  high  compared  to  the  desired  PCM  output  voltage 
of  700  V.  It  also  unnecessarily  stresses  the  PCM  devices  and  produces  a  lower  than  optimal  duty 
cycle.  It  is  concluded  that  connecting  some  of  the  generator  windings  in  parallel  or  using  a  lower 
number  of  turns  to  produce  a  lower  output  voltage  might  be  viable  options  to  alleviate  this 
problem.  However,  care  must  be  taken  to  account  for  possible  circulating  currents  in  the  parallel 
windings  in  the  event  of  a  fault. 

7. 1 . 1 .6  Optimization  Methods 

When  designing  the  PM  generator,  the  goal  was  to  produce  an  optimal  design  that  met 
the  required  performance  criteria  while  minimizing  the  rotor  losses,  the  size,  and  the  weight  of 
the  machine.  Optimizing  the  design  involves  trying  to  simultaneously  maximize  both  the 
electrical  loading  (current  density)  and  magnetic  loading  (air  gap  flux  density).  It  is  concluded 
that  it  is  extremely  difficult  to  do  both  at  the  same  time  since  there  are  competing  factors  such  as 
saturation  limits,  cooling  limits,  and  rotational  stress  limits. 

In  this  paper,  each  of  the  80  machines  is  optimized  by  trying  to  maximize  the  current 
density  (within  the  limits)  while  ensuring  an  adequate  air  gap  flux  density.  This  is  done  because 
it  seems  to  produce  the  smallest,  high-performance  machines.  It  is  therefore  recommended  in 
any  PM  generator  design  process  that  careful  consideration  be  given  to  both  the  electrical  and 
magnetic  loading. 

7. 1 .2  Power  Electronics 

The  power  electronics  module  (PEM)  designed  in  conjunction  with  the  PM  generator 
uses  one  of  many  possible  topologies.  The  combination  of  the  passive  rectifier  and  buck 
converter  is  a  stable,  well-proven  design.  It  is  optimized  to  match  the  generator  by  selecting  the 
passive  components  properly.  However,  the  design  does  have  limitations  because  it  uses  hard- 
switching  and  large  passive  components. 
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7. 1.2.1  Rectification 

Passive  rectification  is  utilized  in  the  PEM  because  of  its  simplicity  and  potential  for 
lower  input  harmonics.  It  is  concluded  that  the  rectification  scheme  is  satisfactory  and  couples 
well  with  the  PM  generator.  However,  it  is  recognized  that  a  more  detailed  analysis  of  the 
rectification  scheme  could  produce  a  better  topology. 

7. 1.2.2  Conversion  and  Switching 

As  discussed  in  Chapter  6,  one  of  the  main  problems  with  the  power  electronics  module 
is  that  the  losses  are  too  high  due  to  the  IGBT  switching  losses.  This  is  because  hard  switching 
is  used  in  the  power  conversion  module  as  opposed  to  soft  switching.  Soft  switching  minimizes 
device  switching  losses,  reduces  EMI,  permits  higher  switching  frequencies,  and  achieves  higher 
efficiencies. 

Analyzing  different  soft-switching  topologies  is  beyond  the  scope  of  this  paper  but  the 
effects  soft  switching  has  on  the  PCM  are  studied.  It  is  concluded  that  if  soft  switching  is 
implemented  in  the  PCM,  the  efficiency  increases  by  10  -  20%  with  the  losses  reduced  by 
overall  1  MW.  Given  the  power  requirements  of  the  power  generation  module,  soft  switching 
should  be  employed.  It  is  therefore  recommended  that  in  any  high-power  system  that  utilizes 
power  conversion  modules,  soft  switching  should  be  examined  to  ascertain  its  usefulness. 

7. 1.2.3  Control 

One  control  scheme  for  providing  voltage  regulation  for  the  power  generation  module  is 
discussed  in  section  4.3.4.  This  control  methodology  utilizing  duty  ratio  control  with  a  PID 
controller  is  relatively  simple  to  implement.  This  ease  of  control  is  possible  because  the  buck 
converter  is  a  straightforward  topology.  It  is  concluded  that  if  a  more  complex,  soft-switching 
design  is  implemented,  more  extensive  control  schemes  would  be  required. 

Even  though  it  is  not  discussed  in  detail,  the  PM  generator  provides  for  easy  control 
because  it  is  a  synchronous  machine.  This  is  a  distinct  advantage  over  an  induction  generation 
which  requires  extremely  complicated  control  algorithms. 

7.2  Power  Generation  Module 

Overall,  80  different  16  MW  machines  and  power  conversion  modules  are  sized, 
designed,  and  analyzed  with  a  final  design  selected.  The  permanent  magnet  generation  module. 
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consisting  of  the  generator  and  associated  power  electronics,  provides  an  excellent  alternative  to 
traditional  wound  rotor  machines  for  naval  applications.  It  is  concluded  that  the  new  design 
offers  significant  reductions  in  both  weight  and  volume.  Specifically,  it  is  estimated  that  the  PM 
generation  module  has  a  7x  reduction  in  volume  and  a  lOx  reduction  in  weight  compared  to 
similarly  rated  wound  rotor  systems.  These  reductions  can  provide  flexibility  to  naval  architects 
since  power,  weight,  and  volume  are  integral  parts  of  the  design  and  construction  processes. 
However,  further  study  is  necessary  to  verify  the  16  MW  PM  generation  module’s  thermal, 
structural,  and  mechanical  performance. 

7.3  Recommendations/Further  Study 

This  paper  discussed  the  electrical  and  magnetic  design  and  analysis  of  a  permanent 
magnet  generation  module  including  sizing,  detailed  machine  design  and  analysis,  analytical 
models  for  different  rotor  loss  mechanisms  as  well  as  with  other  machine  losses,  power 
electronics  and  conversion  for  connecting  the  high-speed  generator  to  a  DC  distribution  system, 
and  in  depth  simulation  of  the  complete  system.  Numerous  design  issues  are  addressed  and 
several  issues  are  raised  about  the  potential  improvements  a  PM  generation  system  can  offer.  A 
proposed  PM  generation  module  design  is  presented  along  with  a  detailed  design  methodology. 

However,  this  paper  does  not  full  cover  all  aspects  associated  with  the  design  and 
analysis  of  a  PM  generation  module.  In  order  to  fully  ascertain  whether  the  proposed  design  is 
technically  sound,  further  study  is  required  in  both  the  machine  and  power  electronics  areas. 
Specifically,  with  the  respect  to  the  generator,  the  following  is  recommended: 

•  Using  the  calculated  losses,  perform  a  detailed  thermal  analysis  of  the  machine  including 
development  of  the  required  cooling  systems 

•  Conduct  an  in-depth  structural  and  mechanical  study  of  the  PM  generator,  including 
vibration  and  acoustic  analysis 

For  the  power  electronics  portions,  the  further  examination  is  recommended,  including: 

•  Examine  different  rectification  schemes  to  see  what  effects  they  would  have  on  the 
whole  system 

•  Explore  soft  switching  topologies  and  perform  detailed  analysis  on  how  to  implement 
them  to  ensure  acceptable  switching  losses  while  matching  performance  with  the 
generator 
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For  the  complete  PM  generation  system,  the  following  is  recommended: 

•  Using  the  developed  PM  generator  and  power  electronics  module,  implement  a  control 
scheme  and  conduct  in-depth  analysis  with  the  system  connected  to  an  entire  distribution 
system  under  changing  load  conditions 
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Rotor  Tip  Speed  (m/s) 

Terminal  Voltage 

Retaining  Can  Loss  per  Unit  Area 

Slot  Depression  Width 

Slot  Average  Width 

Slot  Bottom  Width 

Slot  Top  Width  (closest  to  rotor) 

Tooth  Width 

Synchronous  Reactance 

Surface  impedance 


126 


Acknowledgements 

The  author  would  like  to  sincerely  thank  Dr.  James  Kirtley  for  all  his  outstanding  help  and 
guidance.  He  consistently  provided  an  enormous  amount  of  assistance  and  tolerated  numerous 
questions  and  problems,  some  of  which  were  probably  trivial  to  him.  Without  his  continual  help, 
this  thesis  could  not  have  been  completed. 

The  author  would  also  like  to  thank  the  following  people: 

•  Dr.  Timothy  McCoy,  for  providing  valuable  assistance,  reference  material,  and 
instruction  throughout  the  thesis  process. 

•  Mr.  David  Clayton  (NAVSEA  05Z)  for  providing  background  electrical  information. 

•  Mr.  Dana  Delisle  (NAVSEA  05Z)  for  providing  background  information  on  different 
types  of  generators. 

•  Mr.  Joseph  Borraccini  (NSWCCD  [Philadelphia])  for  providing  valuable  information  on 
DC-DC  converters. 


127 


Page  Intentionally  Left  Blank 


128 


List  of  References 


[1]  T.  F.  Chan,  L.  L.  Lai,  &  Lie-Tong  Yan,  “Performance  of  a  Three-Phase  AC  Generator 
with  Inset  NdFeB  Permanent  Magnet  Rotor,”  IEEE  Transactions  on  Energy  Conversion.  Vol.  19, 
No.  1,  March  2004. 

[2]  Surong  Huang,  Jian  Luo,  Franco  Leonardi,  &  Thomas  A.  Lipo,  “A  General  Approach  to 
Sizing  and  Power  Density  Equations  for  Comparison  of  Electric  Machines,”  IEEE  Transactions 
on  Industry  Applications,  Vol.  14,  No.  1,  Jan/Feb  1998. 

[3]  Anteon  Corporation,  “Future  Shipboard  Power  System  Architecture  Analysis  Interim 
Report:  Shipboard  Electrical  Distribution  System  Characteristics,”  December  2001. 

[4]  LCDR  Timothy  McCoy,  Matthew  Stauffer,  &  Edward  Harvey,  “Full  Scale  Land  Based 
Testing  of  the  U.S.  Navy’s  Integrated  Power  System  (IPS),”  Naval  Sea  System  Command, 

March  2000. 

[5]  GMS,  Inc.  &  Draper  Laboratory,  Electric  Power  Applications  to  Naval  Ships.  A 
Technology  Review.  May  1997. 

[6]  Timothy  J.  Doyle,  Howard  O.  Stevens,  &  Henry  Robey,  “An  Historical  Overview  of 
Navy  Electric  Drive.” 

[7]  Ronald  O’Rourke,  Electric  Drive  Propulsion  for  U.S.  Navy  Ships:  Background  and  Issues 
for  Congress.  Congressional  Research  Service,  The  Library  of  Congress,  July  2000. 

[8]  Jason  Pepi  &  Dr.  Peter  Mongeau,  “High  Power  Density  Permanent  Magnet  Generators,” 
DRS  Electric  Power  Technologies,  Inc.,  2004. 

[9]  Wang  Fengxiang,  Zheng  Wenpeng,  Zong  Ming,  &  Wang  Baoguo,  “Design 
Considerations  of  High-Speed  PM  Generators  for  Micro  Turbines,”  IEEE  Transactions.  2002. 

[10]  Thomas  C.  Gillmer  &  Bruce  Johnson,  Introduction  to  Naval  Architecture.  Naval  Institute 
Press,  Dec  1982. 

[11]  Dr.  Timothy  McCoy,  “Contemporary  Naval  Electric  Power  Systems,”  Course  13.412 
Lecture.  Massachusetts  Institute  of  Technology,  Fall  2003. 

[12]  Richard  J.  Newman,  “For  the  Navy,  an  Engine  of  Change,”  U.  S.  News  &  World  Report. 
December  1999. 

[13]  LCDR  Andy  Gish,  LCDR  Bill  Hardman,  LCDR  Steve  Ramsey,  Concept  Design  of  a 
Reduced  Displacement  Fast  Attack  Submarine.  MIT  13  A,  May  2004. 


129 


[14]  Dr.  Makhlouf  Benatmane,  LCDR  Timothy  McCoy,  &  Ms.  Tracey  L.  Cooper,  “Electric 
Propulsion  Full-Scale  Development  in  U.S.  Navy  Surface  Ships,”  October  1998. 

[15]  Science  Applications  International  Corporation,  Superconducting  Generator  Study  Final 
Report.  Dec  2002. 

[16]  James  L.  Kirtley,  “Course  6.685:  Electric  Machines,  Class  Notes  1,”  Massachusetts 
Institute  of  Technology,  2003. 

[17]  Dr.  Duane  Hanselman,  Brushless  Permanent  Magnet  Motor  Design.  The  Writer’s 
Collective,  2003. 

[18]  J.  R.  Hendershot,  Jr.  &  T.  J.  E.  Miller,  Design  of  Brushless  Permanent  Magnet  Motors. 
Magna  Physics  Publishing  and  Clarendon  Press,  1994. 

[19]  James  L.  Kirtley,  “Course  6.685:  Electric  Machines,  Class  Notes  7,”  Massachusetts 
Institute  of  Technology,  2003. 

[20]  E.  Spooner  &  A.  C.  Williamson,  “Direct  Coupled,  Permanent  Magnet  Generators  for 

Wind  Turbine  Applications,”  IEEE  Proceedings.  Electrical  Power  Applications,  Vol.  143,  No.  1, 
Jan  1996.  '”~~ 

[21]  Johannes  J.  H.  Paulides,  Geraint  W.  Jewell,  &  David  Howe,  “An  Evaluation  of 
Alternative  Stator  Lamination  Materials  for  a  High-Speed,  1.5  MW,  Permanent  Magnet 
Generator,”  IEEE  Transactions  on  Magnetics.  Vol.  40,  No.  4,  July  2004. 

[22]  Mark  Rippy,  An  Overview  Guide  for  the  Selection  of  Lamination  Materials,  Proto 
Laminations,  Inc.,  2004. 

[23]  D.  Pavlik,  V.  K.  Garg,  J.  R.  Repp,  &  J.  Weiss,  “A  Finite  Element  Technique  for 
Calculating  the  Magnet  Sizes  and  Inductances  of  Permanent  Magnet  Machines,”  IEEE 
Transactions  of  Energy  Conversion.  Vol.  3,  No.  1,  March  1988. 

[24]  Do  Hyun  Kang,  Paul  Curiae,  Yeon  Ho  Jung,  &  Soo  Jin  Jung,  “Prospects  for 
Magnetization  of  Large  PM  Rotors:  Conclusions  from  a  Development  Case  Study,”  IEEE 
Transactions  on  Energy  Conversion.  Vol.  18,  No.  3,  Sept  2003. 

[25]  Website  for  Magnetic  Component  Engineering,  Inc.,  2004. 

[26]  James  L.  Kirtley,  “Course  6.685:  Electric  Machines,  Class  Notes  5,”  Massachusetts 
Institute  of  Technology,  2003. 

[27]  N.  Bianchi  &  A.  Lorenzoni,  “Permanent  Magnet  Generators  for  Wind  Power  Industry: 

An  Overall  Comparison  with  Traditional  Generators,”  Opportunities  and  Advances  in 
International  Power  Generation.  Conference  Publication  No.  419,  1996. 


130 


[28]  Co  Huynh,  Larry  Hawkins,  Ali  Farahani,  &  Patrick  McMullen,  “Design  and 
Development  of  a  Two-Megawatt,  High-Speed  Permanent  Magnet  Alternator  fro  Shipboard 
Application,”  Calnetix,  Inc.,  2004. 

[29]  Website  www.consult-g2.com/course/chapter7/chapter.html. 

[30]  M.  A.  Rahman  &  G.  R.  Slemon,  “Promising  Applications  of  Neodymium  Iron  Boron  Iron 
Magnets  in  Electrical  Machines,”  IEEE  Transactions  on  Magnetics,  Vol.  No.  5,  Sept  1985. 

[31]  Dr.  Peter  Mongeau,  “High  Torque  Density  Propulsion  Motors,”  DRS  Electric  Power 
Technologies,  Inc.,  2004. 

[32]  O.  Aglen  &  A.  Andersson,  “Thermal  Analysis  of  a  High-Speed  Generator,”  IEEE 
Transactions.  2003. 

[33]  S.  M.  Abu  Sharkh,  M.  R.  Harris,  &  N.  Taghizadeh  Irenji,  “Calculation  of  Rotor  Eddy- 
Current  Loss  in  High-Speed  PM  Alternators,”  IEE  Conference  Publication  No.  444.  1997. 

[34]  D.  K.  Yue,  “Course  13.021:  Marine  Hydrodynamics,  Lecture  Notes  17,”  Massachusetts 
Institute  of  Technology,  2002. 

[35]  H.  Polinder  &  M.  J.  Hoeijmakers,  “Eddy-Current  Losses  in  the  Segmented  Surface- 
Mounted  Magnets  of  a  PM  Machine,”  IEE  Proceedings,  Electrical  Power  Applications,  Vol.  146, 
No.  3,  May  1999. 

[36]  James  L.  Kirtley,  “Course  6.685:  Electric  Machines,  Class  Notes  3,”  Massachusetts 
Institute  of  Technology,  2003. 

[37]  Mihai  Comanescu,  Ali  Keyhani,  &  Min  Dai,  “Design  and  Analysis  of  42V  Permanent 
Magnet  Generator  for  Automotive  Applications,”  IEEE  Transactions  on  Energy  Conversion. 

Vol.  18,  No.  1,2003. 

[38]  A.  Munoz-Garcia  &  D.  W.  Novotny,  “Utilization  of  Third  Harmonic-Induced-Voltages  in 
PM  Generators,”  IEEE,  1996. 

[39]  R.  W.  Erickson,  Fundamentals  of  Power  Electronics,  Prentice-Hall,  1997. 

[40]  Arthur  Kelley  &  William  F.  Yadusky,  “Rectifier  Design  for  Minimum  Line-Current 
Harmonics  and  Maximum  Power  Factor,”  IEEE  Transactions  on  Power  Electronics.  Vol.  7,  No. 

2,  1992. 

[41]  John  G.  Kassakian,  Martin  F.  Schlecht,  &  George  C.  Verghese,  Principles  of  Power 
Electronics.  Addison-Wesley  Publishing  Company,  Inc.,  1991. 

[42]  Ned  Mohan,  Tore  M.  Undeland,  &  William  P.  Robbins,  Power  Electronics:  Converters. 
Applications,  and  Design.  John  Wiley  &  Sons,  Inc.,  1995. 


131 


[43]  Dr.  David  Perreault,  Course  6.334:  Power  Electronics  Notes,  Massachusetts  Institute  of 
Technology,  2004. 

[44]  Website  http://henry.fbe.fh-darmstadt.de/smps_e/ 

[45]  T.  K.  Phelps  &  W.  S.  Tate,  “Optimizing  Passive  Input  Filter  Design,”  Proceedings  of  the 
6th  National  Solid  State  Power  Conversion  Conference.  1979. 

[46]  Website  www.electronicstalk.com,  March  2002. 

[47]  Website  www.pwrx.com. 

[48]  Website  www.pi.hitachi.co/jp/pse/index.html. 

[49]  Website  www.ThinGap.com. 

[50]  Website  www.pulseeng.com. 

[51]  Ming-Jen  Pan,  Barry  Bender,  Roy  Rayne,  &  Michael  Lanagan,  Development  of 
Advanced  Dialectrics  for  High  Energy  Density  Capacitors.  2004. 

[52]  Troy  W.  Barbee,  Gary  W.  Johnson,  High  Energy  Density  Capacitors  for  Power 
Electronic  Applications  using  Nano-Structure  Multilayer  Technology.  Lawrence  Livermore 
National  Laboratory,  1995. 

[53]  Technology  for  the  United  States  Navy  and  Marine  Corps,  2000-2035:  Becoming  a  21st 
Century  Force.  Vol.  2:  Technology.  2004. 

[54]  E.  R.  Brown,  Megawatt  Solid-State  Electronics,  DARPA  Microsystems  Technology 
Office,  1998. 

[55]  Yahya  Shakweh,  “MV  Inverter  Stack  Toplogies,”  Power  Engineering  Journal.  June  2001. 

[56]  MIL-STD-1399,  Interface  Standard  for  Shipboard  Systems.  Section  300A:  Electric 
Power.  Alternating  Current,  1987. 

[57]  Irving  H.  Shames  &  James  M.  Pitarresi,  Introduction  to  Solid  Mechanics,  3rd  Ed., 

Prentice  Hall,  2000. 

[58]  Website  www.specialmetals.com. 

[59]  Website  home.san.rr.com/nessengr/techdata/metalresis.html. 

[60]  Website  www.eddy-current.com/condres.htm. 


132 


[61]  Elsevier.  Carbon  Fiber  Composites.  1994. 

[62]  Website  www.grc.nasa.gov/www/epbranch/other/interlabs.htm. 

[63]  Michael  F.  Ashby  &  David  R.  H.  Jones,  Engineering  Materials.  Pergamon  Press,  1991. 

[64]  Wole  Soboyejo,  Mechanical  Properties  of  Engineering  Materials.  2003. 

[65]  James  L.  Kirtley  &  Edward  C.  Lovelace,  “Drag  Loss  in  Retaining  Rings  of  Permanent 
Magnet  Motors,”  SatCon  Technology  Corporation,  March,  2003. 

[66]  H.  Polinder  &  M.  J.  Hoeijmakers,  “Eddy-Current  Losses  in  the  Permanent  Magnets  of  a 
PM  Machine,”  EMD97,  Conference  Publication  No.  444,  1997. 

[67]  Website  www.powerdesigners.com/InfoWeb/design_center/articles/SSI/ssi.shtm. 

[68]  Kunrong  Wang,  Fred  C.  Lee,  Guichao  Jua,  and  Dusan  Borojevic,  “A  Comparative  Study 
of  Switching  Losses  of  IGBTs  Under  Hard-Switching,  Zero-Voltage-Switching,  and  Zero- 
Current-Switching,”  IEEE  Transactions.  1994. 

[69]  IEEE  Standards  Board,  IEEE  Recommended  Practice  for  Monitoring  Electric  Power 
Quality,  1995. 


133 


Page  Intentionally  Left  Blank 


134 


Appendix  A.  Detailed  Power  Requirements 


In  order  to  properly  size  the  PM  generator  being  designed,  a  typical  load  list  is  developed 
for  an  IPS  naval  ship.  For  the  propulsion  portion,  it  is  assumed  that  the  ship  contains  a  single 
shaft  requiring  30,000  SHP.  Shaft  horsepower  (SHP)  is  the  power  measured  in  the  shafting 
within  the  ship  by  a  torsion  meter  as  close  as  possible  to  the  propeller  or  stem  tube  [10].  This 
SHP  then  translates  to  32,552.1  BHP  (24,275  kW)  assuming  a  96%  efficient  electric  propulsion 
motor  and  a  96%  efficient  power  control  module  for  the  motor.  Brake  horsepower  (BHP)  is  the 
power  required  at  the  engines/generators.  It  is  also  assumed  that  the  ship  contains  two  equal¬ 
sized  generators.  Overall,  Table  21  lists  the  systems  and  their  associated  connected  (maximum) 
electrical  power  requirements. 


Table  21:  Ship  Connected  Loads 


na 

1 . — | 

Ilf  Engineering  Systems  j| 

Hydraulics 

150 

Engine  Room  Fresh  Water 

50 

Main  Sea  Water 

300 

Electronics  Fresh  Water 

25 

Aux  Sea  Water 

80 

Reactor  Cooling 

2,000 

Air  Conditioning 

300 

Reactor  Fresh  Water 

50 

Chill  Water 

50 

Reactor  Heating 

350 

Distilling 

20 

Steam  Feed  System 

300 

Lube  Oil 

50 

Steam  Condensing  System 

80  i 

Miscellaneous 

25 

||  Auxiliary  Systems  III 

Atmosphere  Control 

150 

Drain  System 

100 

Refrigeration 

10 

Ventilation 

250 

High  Pressure  Air 

100 

Damage  Control 

20 

Low  Pressure  Air 

25 

Sanitary  System 

10 

T rim  System 

100 

III  Combat  Systems  If 

IlSonar 

250 

Radio/ESM 

40  1! 

IlCombat  Control 

250 

Navigation 

20  || 

11  Hotel 

Loads  fi 

IlLighting 

60 

. . . . 

50  l! 

||Laundry 

20 

Galley 

70  ! 

III  .  ,  Total  Non-Propulsion 

III  Propulsion  ||| 

Main  Motor 

24,275 

Auxiliary  Propulsion 

300 

1 

SHSPBIII 
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Classic  U.S.  Navy  design  practices  (DDS  310-1)  determine  generator  sizes  by  first 
estimating  the  expected  load  during  the  preliminary  design  phases.  Propulsion  load  is  a  function 
of  the  hull  form,  maximum  speed,  and  drive  technology.  Ship  service  loads  are  a  function  of  the 
ship’s  mission,  crew  size,  and  ship’s  size.  Once  this  is  complete,  the  estimates  are  adjusted  for 
anticipated  construction  and  life  cycle  growth.  Propulsion  load  is  not  expected  to  change  but 
ship  service  loads  are  given  a  20%  margin  for  construction  growth  and  a  20%  margin  for  life 
cycle  growth  [11].  The  generators  are  then  sized  using  Eqn  A- 1,  where  Pgen  is  the  power 
provided  by  one  generator. 

_  (1.21.2  MaxLoad  ) 

Pgen  "  [0.9- (n  -  1)] 

where  n  =  number  of  generators 

Eqn  A- 1 

Applying  this  method  to  the  total  electrical  load  from  Table  21  produces  two  generators  that  are 
each  capable  of  47.9  MW. 

Since  IPS  ships  require  electrical  power  for  ship  loads  and  propulsion,  efficient  power 
management  should  be  implemented.  This  means  that  power'is  shared  among  all  systems  so  that 
at  top  speed,  more  power  is  diverted  to  the  propulsion  motor,  and  at  lower  speeds  power  is 
available  for  other  uses.  It  is  assumed  that  rated  ship  service  power  and  full  propulsion  power 
are  possible  at  the  same  time  only  if  both  PM  generators  are  available.  If  one  generator  is  lost, 
the  rated  ship  service  load  is  still  provided  with  the  remaining  power  sent  to  the  propulsion 
motor.  In  addition,  full  power  for  both  the  main  motor  and  auxiliary  propulsion  is  not  required 
simultaneously.  Similar  to  the  classical  method,  a  20%  margin  for  construction  growth  and  a 
20%  margin  for  life  cycle  growth  are  allotted.  Based  on  these  ideas,  Eqn  A-  2  is  developed  to 
size  the  generators,  where  Ptotai  is  the  total  generation  capacity  (both  generators). 

^  _  (1.21.2MaxNonPropulsionLoad)  +  MainMotorLoad 

Pto,al  =  ~ 

where  n  =  number  of  generators 

Eqn  A-  2 

Utilizing  this  process  produces  a  total  generation  capacity  of  3 1 .99  MW  or  two  generators  that 
are  each  capable  of  16  MW. 
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A  standard  U.S.  Navy  ship  spends  95%  of  its  time  at  less  than  full  speed  [12].  Because 
of  this,  the  classical  method  unnecessarily  over-estimates  the  power  requirements  for  the  ship 
and  results  in  generator  sets  that  are  too  large.  If  the  average  demand  for  each  system  is  also 
considered,  then  the  disparity  becomes  even  worse.  Table  22  lists  each  system  along  with  a 
demand  factor  which  represents  the  average  amount  of  power  required  for  normal  operation  of 
the  systems.  Based  on  the  9.8  MW  total  average  loads,  the  classical  method  generators  would 
each  be  loaded  to  only  10.2  %  capacity  which  is  extremely  inefficient. 

The  improved  generator  sizing  method  does  a  better  job  of  ensuring  efficient  operation 
while  not  excessively  limiting  the  ship  if  one  generator  is  damaged  or  down  for  maintenance.  If 
only  one  of  the  16  MW  generators  is  available,  the  connected  ship  service  load  (5.4  MW)  can  be 
supplied  while  still  having  10.6  MW  for  propulsion,  giving  44  %  propulsion  capability  which  is 
acceptable  when  operating  in  reduced  status. 
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Table  22:  Ship  Average  Loads 


I  ■'  ■  ;  '  V  ‘ 


Hydraulics 


Main  Sea  Water 


Aux  Sea  Water 


Chill  Water 


Distillin 


Lube  Oil 


Engine  Room  Fresh  Water 


Electronics  Fresh  Water 


Reactor  Coolin 


Reactor  Fresh  Water 


Reactor  Heatin 


Steam  Feed  System 


Steam  Condensing  System 


Miscellaneous 


Atmosphere  Control 


Refrigeration 


High  Pressure  Air 


Low  Pressure  Air 


T rim  System 


Drain  System 


Ventilation 


Damage  Control 


Sanitary  System 


Sonar 


Combat  Control 


Radio/ESM 


Navigation 


,■  "  'H 


Engineering  Systems 


150 

1.00 

150 

300 

0.20 

60 

80 

0.20 

16 

300 

0.75 

225 

50 

1.00 

50 

20 

!  0.50 

10 

50 

1.00 

50 

50 

1.00 

50 

25 

1.00 

25 

2,000 

0.25 

500 

50 

1.00 

50 

350 

0.25 

88 

300 

0.75 

225 

80 

1.00 

80 

25 

1.00 

25 

Liqhtin 


Laund 


Water  Heatin 


Galley 


Main  Motor 


Auxiliary  Propulsion 


Auxiliary  Systems 


150 


10 


100 


25 


100 


100 


250 


20 


10 


Combat  Systems 


250 


250 


40 


20 


Hotel  Loads 


60 


:.T~rrrr.,  i 


Propulsion 


24,275 


300 
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Appendix  B.  MATLAB  Code:  Basic  Sizing  Method 


Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  pmlbasic 

%  Program  performs  basic  sizing  and  parameter  calculations 
%  for  generators. 


%  Definition  of  variables 
%  Name 

%  General  variables 
%  Pwr 
%  rpm 
%  psi 
%f 


%  omega 
%  vtip 
%  LovD 
%  stress 

%  Rotor  variables 


Variable 


Required  power 
Speed  (RPM) 

Power  factor  angle 
Electrical  frequency  (Hz) 
Electrical  frequency  (rad/sec) 
Tip  speed  (m/s) 

L/D  ratio 

Gap  shear  stress  (psi) 


%  R 
%  D 
%  Lst 

%p 

%  Bg 

%  Stator  variables 
%  Kz 
%  Jz 
%  hs 


Rotor  radius  (m) 

Rotor  diameter  (m) 

Rotor  stack  length  (m) 

Number  of  pole  pairs 
Expected  air  gap  flux  density  (T) 

Surface  current  density  (A/m) 
Current  density  (A/m2) 

Slot  height  (m) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


clear; 

%  Constants  &  conversion  factors 


hs  =  .015; 
lams  =  0.5; 
convl  =9.81; 
conv2  =  703.0696; 


>  Assume  slot  depth  of  15  mm 

>  Assume  slot  fill  fraction 
-  9.81  W  per  Nm/s 

i  703.0696  N/m2  per  psi 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  INPUTS 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

Pwr  =  16e6;  %  Required  power 

vtip  =  200;  %  Max  tip  speed  (m/s) 

LovD  =  2.851;  %  Wound  rotor  usually  0.5-1.0,  PM  1.0-3.0  N 

%  Shear  stress  usually  1-10  psi  small  machines,  10-20  large  liquid 
%  liquid  cooled  machines 
stress  =15; 

p  =  3;  %  Pole  pairs 

Bg  =  0.8;  %  Tesla 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Calculations 
%  Size 
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%  Initially  use  Pwr  =  2*pi*R*Lst*stress*vtip 
%  Lst  =  2*LovD*R 
hscm  =  hs*100; 

R  =  sqrt(Pwr/(2*pi*(LovD*2)*vtip*stress*convl*conv2)); 

D  =  2*R; 

Lst  =  LovD*D; 

%  Speed 

omega  =  (p*vtip)/R; 
f  =  omega/(2*pi); 
rpm  =  (60*f)/p; 

%  Current  densities 

Kz  =  (stress*conv2)/(Bg*100); 

Ja  =  10*Kz/(hscm*lams); 

%  Output 

fprintf('Basic  Machine  DesignVn’); 
fprintf('Input  Parameters:\n'); 

fprintf('Power  =  %10.1f  kW  Shear  Stress  =  %10.1f  psi\n\Pwr/le3,  stress); 
f£>rintf('L/D  Ratio  =%10.2f  Tip  Speed  =  %10.1f  rn/s\n',LovD,vtip); 
f^>rintf(’Pole  Pairs  =  %10.1f  Air  Gap  Bg  =  %10.1f  TVn',p,Bg); 
fprintf('Output:\n'); 

fprintfORotor  Radius  =  %10.3f  m  Stack  Length  =  %10.3f  m\n',R,Lst); 
fprintf('Speed  =  %10.0f  RPM  Frequency  =  %  10. If Hz\n',rpm,f); 
fprintfCJa  =  %10.2f  A/cm2W,Ja); 
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Appendix  C.  PM  Machine  Database 
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Appendix  D.  MATLAB  Code:  Sizing  Method  1 


%  Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  pm  1  input 
%  Program  used  as  input  file  for  pmlcalc 
%  All  necessary  input  parameters  entered  here. 


clear; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Definition  &  Entry  of  variables 
%  General  variables 

Pwr  =  16e6;  %  Required  power  (W) 

rpm  =  1 3000;  %  Speed  (RPM) 

psi  =  0;  %  Power  factor  angle 

%  Rotor  variables 


R  =  0.147; 
hm  =  0.025; 

Lst  =  0.838; 
p  =  3; 

Br=  1.2; 
thm  =  50; 
thsk  =10; 

%  Stator  variables 
q  =  3; 

Ns  =  36; 

Nsp  =  1 ; 
g  =  .004; 
tfrac  =  0.5; 
hs  =  .025; 
hd  =  .0005; 
wd  =  le-6; 
syrat  =  0.7; 

Nc  =  1; 
lams  =  0.5; 
sigst  =  6.0e+7; 

%  Densities 
rhos  =  7700; 
rhom  =  7400; 
rhoc  =  8900; 


%  Rotor  radius  (m) 

%  Magnet  thickness  (m) 

%  Rotor  stack  length  (m) 

%  Number  of  pole  pairs 
%  Magnet  remnant  flux  density  (T) 

%  Magnet  physical  angle  (deg) 

%  Magnet  skew  angle  (actual  deg) 

%  Number  of  phases 
%  Number  of  slots 
%  Number  of  slots  short  pitched 
%  Air  gap  (m) 

%  Peripheral  tooth  fraction 
%  Slot  depth  (m) 

%  Slot  depression  depth  (m) 

%  Slot  depression  width  (m) 

%  Stator  back  iron  ratio  (yoke  thick/rotor  radius) 
%  Turns  per  coil 
%  Slot  fill  fraction 
%  Stator  winding  conductivity 

%  Steel  density  (kg/m3) 

%  Magnet  density  (kg/m3) 

%  Conductor  density  (kg/m3) 


%  Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  pmlcalc 

%  Program  performs  sizing  and  parameter  calculations 
%  for  permanent  magnet  machines  with  surface  magnets  and 
%  slotted  stators. 

%  Program  developed  from  J.L.  Kirtley  script  with  permission 
%  MUST  RUN  pml input  PRIOR  TO  RUNNING  pmlcalc 
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%  Definition  of  variables 
%  Name 

%  General  variables 
%  Pwr 
%  rpm 
%  psi 


%  omega 
%  vtip 
%  lambda 
%  Ea 

%  Rotor  variables 


Variable 


Required  power  (W) 

Speed  (RPM) 

Power  factor  angle 
Electrical  frequency  (Hz) 
Electrical  frequency  (rad/sec) 
Tip  speed  (m/s) 

Flux  linkage 

RMS  Internal  voltage  (V) 


%  R  Rotor  radius  (m) 

%  hm  Magnet  thickness  (m) 

%  Lst  Rotor  stack  length  (m) 

%  p  Number  of  pole  pairs 

%  Br  Magnet  remnant  flux  density  (T) 

%  thm  Magnet  physical  angle  (deg) 

%  thsk  Magnet  skew  angle  (actual  deg) 

%  Stator  variables 


%  q 
%  m 
%  Ns 
%  Nsp 

%g 
%  ge 
%  tfrac 
%  hs 
%  hd 
%  wd 
%  syrat 
%Nc 
%  lams 
%  sigst 
%  Kc 

%  Loss  Models 
%  PO  Base  power  for  core  losses 

%  FO  Base  frequency  for  core  loss 

%  BO  Base  flux  density 

%  epsb  Flux  density  exponent 

%  epsf  Frequency  exponent 

%  rhos  Steel  density 

%  rhom  Magnet  density 

%  rhoc  Conductor  density 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


Number  of  phases 

Slots  per  pole  per  phase 

Number  of  slots 

Number  of  slots  short  pitched 

Air  gap  (m) 

Effective  air  gap  (m) 

Peripheral  tooth  fraction 
Slot  depth  (m) 

Slot  depression  depth  (m) 

Slot  depression  width  (m) 

Stator  back  iron  ratio  (yoke  thick/rotor  radius) 

Turns  per  coil 

Slot  fill  fraction 

Stator  conductivity 

Carter  coefficient 


%  Constants  to  be  used 
muO  =  4*pi*le-7; 
tol  =  le-2; 
cpair=  1005.7; 
rhoair  =  L205  ; 
nuair=  1.5e-5; 

PO  =  36.79; 

FO  =  1000; 

B0  =1.0; 


%  Free  space  permeability 
%  Tolerance  factor 

%  Specific  heat  capacity  of  air  (J/kg*C) 

%  Density  of  air  at  20  C  (kg/m3) 

%  Kinematic  viscosity  of  air  at  20  C  (m2/s) 
%  Base  Power  Losss,  W/lb 
%  Base  freuency,  60  Hz 
%  Base  flux  density,  1.0  T 
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epsb  =  2.12; 
epsf=1.68; 

%  Generate  geometry  of  machine 

%  Number  of  slots/pole/phase 
m  =  Ns/(2*p*q); 

%  Number  of  armature  turns  (each  slot  has  2  half  coils) 

Na  =  2*p*m*Nc; 

%  Tooth  width 

wt=  2*pi*(R+g+hm+hd)*tfrac/Ns; 

%  Slot  top  width  (at  air  gap) 

wst  =  2*pi*(R+g+hm+hd)*(l-tfrac)/Ns; 

%  Slot  bottom  width 

wsb  =  wst*(R+g+hd+hs)/(R+g+hm+hd); 

%  Stator  core  back  iron  depth  (as  p  increases,  dc  decreases) 
dc  =  syrat*R/p; 

%  Full-pitch  coil  throw 
Nsfp  =  floor(Ns/(2*p)); 

%  Actual  coil  throw 
Nsct  =  Nsfp  -  Nsp; 

%  Estimate  end  turn  length 
%  End  turn  travel  (one  end) 
laz  =  pi*(R+g+hm+hd+0.5*hs)*Nsct/Ns; 

%  End  length  (half  coil) 
le2  =  pi*laz; 

%  End  length  (axial  direction) 
lei  =  2*le2/(2*pi); 

%  Calculate  electrical  frequency  &  surface  speed 
f  =  p*rpm/60; 
omega  =  2*pi*f; 
vtip  =  R*omega/p; 

%  Winding  &  skew  factors 

gama  =  2*pi*p/Ns; 

alfa  =  pi*Nsct/Nsfp; 

kp  =  sin(pi/2)*sin(alfa/2); 

kb  =  sin(m*gama/2)/(m*sin(gama/2)); 

kw  =  kp*kb; 

ths  =  ((p*thsk)+le-6)*(pi/180);  %  skew  angle  (elec  rad) 

ks  =  sin(ths/2)/(ths/2); 

%  Calculate  magnetic  gap  factor 
Rs  =  R+hm+g; 

Ri  =  R; 

R1  =  R; 

R2  =  R+hm; 

kg  =  ((RiA(p- 1  ))/(RsA(2*p)-Ri  A(2*p)))  *((p/(p+ 1  ))*(R2A(p+ 1  )-R  1  A(p+ 1 )). . . 
+(p*RsA(2*p)/(p- 1  ))*(R  1  A(  1  -p)-R2A(  1  -p))); 

%  Calculate  air  gap  magnetic  flux  density 

%  Account  for  slots,  reluctance,  and  leakage 
ws  =  (wst+wsb)/2;  %  Average  slot  width 

taus  =  ws  +  wt;  %  Width  of  slot  and  tooth 
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Kc  =  l/(l-(l/((taus/ws)*((5*g/ws)+l)))); 
ge  =  Kc*g; 

Cphi  =  (p*thm)/180;  %  Flux  concentration  factor 

K1  =  0.95;  %  Leakage  factor 

Kr  =  1.05;  %  Reluctance  factor 

murec  =1.05;  %  Recoil  permeability 

PC  =  hm/(ge*Cphi);  %  Permeance  coefficient 

Bg  =  ((Kl*Cphi)/(l+(Kr*murec/PC)))*Br; 

%  Calculate  magnetic  flux  and  internal  voltage 
thmrad  =  thm*(pi/180); 

B1  =  (4/pi)*Bg*kg*sin(p*thmrad/2); 
lambda  =  2*Rs*Lst*Na*kw*ks*B  1/p; 

Ea  =  omega*lambda/sqrt(2);  %  RMS  back  voltage 

%  Calculation  of  inductances/reactances 

%  Air-gap  inductance 

Lag  =  (q/2)*(4/pi)*(muO*NaA2*kwA2*Lst*Rs)/(pA2*(g+hm)); 

%  Slot  leakage  inductance 

perm  =  mu0*((l/3)*(hs/wst)  +  hd/wst); 

Las  =  2*p*Lst*perm*(4*NcA2*(m-Nsp)+2*Nsp*NcA2); 

Lam  =  2*p*Lst*Nsp*NcA2*perm; 

if  q  =  3 

Lslot  =  Las  +  2*Lam*cos(2*pi/q);  %  3  phase  equation 

else 

Lslot  =  Las  -  2*Lam*cos(2*pi/q);  %  multiple  phases 

end 

%  End-turn  inductance  (Hanselman) 

As  =  ws*hs;  %  Slot  area 

Le  =  ((Nc*muO*(taus)*NaA2)/2)*log(wt*sqrt(pi)/sqrt(2*As)); 

%  Total  inductance  and  reactance 
Ls  =  Lag+Lslot+Le; 

Xs  =  omega*Ls; 

%  Lengths,  Volumes,  and  Weights 

%  Armature  conductor  length 
Lac  =  2*Na*(Lst+2*le2); 

%  Armature  conductor  area  (assumes  form  wound) 

Aac  =  As*lams/(2*Nc); 

%  Mass  of  armature  conductor 
Mac  =  q*Lac*Aac*rhoc; 

%  Overall  machine  length 
Lmach  =  Lst+2*lel; 

%  Core  inside  radius 
Rci  =  R+hm+g+hd+hs; 

%  Core  outside  radius 
Rco  =  Rci+dc; 

%  Overall  diameter 
Dmach  =  2*Rco; 

%  Core  mass 

Mcb  =  rhos*pi*(RcoA2-RciA2)*Lst;  %  Back  iron 

Met  =  rhos*Lst*(Ns*wt*hs+2*pi*R*hd-Ns*hd*wd);  %  Teeth 

Me  =  Mcb  +  Met; 

%  Magnet  mass 
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Mm  =  0.5*(p*thmrad)*((R+hm)A2-RA2)*Lst*rhom; 

%  Shaft  mass 

Ms  =  pi*RA2*Lst*rhos; 

%  15%  service  fraction 

Mser  =  0.15*(Mc+Ms+Mm+Mac); 

%  Total  mass 

Mtot  =  Mser+Mc+Ms+Mm+Mac; 

%  Stator  resistance 
Ra  =  Lac/(sigst*Aac); 

%  Core  Loss  Calculations 

%  Tooth  Flux  Density 
Bt  =  Bg/tfrac; 

%  Back  iron  flux  density  (Hanselman) 

Bb  =  Bg*R/(p*dc); 

%  Core  back  iron  loss 

Pcb  =  Mcb*PO*abs(Bb/BO)Aepsb*abs(f/FO)Aepsf; 

%  Teeth  Loss 

Pet  =  Mct*PO*abs(Bt/BO)Aepsb*abs(f/FO)Aepsf; 

%  Total  core  loss 
Pc  =  Pcb  +  Pet; 

%  Start  loop  to  determine  terminal  voltage  and  current 
notdone  =  1; 

i  =  0; 

la  =  Pwr/(q*Ea); 
while  notdone  ==1 
i  =  i  +  1; 
xa  =  Xs*Ia/Ea; 

%  Conductor  losses 
Pa  =  q*IaA2*Ra; 


%  Gap  friction  losses 
%  Reynold's  number  in  air  gap 
omegam  =  omega/p; 

Rey  =  omegam*R*g/nuair; 

%  Friction  coefficient 
Cf  =  .0725/ReyA  2; 

%  Windage  losses 

Pwind  =  Cf*pi*rhoair*omegamA3*RA4*Lst; 

%  Get  terminal  voltage 

Va  =  sqrt(EaA2-((Xs+Ra)*Ia*cos(psi))A2)-(Xs+Ra)*Ia*sin(psi); 

Ptemp  =  q*Va*Ia*cos(psi)-Pwind; 
error  =  Pwr/Ptemp; 
err(i)  =  error; 
if  abs(error-l)  <  tol 
notdone  =  0; 
else 

la  =  Ia*error; 
end 
end 
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%  Remaining  performance  parameters 

%  Current  density 
Ja  =  Ia/Aac; 

%  Power  and  efficiency 
Pin  =  Pwr+Pc+Pa+Pwind; 
eff  =  Pwr/Pin; 
pf  =  cos(psi); 

fprintfOpmlcalc  complete:  Ready  An’); 


%  Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  pm  1  output 
%  Program  outputs  values  from  pm  1  calc. 

%  Program  developed  from  J.L.  Kirtley  script  with  permission 

%  MUST  RUN  pmlinput  and  pmlcalc  PRIOR  TO  RUNNING  pmloutput 


%  Variables  for  output  display 
Pout  =  Pwr/le3; 

Jao  =  Ja/le4; 

Pco  =  Pc/le3; 

Pwindo  =  Pwind/le3; 

Pao  =  Pa/le3; 
wso  =  ws*1000; 
hso  =  hs*1000; 
wto  =  wt*1000; 
dco  =  dc*1000; 

Lso  =  Ls*1000; 
hmo  =  hm*1000; 
go  =  g*1000; 


%  Output  Section: 

fprintf('\nPM  Machine  Design,  Version  1:  Surface  Magnet,  Slotted  StatorW); 


fprintf(’Machine  SizeAn’); 
fprintf('Machine  Diameter  =  %8.3f  m 
fprintf(’Rotor  radius  =  %8.3f  m 

fprintf('Slot  Avg  Width  =  %8.3f  mm 
fprintf('Back  Iron  Thick  =  %8.3f  mm 

fprintf('Machine  Ratings:\n'); 
fprintf('Power  Rating  =  %8.1f  kW 

fprintf('Va  (RMS)  =  %8.0f  V 

fprintf('Ea  (RMS)  =  %8.0f  V 

fprintf(fSynch  Reactance  =  %8.3f  ohm 
fprintf('Stator  Cur  Den  =  %8.1f  A/cm2 

fJ)rintfCEfficiency  =  %8.3f 

fprintf('Phases  =  %8.0f 

fprintf('Stator  ParametersAn’); 
fprintf('Number  of  Slots  =  %8.0f 


Machine  Length  =  %8.3f  m\n',Dmach,Lmach); 
Active  length  =  %8.3f  m\n',R,Lst); 

Slot  Height  =  %8.3f  mm\n',wso,hso); 

Tooth  Width  =  %8.3f  mm\n',dco,wto); 


Speed  =  %8.0f  RPM\n’,  Pout,rpm); 

Current  =  %8.1f  A\n',  Va,Ia); 

Arm  Resistance  =  %8.5f  ohm\n’,Ea,Ra); 
Synch  Induct  =  %8.3f  mH\n',Xs,Lso); 
Tip  Speed  =  %8.0f  m/s\n’,  Jao,vtip); 
Power  Factor  -  %8.3f\n',  eff,pf); 
Frequency  =  %8.  If  Hz\n’,q,f); 


Num  Arm  Turns  =  %8.0f  \n',Ns,Na); 
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fprintf('Breadth  Factor  =  %8.3f 

fprintf(Tooth  Flux  Den  =  %8.2f  T 

fprintf('Slots/pole/phase  =  %8.2f\n',m); 

fprintf('Rotor  Parameters  :\n'); 
f^rintf(’Magnet  Height  =  %8.2f  mm 

fprintf('Air  gap  =  %8.2f  mm 

fprintf('Magnet  Remanence  =  %8.2fT 
fprintf('Magnet  Factor  =  %8.3f 

fprintf('Machine  Losses:\n'); 
fprintf('Core  Loss  =  %8.1f  kW 

fprintf('Windage  Loss  =  %8.1f  kW 

fprintf('Machine  Weights  :\n’); 
fprintf('Core  =  %8.2fkg 

fprintf('Magnet  =  %8.2f  kg 

fprintf(’Services  =  %8.2f  kg 


Pitch  Factor  =  %8.3f  W,  kb,kp); 
Back  Iron  =  %8.2f  T\n’,  Bt,Bb); 


Magnet  Angle  =  %8.1f  degm\n',hmo,thm); 
Pole  Pairs  =  %8.0f  \n\go,p); 

Aig  Gap  Bg  =  %8.2f  T\n’,Br,Bg); 

Skew  Factor  =  %8.3f  \n‘,kg,ks); 


Armature  Loss  =  %8.1f  kWW,  Pco,Pao); 
Rotor  Loss  =  TBD  kWVn',  Pwindo); 


Shaft  =  %8.2f  kg\n’,Mc,Ms); 

Armature  =  %8.2f  kg\n’,Mm,Mac); 

Total  =  %8.2f  kg\n',Mser,Mtot); 
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Appendix  E.  MATLAB  Code:  Sizing  Method  2 


%  Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  pm2input 
%  Program  used  as  input  file  for  pm2calc 
%  All  necessary  input  parameters  entered  here. 

clear; 


%  Definition  &  Entry  of  variables 
%  General  variables 

Pwr  =  16e6;  %  Required  power  (W) 

rpm  =  1 3000;  %  Speed  (RPM) 

psi  =  0;  %  Power  factor  angle 

Bsat  =  1.65;  %  Stator  saturation  flux  density 

%  Rotor  variables 

vtip  =  200;  %  Tip  speed  limit  (m/s) 

p  =  3;  %  Number  of  pole  pairs 

Br  =  1 .2;  %  Magnet  remnant  flux  density  (T) 

thsk  =  10;  %  Magnet  skew  angle  (elec  deg) 

PC  =  5.74;  %  Permeance  coefficient  for  magnets 


%  Stator  variables 
Ja  =  2200; 
q  =  3; 
m  =  2; 

Nsp  =  1; 
g  =  .004; 
hs  =  .025; 
hd  =  .0005; 
wd  =  le-6; 
ws  =  .016; 

Nc  =  1; 
lams  =  0.5; 
sigst  =  6.0e+7; 

%  Densities 
rhos  =  7700; 
rhom  =  7400; 
rhoc  =  8900; 


%  Initial  current  density  (A/cm2) 
%  Number  of  phases 
%  Slots/pole/phase 
%  Number  of  slots  short  pitched 
%  Air  gap  (m) 

%  Slot  depth  (m) 

%  Slot  depression  depth  (m) 

%  Slot  depression  width  (m) 

%  Avg  slot  width  (m) 

%  Turns  per  coil 
%  Slot  fill  fraction 
%  Stator  winding  conductivity 

%  Steel  density  (kg/m3) 

%  Magnet  density  (kg/m3) 

%  Conductor  density  (kg/m3) 


%  Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  pm2calc 

%  Program  performs  sizing  and  parameter  calculations 
%  for  permanent  magnet  machines  with  surface  magnets  and 
%  slotted  stators. 

%  MUST  RUN  pm2input  PRIOR  TO  RUNNING  pm2calc 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Definition  of  variables 
%  Name  Variable 


151 


%  General  variables 

%  Pwr 

Required  power  (W) 

%  rpm 

Speed  (RPM) 

%  psi 

Power  factor  angle 

%f 

Electrical  frequency  (Hz) 

%  omega 

Electrical  frequency  (rad/sec) 

%  vtip 

Tip  speed  (m/s) 

%  lambda 

Flux  linkage 

%  Ea 

%  Rotor  variables 

RMS  Internal  voltage  (V) 

%  R 

Rotor  radius  (m) 

%  hm 

Magnet  thickness  (m) 

%  Lst 

Rotor  stack  length  (m) 

%p 

Number  of  pole  pairs 

%  Br 

Magnet  remnant  flux  density  (T) 

%  thm 

Magnet  physical  angle  (deg) 

%  thsk 

%  Stator  variables 

Magnet  skew  angle  (actual  deg) 

%q 

Number  of  phases 

%  m 

Slots  per  pole  per  phase 

%  Ns 

Number  of  slots 

%  Nsp 

Number  of  slots  short  pitched 

%g 

Air  gap  (m) 

%  ge 

Effective  air  gap  (m) 

%  tfrac 

Peripheral  tooth  fraction 

%  hs 

Slot  depth  (m) 

%hd 

Slot  depression  depth  (m) 

%  wd 

Slot  depression  width  (m) 

%  syrat 

Stator  back  iron  ratio  (yoke  thick/rotor  radius) 

%Nc 

Turns  per  coil 

%  lams 

Slot  fill  fraction 

%  sigst 

Stator  conductivity 

%  Kc 

%  Loss  Models 

Carter  coefficient 

%P0 

Base  power  for  core  losses 

%F0 

Base  frequency  for  core  loss 

%  BO 

Base  flux  density 

%  epsb 

Flux  density  exponent 

%  epsf 

Frequency  exponent 

%  rhos 

Steel  density 

%  rhom 

Magnet  density 

%  rhoc 

Conductor  density 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Constants  to  be  used 

muO  =  4*pi*le-7; 

%  Free  space  permeability 

tol  =  le-2; 

%  Tolerance  factor 

cpair  =  1005.7; 

%  Specific  heat  capacity  of  air  (J/kg*C) 

rhoair=  1.205  ; 

%  Density  of  air  at  20  C  (kg/m3) 

nuair  =  1.5e-5; 

%  Kinematic  viscosity  of  air  at  20  C  (m2/s) 

P0  =  36.79; 

%  Base  Power  Losss,  W/lb 

F0  = 1000; 

%  Base  freuency,  60  Hz 

B0  =  1.0; 
epsb  =  2.12; 
epsf  =  1.68; 

%  Base  flux  density,  1.0  T 
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%  Calculate  electrical  frequency  &  rotor  radius 
f  =  p*rpm/60; 
omega  =  2*pi*f; 

R  =  p*vtip/omega; 

%  Winding  &  skew  factors 

Ns  =  floor(2*q*p*m);  %  Number  of  slots 

gama  =  2*pi*p/Ns; 

Nsfp  =  floor(Ns/(2*p)); 

Nsct  =  Nsfp  -  Nsp; 

alfa  =  pi*Nsct/Nsfp; 

kp  =  sin(pi/2)*sin(alfa/2); 

kb  =  sin(m*gama/2)/(m*sin(gama/2)); 

kw  =  kp*kb; 

ths  =  ((p*thsk)+le-6)*(pi/180);  .  %  skew  angle  (elec  rad) 

ks  =  sin(ths/2)/(ths/2); 

%  Calculate  magnet  dimensions,  tooth  width,  &  air  gap  flux  density 
thme  =1;  %  Initial  Magnet  angle  (deg  e) 

notdone  =  1 ; 

ge  =  g;  %  Initial  effective  air  gap 

while  notdone  ==  1 

alpham  =  thme/180;  %  Pitch  coverage  coefficient 

Cphi  =  (2*alpham)/(l+alpham);  %  Flux  concentration  factor 
hm  =  ge*Cphi*PC;  %  Magnet  height 

Ds  =  2*(R+hm+g);  %  Inner  stator/air  gap  diameter 

K1  =  0.95;  %  Leakage  factor 

Kr  =  1.05;  %  Reluctance  factor 

murec  =1.05;  %  Recoil  permeability 

Bg  =  ((Kl*Cphi)/(l+(Kr*murec/PC)))*Br; 
wt  =  ((pi*Ds)/Ns)*(Bg/Bsat);  %  Tooth  width 
taus  =  ws  +  wt;  %  Width  of  slot  and  tooth 

Kc=  l/(l-(l/((taus/ws)*((5*g/ws)+l))));  %  Carter's  coefficient 
ge  =  Kc*g; 
eratio  =  ws/wt; 
if  abs(eratio  -  1)  <  tol 
notdone  =  0; 
else 

thme  =  thme  +  1 ; 
end 
end 

%  Set  final  values 
thm  =  thme/p; 
thmrad  =  thm*(pi/180); 
hm  =  ge*Cphi*PC; 

Ds  =  2*(R+hm+g); 

%  Generate  geometry  of  machine 

%  Peripheral  tooth  fraction 
tfrac  =  wt/(wt+ws); 

%  Slot  top  width  (at  air  gap) 
wst=  2*pi*(R+g+hm+hd)*tfrac/Ns; 

%  Slot  bottom  width 

wsb  =  wst*(R+g+hm+hd+hs)/(R+g+hm+hd); 


%  Magnet  physical  angle 

%  Magnet  height 
%  Inner  stator/air  gap  diameter 
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%  Stator  core  back  iron  depth 

dc  =  (pi*Ds*thmrad/(4*p))*(Bg/Bsat); 

%  Core  inside  radius 
Rci  =  R+hm+g+hd+hs; 

%  Core  outside  radius 
Rco  =  Rci+dc; 

%  Slot  area 
As  =  ws*hs; 

%  Estimate  end  turn  length 
%  End  turn  travel  (one  end) 
laz  =  pi*(R+g+hm+hd+0.5*hs)*Nsct/Ns; 

%  End  length  (half  coil) 
le2  =  pi*laz; 

%  End  length  (axial  direction) 
lei  =  2*le2/(2*pi); 

%  Calculate  magnetic  gap  factor 
Rs  =  R+hm+g; 

Ri  =  R; 

R1  =  R; 

R2  =  R+hm; 

kg  =  ((Ri  A(p- 1  ))/(Rs  A(2  *p)-Ri  A(2  *p)))  *((p/(p+ 1  ))*(R2  A(p+ 1  )-R  1  A(p+ 1 )) . . . 
+(p*RsA(2*p)/(p- 1  ))*(R  1  A(  1  -p)-R2A(  1  -p))) ; 

%  Core  loss  calculations  (per  length) 

%  Core  mass  per  length 
McbperL  =  rhos*pi*(RcoA2-RciA2); 

MctperL  =  rhos*(Ns*wt*hs+2*pi*R*hd~Ns*hd*wd); 

McperL  =  McbperL  +  MctperL; 

%  Tooth  Flux  Density 
Bt  =  Bg/tfrac; 

%  Back  iron  flux  density  (Hanselman) 

Bb  =  Bg*R/(p*dc); 

%  Core  back  iron  loss  per  length 

PcbperL  =  McbperL*P0*abs(Bb/B0)Aepsb*abs(CT?0)Aepsf; 

%  Teeth  Loss  per  length 

PctperL  =  MctperL*PO*abs(Bt/BO)Aepsb*abs(f/FO)Aepsf; 

%  Total  core  loss  per  length 
PcperL  =  PcbperL  +  PctperL; 

%  Current  and  surface  current  density 
%  Armature  turns  (each  slot  has  2  half  coils) 

Na  =  2*p*m*Nc; 

%  Arm  cond  area  (assumes  form  wound) 

Aac  =  (As*lams)/(2*Nc); 

%  Power  &  Current  waveform  factors  (Lipo) 
ke  =  0.52; 
ki  =  sqrt(2); 

%  Initial  terminal  current 
la  =  Ns*lams*As*Ja*le4/(2*q*Na); 
notfin  =  1; 

Lst  =  0. 1 ;  %  Initial  stack  length 

i  =  l; 


%  Back  iron 
%  Teeth 
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%  Start  loop  to  determine  Lst,  Ea,  Va,  and  la 
notdone  =  1 ; 
k  =  0; 

while  notdone  ==  1 
k  =  k+  1; 

%  Surface  current  density 
A  =  2*q*Na*Ia/(pi*Ds); 


%  Calculate  stack  length  of  machine 
%  Loop  to  get  stack  length 
while  notfin  ==  1 
%  Gap  power 

Pgap  =  4*pi*ke*ki*kw%s*kg*sin(thmrad)*(f/p)*A*Bg*(DsA2)*Lst; 
%  Length  of  conductor 
Lac  =  2*Na*(Lst+2*le2); 

%  Stator  resistance 
Ra  =  Lac/(sigst*Aac); 

%  Copper  Loss 
Pa  =  q*IaA2*Ra; 

%  Core  losses 
Pc  =  PcperL*Lst; 

%  Iterate  to  get  length 
Ptempl  =  Pgap-Pa-Pc; 
error  =  Pwr/Ptempl; 
err(i)  =  error; 
if  abs(error-l)  <  toi 
notfin  =  0; 
else 

Lst  =  Lst*error; 
i  =  i  +  1; 
end 
end 

%  Calculate  magnetic  flux  and  internal  voltage 
thmrad  =  thm*(pi/180); 

B1  =  (4/pi)*Bg*kg*sin(p*thmrad/2); 
lambda  =  2*Rs*Lst*Na*kw*ks*Bl/p; 

Ea  =  omega*lambda/sqrt(2);  %  RMS  back  voltage 

%  Calculation  of  inductances/reactances 

%  Air-gap  inductance 

Lag  =  (q/2)*(4/pi)*(muO*NaA2*kwA2*Lst*Rs)/(pA2*(g+hm)); 

%  Slot  leakage  inductance 

perm  =  mu0*((l/3)*(hs/wst)  +  hd/wst); 

Las  =  2*p*Lst*perm*(4*NcA2*(m-Nsp)+2*Nsp*NcA2); 

Lam  =  2*p*Lst*Nsp*NcA2*perm; 
if  q  ==  3 

Lslot  =  Las  +  2*Lam*cos(2*pi/q);  %  3  phase  equation 

else 

Lslot  =  Las  -  2*Lam*cos(2*pi/q);  %  multiple  phases 

end 

%  End-turn  inductance  (Hanselman) 

taus  =  ws  +  wt;  %  Width  of  slot  and  tooth 

Le  =  ((Nc*muO*(taus)*NaA2)/2)*log(wt*sqrt(pi)/sqrt(2*As)); 
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%  Total  inductance  and  reactance 
Ls  =  Lag+Lslot+Le; 

Xs  =  omega*Ls; 

%  Lengths,  Volumes,  and  Weights 

%  Armature  conductor  length 
Lac  =  2*Na*(Lst+2*le2); 

%  Mass  of  armature  conductor 
Mac  =  q*Lac*Aac*rhoc; 

%  Overall  machine  length 
Lmach  =  Lst+2*lel; 

%  Overall  diameter 
Dmach  =  2*Rco; 

%  Core  mass 
Me  =  McperL*Lst; 

%  Magnet  mass 

Mm  =  0.5*(p*thmrad)*((R+hm)A2-RA2)*Lst*rhom; 

%  Shaft  mass 

Ms  =  pi*RA2*Lst*rhos; 

%  15%  service  fraction 

Mser  =  0.15*(Mc+Ms+Mm+Mac); 

%  Total  mass 

Mtot  =  Mser+Mc+Ms+Mm+Mac; 

%  Gap  friction  losses 
%  Reynold's  number  in  air  gap 
omegam  =  omega/p; 

Rey  =  omegam*R*g/nuair; 

%  Friction  coefficient 
Cf  =  .0725/ReyA  2; 

%  Windage  losses 

Pwind  =  Cf*pi*rhoair*omegamA3*RA4*Lst; 

%  Get  terminal  voltage 
xa  =  Xs*Ia/Ea; 

Va  =  sqrt(EaA2-((Xs+Ra)*Ia*cos(psi))A2)-(Xs+Ra)*Ia*sin(psi); 

Ptemp  =  q*Va*Ia*cos(psi)-Pwind; 

Perror  =  Pwr/Ptemp; 

Perr(k)  =  Perror; 
if  abs(Perror-l)  <  tol 
notdone  =  0; 
else 

la  =  Ia*Perror; 
end 
end 

%  Remaining  performance  parameters 

%  Current  density 
Ja  =  Ia/Aac; 

%  Power  and  efficiency 
Pin  =  Pwr+Pc+Pa+Pwind; 
eff  =  Pwr/Pin; 
pf  =  cos(psi); 
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fprintf('pm2calc  complete:  Ready  An'); 


%  Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  pm2output 
%  Program  outputs  values  from  pm2calc. 

%  Program  developed  from  J.L.  Kirtley  script  with  permission 

%  MUST  RUN  pm2input  and  pm2calc  PRIOR  TO  RUNNING  pm2output 

%  Variables  for  output  display 
Pout  =  Pwr/le3; 

Jao  =  Ja/le4; 

Pco  =  Pc/le3; 

Pwindo  =  Pwind/le3; 

Pao  =  Pa/le3; 
wso  =  ws*1000; 
hso  =  hs*1000; 
wto  =  wt*1000; 
dco  =  dc*1000; 

Lso  =  Ls*1000; 
hmo  =  hm*1000; 
go  =  g*1000; 

%  Output  Section: 

fprintf('\nPM  Machine  Design,  Version  2:  Surface  Magnet,  Slotted  StatorW); 
fprintf(Machine  Size:\n'); 

fprintf(Machine  Diameter  =  %8.3f  m  Machine  Length  =  %8.3f  m\n\Dmach,Lmach); 

fprintf('Rotor  radius  =  %8.3f  m  Active  length  =  %8.3f  m\n’,R,Lst); 

fprintf('Slot  Avg  Width  =  %8.3f  mm  Slot  Height  =  %8.3f  mm\n’,wso,hso); 

fprintf('Back  Iron  Thick  =  %8.3f  mm  Tooth  Width  =  %8.3f  mm\n',dco,wto); 

fprintf(Machine  Ratings:\n'); 

f^rintf('Power  Rating  =  %8. If  kW  Speed  =  %8.0f  RPM\n',  Pout,rpm); 

fprintf('Va  (RMS)  =  %8.0fV  Current  =  %8.1f  A\n’,  Va,Ia); 

fprintf('Ea  (RMS)  =  %8.0f  V  Arm  Resistance  =  %8.5f  ohm\n',Ea,Ra); 

fprintf(’Synch  Reactance  =  %8.3f  ohm  Synch  Induct  =  %8.3f  mH\n’,Xs,Lso); 
fprintf(’Stator  Cur  Den  =  %8.1fA/cm2  Tip  Speed  =  %8.0f  m/sVn’,  Jao,vtip); 

fprintf(’Efficiency  =  %8.3f  Power  Factor  =  %8.3f\n',  eff,pf); 

fprintf(’Phases  =  %8.0f  Frequency  =  %8.1f  Hz\n',q,f); 

fprintfOStator  Parameters  :\n’); 

fprintf('Number  of  Slots  =  %8.0f  Num  Arm  Turns  =  %8.0f  \n’,Ns,Na); 

fprintf('Breadth  Factor  =  %8.3f  Pitch  Factor  =  %8.3f  \n',  kb,kp); 

fprintfCTooth  Flux  Den  =  %8.2fT  Back  Iron  =  %8.2f  T\n’,  Bt,Bb); 

fprint^’Slots/pole/phase  =  %8.2f\n’,m); 

fprintf('Rotor  Parameters:W); 

fprintf(Magnet  Height  =  %8.2f  mm  Magnet  Angle  =  %8.1f  degm\n',hmo,thm); 

fprintf('Air  gap  =  %8.2f  mm  Pole  Pairs  =  %8.0f  W,go,p); 

fprintf(Magnet  Remanence  =  %8.2fT  AigGapBg^  %8.2f  T\n',Br,Bg); 
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fprintf('Magnet  Factor  =  %8.3f 


Skew  Factor  =  %8.3f  Vn',kg,ks); 


fprintf('Machine  Losses  :\n'); 

fprintf('Core  Loss  =  %8.1fkW  Armature  Loss  =  %8.1f  kWW,  Pco,Pao); 

fprintf('Windage  Loss  =  %8.1f  kW  Rotor  Loss  =  TBD  kWW,  Pwindo); 

fprintf(’Machine  Weights:\n'); 


fprintf(’Core  = 

%8.2f  kg 

Shaft  = 

%8.2f  kg\n',Mc,Ms); 

fprintf('Magnet  = 

%8.2f  kg 

Armature  = 

%8.2f  kg\n\Mm,Mac); 

fprintf(’Services  = 

%8.2f  kg 

Total  = 

%8.2f  kg\n\Mser,Mtot); 
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Appendix  F.  MATLAB  Code:  Bode  Plot 


%  Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  Buckfilter 

%  Program  calculates  transfer  function  and  outputs 
%  Bode  plot  for  buck  converter  input  filter 

clear; 

%  Input  parameters 
R  =  4.79e-3; 

Cf=2.84e-3; 

Cb  =  28.4e-3; 

Lf  =  1.415e-6; 

%  Set  up  transfer  function 
num=[R*Cb  1]; 

den  =  [Lf*R*Cf*Cb  Lf*(Cf+Cb)  R*Cb  1]; 

H  =  tf(num,den); 
bode(H) 
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Appendix  G. 

MATLAB  Code:  PM  Generator  Waveforms 

%  Jonathan  Rucker,  MIT  Thesis 

%  May  2005 

%  Program:  pmwave 

%  Program  calculates  and  outputs  different  waveforms, 

%  calculates  THD,  and  computes  the  harmonic  content 

%  for  permanent  magnet  machines  with  surface  magnets  and 

%  slotted  stators. 

%  MUST  RUN  pm  1  input  and  pm  1  calc  PRIOR  TO  RUNNING  pmwave 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Definition  of  variables 

%  Name 

Variable 

%  General  variables 

%  Pwr 

Required  power  (W) 

%  rpm 

Speed  (RPM) 

%  psi 

Power  factor  angle 

%f 

Electrical  frequency  (Hz) 

%  omega 

Electrical  frequency  (rad/sec) 

%  vtip 

Tip  speed  (m/s) 

%  lambda 

Flux  linkage 

%  Ea 

RMS  Internal  voltage  (V) 

%  Rotor  variables 

%  R 

Rotor  radius  (m) 

%  hm 

Magnet  thickness  (m) 

%  Lst 

Rotor  stack  length  (m) 

%p 

Number  of  pole  pairs 

%  Br 

Magnet  remnant  flux  density  (T) 

%  thm 

Magnet  physical  angle  (deg) 

%  thsk 

Magnet  skew  angle  (actual  deg) 

%  Stator  variables 

%q 

Number  of  phases 

%  m 

Slots  per  pole  per  phase 

%  Ns 

Number  of  slots 

%  Nsp 

Number  of  slots  short  pitched 

%g 

Air  gap  (m) 

%  ge 

Effective  air  gap  (m) 

%  tfrac 

Peripheral  tooth  fraction 

%  hs 

Slot  depth  (m) 

%  hd 

Slot  depression  depth  (m) 

%  wd 

Slot  depression  width  (m) 

%  syrat 

Stator  back  iron  ratio  (yoke  thick/rotor  radius) 

%  Nc 

Turns  per  coil 

%  lams 

Slot  fill  fraction 

%  sigst 

Stator  conductivity 

%  Kc 

Carter  coefficient 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Constants  to  be  used 

muO  =  4*pi*le-7; 

%  Free  space  permeability 

tol  =  le-2; 

%  Tolerance  factor 
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%  Harmonics  to  be  evaluated 
n=  1:2:35; 

np  =  p  .*  n;  %  Use  in  kgn  equation 

w  =  n  .*  omega;  %  Harmonic  angular  frequencies 

ffeq  =  w  ./  (2*pi);  %  Harmonic  frequencies 

%  Harmonic  winding  and  skew  factors 

gama  =  2*pi*p/Ns; 

alfa  =  pi*Nsct/Nsfp; 

kpn  =  sin(n  .*  pi/2)  .*  sin(n  .*  alfa/2); 

kbn  =  sin(n  .*  m*gama/2) ./  (m*sin(n  .*  gama/2)); 

kwn  =  kpn  .*  kbn; 

ths  =  ((p*thsk)+le-6)*(pi/180);  %  skew  angle  (elec  rad) 

ksn  =  sin(n  .*  ths/2) ./  (n  .*  ths/2); 

%  Calculate  magnetic  gap  factor 
Rs  =  R+hm+g; 

Ri  =  R; 

R1  =  R; 

R2  =  R+hm; 

kgn  =  ((Ri .  A(np- 1  ))./(Rs.  A(2.  *np)-Ri.  A(2.  *np))).*((np./(np+ 1 )). .. 
.*(R2.A(np+l)-Rl.A(np+l))+(np.*Rs.A(2.*np)./(np-l))... 

.*(R1  .A(  1  -np)-R2.A(l  -np))); 

%  Calculate  magnetic  flux  and  internal  voltage 
thmrad  -  thm*(pi/180); 
thmerad  =  p*thmrad; 

Bn  =  Bg.*((4/pi)./n).*kgn.*sin(n.*thmerad/2).*sin(n.*pi/2); 
lambdan  =  ((2*Rs*Lst*Na).*kwn.*ksn.*Bn)./p; 

Ean  =  (omega. *lambdan);  %  Peak  back  voltage 

%  Normalized  values  for  plotting 
Eanorm  =  abs(Ean) ./  Ean(l); 

%  Voltage  THD 
Eah  =  0; 

for  r  =  2:length(n) 

Eah  =  Eah  +  Ean(r)A2; 
end 

THD=  100*sqrt(Eah/(Ean(l)A2)); 

%  Generate  waveforms 

%  Rotor  physical  angle  goes  from  0  to  2*pi  -  electrical  to  2*p*pi 
ang  =  0:pi/100:2*pi; 
angp  =  p*ang; 

Bout  =  zeros(size(angp)); 

Eaout  =  zeros(size(angp)); 
for  i  =  l:length(n) 

Bout  =  Bout  +  Bn(i).*sin(n(i).*angp); 

Eaout  =  Eaout  +  Ean(i).*sin(n(i).*angp); 
end 

%  Plot  waveforms 

figure(l) 

plot(ang,Bout); 

title('PM  Generator:  Flux  Density'); 
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ylabel('B  (Tesla)1); 
xlabel('Rotor  Angle  (rad)'); 

figure(2) 

plot(ang,Eaout); 

title('PM  Generator:  Back  EMF’j; 
ylabelCPeak  Voltage  (V)1); 
xlabel('Rotor  Angle  (rad)'); 

figure(3) 
hold  on 

title(['PM  Generator:  EMF  Harmonics,  THD  =  ',num2str(THD),‘  %']); 
ylabel('Normalized  Back  EMF'); 
xlabel('Harmonic  Number'); 

text(20, 0.7, ’Dark:  Above  10%  of  Fundamental', 'FontSize',  10); 
text(20, 0.65, ’Light:  Below  10%  of  Fundamental', ‘FontSize’,  10); 
for  z  =  l:length(Eanorm) 
if  Eanorm(z)  <0.10 
bar(n(z),Eanorm(z),'c'); 
else 

bar(n(z),Eanorm(z),'b’); 

end 

end 

hold  off 
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Appendix  H.  MATLAB  Code:  Retaining  Sleeve  Stress  Calculations 


%  Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  pmcanstress 

%  Program  calculates  and  outputs  retaining  can  stress 
%  for  permanent  magnet  machines  with  surface  magnets  and 
%  slotted  stators. 

%  MUST  RUN  pmlinput,  pmlcalc,  and  pmwave  PRIOR  TO  RUNNING  pmcanstress 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Calculate  retaining  sleeve  stresses 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Conversion 

Patopsi  =  1.45038e-4;  %  psi  per  Pa 

%  Material  yield  stresses  (ksi) 

Stain_str  =  90; 

Alum„str  =  75; 

Titan_str  =  110; 

CarFib_str  =  100; 

Inconel_str  =  132; 

%  Safety  factor 
SF=  1.2; 

%  Force  on  magnets/sleeves  is  centrifugal  force 
%  Magnet  tangential  velocity 
vmag  =  ((R+hm)*omega)/p; 

%  Centrifugal  force 

Fm  =  (Mm*vmagA2)/(R+hm); 

%  Outward  pressure 

Phoop  =  Fn^(2*pi*(R+hm)*Lst); 

%  Hoop  Stress  (in  general,  str  =  P*R/t) 
stop  =  22; 
for  i  =  l:stop 

t(i)  =  i*.0005;  %  sleeve  thickness  t 

slev(i)  =  t(i)*1000; 

Sthoop(i)  =  (Phoop*(R+hm)/t(i))*Patopsi/1000; 

SFHoop(i)  =  Sthoop(i)*SF; 
end 

%  Output  results 

fprintf('Retaining  Sleeve  Stress:\n'); 
fprintf(’Stress  LimitsAn’); 

fprintf('Stainless  Steel  =  %6.1f  ksi  Aluminum  Alloy  =  %6.1f  ksi\n’,Stain_str,Alum_str); 
fprintf(’Titanium  Alloy  =  %6.1fksi  Carbon  Fiber  =  %6.1f  ksi\n,,Titan__str,CarFib_str); 
fprintf(’Inconel  =  %6. 1  f  ksi\n',Inconel_str); 

fprintf(Actual  Sleeve  Stress:\n’); 

fprintf('Sleeve  Thickness  Actual  Stress  SF  StressW); 
for  i  =  l:stop 

fprintf('  %5.2fmm  %6.1fksi  %6.1f ksi\n',... 
slev(i),Sthoop(i),SFHoop(i)); 
end 
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Appendix  I.  MATLAB  Code:  Rotor  Losses  from  Winding  Time  and 
Space  Harmonics 

%  Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  pmharmloss 

%  Program  performs  rotor  loss  calculations  caused  by 
%  winding  time  and  space  harmonics  for  permanent  magnet 
%  machines  with  surface  magnets  and  slotted  stators. 

%  MUST  RUN  pmlinput,  pmlcalc,  pmloutput,  and  get  harmonic 
%  current  data  from  PSIM  prior  to  running  pmharmloss 

%  Constants  to  be  used 

muO  =  4*pi*le-7;  %  Free  space  permeability 

tol  =  le-2;  %  Tolerance  factor 

%  Retaining  sleeve/magnet  material  resistivity  (ohm-m) 

Stain_res  =  0.72e-6; 

Titan_res  =  0.78e-6; 

CarFib_res  =  9.25e-6; 

Inconel_res  =  0.98e-6; 

Magnet_res  =  1.43e-6; 

%  Retaining  sleeve  thickness  set  at  0.5mm  less  than  air  gap 
h_sl  =  g  -  0.0005;  %  Sleeve  thickness 

g_act  =  g  -  h_sl;  %  Actual  air  gap 

%  Retaining  sleeve  conductivities  (S/m) 

cond„s  =  1/Stain_res; 

cond_t  =  1/Titanjres; 

cond„c  =  l/CarFib_res; 

condji  =  l/Inconel_res; 

%  Magnet  &  actual  sleeve  cond  (S/m) 
cond_m  =  1/Magnet_res; 
cond_sl  =  cond_s; 

%  Input  time  harmonic  peak  currents  from  PSEM 
II  =2895; 

13  =  0; 

15  =  209; 

17  =  89.2; 

19  =  0; 

111=39.2; 

113  =  27.6; 

115=0; 

117  =  17.0; 

119=12.8; 

121  =0; 

123  =  7.3; 

125  =  6.3; 

127  =  0; 
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129  =  4.8; 

131=3.9; 

%  Put  currents  in  array 

Iharm  =  [II  13  15  17  19  II 1  113  115  117  119 121  123  125... 

127  129  131]; 

%  Calculate  current  THD 
Iah  =  0; 

for  r  =  2:length(Iharm) 

Iah  =  Iah  +  Iharm(r)A2; 
end 

THDi  =  1 00*sqrt(lah/(lharm(  1 )  A2)) ; 

%  Calculate  current  densities 
Iz  =  (l/sqrt(2)).*Iharm; 

Kz  =  ((q/2)*(Na/(2*pi*Rs))).*Iz; 

Iz _ 1  =  (l/sqrt(2)).*Il; 

Kz_l  =  ((q/2)*(Na/(2*pi*Rs))).*Iz_l; 

%  Harmonics  to  be  evaluated 
n=  1:2:31;  " 
w  =  n  .*  omega; 
freq  =  w  ./  (2*pi); 
lam  =  (2*(2*pi/(2*p)))./n; 
k  =  (2*pi)./lam; 

%  Eta  values 
eta_m  =  sqrt((j*muO*cond_m).*w  +  (k.A2)); 
eta_s  =  sqrt((j*muO*cond_sl).*w  +  (k.A2)); 

%  Surface  coefficient  at  top  of  magnet  layer 
alpha_m  =  j .  *(k./eta_m) .  *coth(eta_m.  *hm) ; 

%  Surface  coefficient  at  top  of  retaining  sleeve 

topi  =  (j.*(k./eta_s).*sinh(eta„s.*h_sl))  +  (alpha_m.*cosh(eta_s.*h_sl)); 
botl  =  (j.*(k./eta_s).*cosh(eta_s.*h„sl))  +  (alpha_m.*sinh(eta_s.*h_sl)); 
alpha_s  =  j .  *(k./eta_s).  *(top  1  ./bot  1 ); 

%  Surface  coefficient  at  surface  of  stator 
top2  =  (j.*sinh(k.*g_act))  +  (alpha_s.*cosh(k.*g_act)); 
bot2  =  (j.*cosh(k.*g_act))  +  (alpha„s.*sinh(k.*g_act)); 
alpha_f  =  j.*(top2./bot2); 

%  Surface  impedance 
Zs  =  (muO.*w./k).*alpha_f; 

%  Calculate  losses  due  to  time  harmonics 
%  Use  only  fundamental  space  harmonic  factors 
Kz_t  =  kw.*Kz; 

Syt  =  0; 

for  i  =  l:length(n) 

Sy_t(i)  =  0.5*(abs(Kz_t(i))A2)*real(Zs(l)); 

Syt  =  Syt  +  Sy_t(i); 
end 


%  Harmonic  angular  frequencies 
%  Harmonic  frequencies 

%  Wavenumbers 


%  Fundamental  RMS  current 
%  Fundamental  current  density 
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%  Calculate  losses  due  to  space  harmonics 
%  Use  only  fundamental  time  harmonic  current 
kpn  =  sin(n  .*  pi/2)  .*  sin(n  .*  alfa/2); 
kbn  =  sin(n  .*  m*gama/2) ./  (m.*sin(n  .*  gama/2)); 
kwn  =  kpn  .*  kbn; 

Kz_s  =  kwn  .*  Kz_l  ./  n; 

Sys  =  0; 

for  i  =  l:length(n) 

Sy_s(i)  =  0.5*(abs(Kz_s(i))A2)*real(Zs(i)); 

Sys  =  Sys  +  Sy„s(i); 
end 

fprintf(’\nRotor  Losses  Caused  by  Harmonics:\n'); 
fprintf('Time  Harmonic  Losses  =  %6.1f  kW\n\Syt/1000); 
f^rintf('Space  Harmonic  Losses  =  %6.1f  kW\n’,Sys/1000); 
fprintf(Total  Losses  =  %6.  If  kW\n’,(Syt+Sys)/1000); 

fprintf('Current  THD  =  %6.2f  %%\n’,THDi); 
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Appendix  J.  MATLAB  Code:  Rotor  Losses  from  Slot  Effects 


%  Jonathan  Rucker,  MIT  Thesis 
%  May  2005 
%  Program:  pmcanloss 

%  Program  calculates  and  outputs  rotor  losses  caused  by 
%  stator  slot  effects  for  permanent  magnet  machines 
%  with  surface  magnets  and  slotted  stators. 

%  MUST  RUN  pmlinput,  pmlcalc,  pmwave,  and  pmcanstress 
%  PRIOR  TO  RUNNING  pmcanloss 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Calculate  retaining  sleeve  losses 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Retaining  sleeve/magnet  material  resistivity  (ohm-m) 

Stain_res  =  0.72e-6; 

Titan_res  =  0.78e-6; 

CarFib_res  =  9.25e-6; 

Inconel_res  =  0.98e-6; 

Magnet_res  =  1.43e-6; 

%  Calculate  Bd  as  function  of  wst  and  wt  (max  10%  of  Bg) 

Bd  =  (wst/wt)*0.1*Bg; 

%  Calculate  flux  variation  parameters 
beta  =  (wst/(2*pi*Rs))*2*pi; 
lamB  =  2*pi/Ns; 

B  =  (Bd/sqrt(2))*sqrt(beta/lamB); 

%  Calculate  geometry  and  can  loss  factor  for  different  rings 
%  k  is  number  of  rings 
for  k=  1:10 

A(k)  =  pi*2*(R+hm)*Lst/k; 

Ks(k)  =  1  -  ((tanh(p*Lst/(k*2*(R4-hm))))/(p*Lst/(k*2*(R+hm)))); 
end 

%  Input  Stainless  Steel  sleeve  thickness  based  on  stress  results 
for  i  =  l:stop 
if  SFHoop(i)  <=  Stain_str 
t_Stain  =  t(i); 
break 

elseif  t(stop)  >  Stain_str 

fprintf('Hoop  Stress  too  high  for  Stainless  Steel.Xn'); 
else 

dummy  =  t(i); 
end 
end 

%  Input  Titanium  sleeve  thickness  based  on  stress  results 
for  i  =  1  :stop 

if  SFHoop(i)  <=  Titan_str 
t_Titan  =  t(i); 
break 

elseif  t(stop)  >  Titan„str 
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fprintf('Hoop  Stress  too  high  for  Titanium  An’); 
else 

dummy  =  t(i); 
end 
end 

%  Input  Carbon  Fiber  sleeve  thickness  based  on  stress  results 
fori=l:stop 

if  SFHoop(i)  <=  CarFib„str 
t_CarFib  =  t(i); 
break 

elseif  t(stop)  >  CarFib_str 
fprintf('Hoop  Stress  too  high  for  Carbon  Fiber  An'); 
else 

dummy  =  t(i); 
end 
end 

%  Input  Inconel  sleeve  thickness  based  on  stress  results 
for  i  =  l:stop 

if  SFHoop(i)  <-  Inconel_str 
t_Inconel  =  t(i); 
break 

elseif  t(stop)  >  Inconel_str 

fprintf('Hoop  Stress  too  high  for  Inconel.Vn'); 
else 

dummy  =  t(i); 
end 
end 

%  Calculate  can  losses 

w_Stain  =  (piA2/3600)*((B*rpm*(R+hm))A2*tJStain)/Stain_res; 
wJTitan  =  (piA2/3600)*((B*rpm*(R+hm))A2*t_Titan)/Titan_res; 
w_CarFib  =  (piA2/3600)*((B*rpm*(R+hm))A2*t_CarFib)/CarFib_res; 
wjnconel  =  (piA2/3600)*((B*rpm*(R+hm))A2*t_Inconel)/Inconel_res; 

for  k=  1:10 

P_Stain(k)  =  k*w_Stain*Ks(k)*A(k)/1000; 

P_Titan(k)  =  k*w_Titan*Ks(k)*A(k)/1000; 

P_CarFib(k)  =  k*w^CarFib*Ks(k)*A(k)/1000; 

P_Inconel(k)  =  k*w„Inconel*Ks(k)*A(k)/1000; 
end 

%  Calculate  magnet  losses  (only  with  carbon  steel) 

%  Calculate  geometry  and  can  loss  factor 
Am  =  pi*2*R*Lst; 

Ksm  =  1  -  ((tanh(p*Lst/(2*R)))/(p*Lst/(2*R»); 

%  Calculate  magnet  losses 

%  Assumes  only  10%  of  magnet  thickness  sees  effects 
w_Magnet  =  (piA2/3600)*((B*rpm*R)A2*0.1*hm)/Magnet_res; 
PJVIagnet  =  w_Magnet*Ksm*  Am/1 000; 

%  Output  Results 
z  =  [l  5  10]; 

fprintf('Retaining  Can  LossesAn'); 
for  i=  1:3 
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k  =  z(i); 

fprintf('%  1  .Of  Rings:\n',k) 

fprintf('Material  Thickness  Can  Loss\n'); 

fprintf('Stainless  Steel  %5.2f  mm  %6.1f  kW\n,,t_Stain*1000,P_Stain(k)); 
fprintf('Titanium  %5.2f  mm  %6.1f  kW\n',t_Titan*1000,P_Titan(k)); 

fprintf( 'Carbon  Fiber  %5.2f  mm  %6.1f  kW\n’,t_CarFib*1000,P„CarFib(k)); 

fprintf(’  Associated  Magnet  Loss  %6. If  kW\n',P_Magnet); 

fprintf(Tnconel  %5.2f  mm  %6.1f  kW\n\n’,t„Inconel*1000,P_Inconel(k)); 

end 
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Appendix  K.  Results  for  PM  Machine  Variants 


Main  Parameters 
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Appendix  L.  Results  for  Power  Conversion  Module  Variants 
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Appendix  M.  Results  for  Power  Module  Losses 
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Appendix  N.  Results  for  Power  Module  Weights 


Power  Conversion  Module 
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Appendix  O.  Rectifier/Input  Filter  Mass  and  Volume  Calculations 
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Appendix  P.  Converter/Output  Filter  Mass  and  Volume  Calculations 
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Enargy  (J)  J  Mass  (kg)  Currant 


BEilEBB  MiiEB 

M7IF7M  WVWM J 


Enargy  (J)  Mass  (kg) 


Mass  (kg)  Modules  [Mass  (kg) |  Mass  (> 


E3EEEB  EEBB  BESB  fck?»i  ggggiB 
BEOSM  M'XiTrf  E CT  B  E3M 

BEEE3B — ■  EM  BEBB  EBEBEI  K  i  | 

BECB  ^ESBB  B EB||  B,|:*IM 

Ega  BBffiBB  BB3SDB  BBEEB  gTOM  fc  ggg  B  MB2BB1 

lKiM»B  BEEBM  MjEM  Btn/if  wi’i&ZM  BBBBBBBBBBBBBB^M 

m^ubb  m&im  mx^m  w&mm  BEgga  giru  i 

^«m.»  mttm  mxfim  ^rrga  bc^pbi  b  v.x*  km  I 

wmm  bmb  mjrim  wksm  w*xmm  b-mb  mmm  musem  i 

mMttm  liTiTM  BigjgB  mzx&fM  b^^B  ME2EM I 

«<MiW  BKEM*1^*  WfFYim  wwvirn  WfgjgjM  ■Ksuti  i 


■'KM| 

■ami  wEmm  bem  iugg 

g  i<;  n  bthtomi  Em  ■>*)-■  BrR?jra 

■cltltfJM  MutTif  MM  ■*t«  liltf 

|  rn  vi*m  BigggB  |M|M^ 
mii.wtu  m&Jcm  |E  SB  B^33| 

m HgB  BPt»^B  Wm^M  B  !••!  B  H  KB 

mmm  BEE  ME 

unr.y’f^^g^rgM^i.-MiELii^ 
m  Mtldl  mkkMM 

m  g  w.w  MEBM  ■tgW»tflil 

wvzxm  MinyitM  m&xim  m>rum  mnxm 

mslem  mmim  mki-mmma-m 

gi.THga  mupfwm  m^>rm  memm  BEEjgM 


BEKa 

_  UwK&swm 

m'.+Am  BD^ng  m&EM  ■»gMgfraii 
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Tola) 

PC.  sg  aJSL.  *£L  %r  ST  JS* 
Can  ^  (mF)  Enacyy  (J)  ^  Curmit  uH) 


I  UBEB  ■  1  gMI  BEM  CEggi  jHH  B5HB  1EE3-1  ■ 


Number 

ot 

Modules 

Total 

Inductor 

Volum. 

(m*) 

Total 

Capacitor 

Voluma 

(m5) 

Total  IQBT 
Voluma 
(m*) 

Total 

Dlod. 

Voluma 

(m*) 

Sarvlc. 

Volum. 

(m*) 

Total 

Volum. 

(m1) 

■gpElMIjlMBjgM|3j^EjjiaBgEgBBgBgl«lB|gpaijXjnpM 
B  ■ran  ^<n»y  ■»>/•■«  ia«y«iB  Egga  KEE1  TO|  ■*i»im  WSSEM B 

■  ME3jl  ■?»*  B  Kft«X»Il  EgM  KEga  ■  El  ■«*‘i  [>:•*  B»T«S  fcM  ■ 

■  Mrffiif  EggW  ■HM  »L»^o:«f  K  EB1  ■'Kin.-a  K»I*»<M  | 

BKgCTBHj5M|MiMIggj5^ga|gMS£MKgElggMEjngjllgjiEIli 
B  mhezm  ^ttttib  ggjggg  Mg  grvrrMMfr^  grnrrrn  ggggl  i 

B  «♦:■»♦■  ntl~M  wm  grrr^i  g^rera  m-eeiem  ggggl  mwzm  mrxKm  gntsTK7;!  b 
I  gg3j»  ^-x-'-r-m  ■«! f?m  wifrTrm  wziesm  WEBM  MEM  gr»irei  nr»rei  b 

I  «gM  mnfrm  rr»:*:»n  Eggi  ■gaB  »  E  Ijga  ■  gag  Eg  ^'-m  Eg  nm  b 

I  KEgjl  m£  M  B*gM  BOxof  Egja  Bggg  KEEB  BEB  Eggj  BfrM  tol  ■ 
I  »+>:♦*«■  pi.-'M  WMUfFM  lOxoni  mz  Eg  B  gjggl  BEB  BEggl  E»iE?il  B 

■  Bssaa  ms aw  bsizb  Essm  eeebi  km  Bgaa  bbebb  eseii  iggEa  ■ 


I  »  wxm  BEX'X’iM  MEEDM  rtfejQxffli 

■I 1  BB  BIBB  1 1  BBI  Mill  H 


BHajgTrny^^— B 

B  mmjM  Eggg  tir.m  ■■v.'.m  mai  Bg^as 
BMEBEigIlB*E?MM<MiMBlEM^Jiai 
■  ■KaBg^MBEEJBm^iBE^lB^ii 


I  ■  iri;jw  B^wifaM  B^l  B3§B  Bfl£3ujffl 


mm 


B  »  eb  r.i>.:«iH  tettm  mbm  UFmm 

■  BES3B  EjS*OB  WRil  MEM  BE EM  BBEEM 


«—  R*r°  -»>»  ss 


Total 
cap^itor 

Velum.  ***  .  ^ 

*  Currwit  [  uH) 


Total  Moduk  Modulo  Number  Total 

Inductor  IQBT  Diodo  Induct) 

Volum.  VoKinta  Volum.  OT  Volum 

(m*)  (m*)  (m*)  Modules  (m^ 


Inductor  C^Sor  ™  ,QBT  SSL  £*»  VT*W 

Volum*  Volum.  Volum.  V#,U**  Y*** 

(m*)  (m*)  (m^  (m*)  <m  >  «■* 


B  gggE|  MEEOBi  PinTiiB  M&gfM  »k!:g  ■  MlnB  ^irrm  iggEjl  B 

■  BP<*-^-M  WSHSM  mt  gap  Mfztm  WIZM  grrrFii  ■ 

B  Eggjl  P^M  BjEJ  ^oura  Ws*:-.-7l  fcoiw  MB  ■!•?!  ■■x>i^:a  KggEl  B 

■  |M[  jM^aw  mmm  M£ -*l*  ■n;j<  ■■K^hT-igM^a  ■ 

B  mgM  mnfrm  I Kt<>:»n  g^a|  |^E||^^a  i 

■  ggaai  M^irnM  bcm  eim  eo  bq^ikj  ■ 

B  fcV>T)l  TT7»  ^T'V.’il  ITF-tI  ■:ri^M  K^THM  ■>•'■  ■iT.TFl  ■ 

b  Egga  bp<^'1  bbe^jm  gy^i  gaga  |e  iru  mrim  jEgngi  ■ 

■  E3E1 BKEBBB  BBGSSfli  EBEffll  BBSIiUl  BBESM  M.  BB  BES5E3I ESE3I M 

B  |^E|  M3M  »ilr7iM  MOSM  gJ.H.M  r*T.IM=l  j 

■  BBEEMI  MgaBI  MSga  Bpsa  BBP^B  gEE3  gEE3 1 

I  fcCT.l  »TP«  f'T.V.TB  Ivj^tl  »:■>!<■  KWT-l  TM  I'T.irr.l  f )T.TkT.l  ■ 


b  ■£*£■  gasggj  wcEum  wbm  mm  BKBa 
B  tm!LiM  ■■i>:oi  ■  MPTiM  BEEDl 


■  BK  EM  IEE3  B  BBQF*1  ■»^:-]Bi  BEEM  H^^BI 

BME3IK!ESIlMi^MWBEMMi2BBBBB6ga 

■  BBDESM  ESSi  1 KEM  ISEl  BBSM  mwim 

Mm*iutn^mm>r-rimm>^:.mmrrmmj^sa 

BB^JESMBKE^MBEMBBI^BBEia 

■  ■KM  BEESMBBEM  WEEM  BB3^BI  MFKnia 
B  MEM  BE 

■  MK-Vill  fctT.ytiM  MF7?M  »T?M  Wg^r?M 


_  .  .  _  .  ,  Total  Total  Modulo  Modulo  Wumhor  Tot,lJ  Total 

Output  °  °  Capacitor  *f*°*  Total  "^t*1  Inductor  IQBT  Oiod.  Inductor  Capacitor 

Po'9*  Voli,.  c*7d“"  fP^“:  Volum.  y»*.  "’*'<*“•<  JCrtw  Volum.  Volum.  Volum.  ®»  Volum.  Volum. 

(mF)  ErvmvylJ)  ^  Cum»l  uH)  En«»M  (m,  Modulo. 


Sarvlc.  Total 

Volum.  Volum. 
(m*)  (m*) 


I  fc<  ■;■■  f'Tn^  BTyx-a  UTii  n.-jr-i  b 

■  KS3EBI  WEESam  BBBKM  EM  EEa  BE  MM  BCTM  M3B  BEEI 53  EEE1 B 

■  BBSMB  EEM  M££M  ■  EC^i  ms'*n  KESCT  M53B  ■■y«it»!  1  r>T»TFcM  | 

B  ■.i.i:^i[!.iiiBf’''Bi)i)ii):;M  E3E3|  B 

■  fcvr.1  ^TTW  f'Trrin  ■?:"«  ar^i  mi;,^  iriirr.!  i.t.^.i  b 

b  mmam  msm  M^p  bm b  b 

■  miotrn  BBI  - 1  \\ m  ■  r*T)Trni  wy*:*M  MEiEM  BBI  *  MjEEEjl  EH  k>1  B 

B  mSEM  miwsm  fEESSM  Igm  |M  ■■!>>  r«M  B 

B  gZSjgl  WESEM  km  Mtzww  m  f;j1  Eg  Eg  n.i  w.m  a 

B  ^p|  BBBQBJI  *«»/<■  ■  Ejgjgj  ME  g  g  ■  »g  »Ji  BBI  ^  1 B 

B  WE2EM  mxsm  MSEM  gKTira  jj^^M  |^M  BKB  wn^m  1 

B  E3E1  MjSM  BB^^M  KE|jM  mizM  fcl3E3  EEE1 1 

■  rai  B^^ap  BBEjEM  ^BH  ECT  KHEB  «5B  BQEDE3  EEE3 1 


B  M  wm  BEEBM  M'^^JB  BE^M  jEQgjj  MP^HH 


■  ■  P'»M  E  »V»i  ■  MF1RJ  MEM  MSEM  W3&m 


i  BB0E31 BEESOI WEEM  BES31  MEM  MBE5BB 


■  ■  EM  EQ3  *  BBE^M  BB3B3B  KgSM  BKEUI 


Output  ToUI  ToUI  Capacitor  Av*r*9#  ToUI 

Voltag*  <^dUnc  Capacitor  ^  Output  toM-~( 
(mF)  Energy  (J)  Current  uH) 


"“T  Numb.1  ^T°“  Total  IQBT  Smvlc.  ToUI 

'°BT  .1”°*  ot  S?"  Cy~r  Volum.  “1*  volum.  volum. 


Volum.  Volum.  Volum. 

(m*)  (m*)  (m*) 


ol  Volum.  Volum. 

Modules  (m»j  (m*j 


pg*T«iE5Mg»iriTr=iB 

giirm  b 

B  ^ElJBB  *»>!■■  grrrm  grrrrrr^  HiiFi  b 

j  BgjjgjB  MESjISMB  W2EEM  grrr<Ta  pi  rv'f j|  giTirm  gggEj  B 

B  *  i  ■  ■  y  ■  1  fltZM  BEgg  WEEM  ■■!■»■■  nr>^i  B 

B  fcj^MB  — <*  1 E  ■  BEl  K»K«l»ni  KK^i  MEM  jjgjgEjl  j 

B  ttSaj  fc  j-T  gjlKiM  E3|gj  B  ^  j 

B  E5E|  MBE5DBB  nrrrm  K^bb^b  nrirrriM  gTtik'.i  a 

BE5E1BI  i  1 

B  M<ll*  V  B«*^»  B'Xv<»j1  Bgga  g*T«?r«:;B  K»M  ■ 

B  >/»>!»■  BEH31  ■»»^1  gBiVtii  BiMtttfa  gB  EgB  B1 1  f M  a 

B  KE5E|  BBEHiaB  M^^B  1  ■.-*■  sm =■  »tsh.^ibB  KjB  TTTTrrri  w*tx  e I 

I  gEM  fcjtif  B3GBI1  B^^M  W^M  B  EEIEi  IMCii  B 

B  grora  mufwm  gyron  WZESM  MSE^M  mwzm  ■■T.Tr.’i  MSS3EM  B 

■  BESEB  MB'  1 1  •  MUBE1EB  EE^II^Pi  PE331  Kg^gM K»T»lt«M  BEHEiB 


B  B  EM  gi«>T»T  MRiM  M3BM 

j  m  iti-»  BOEB  BjHjB  ^EEH 

■BEajigS^i^BBE^iBgEgiEiai 


BBEMBEE5IBM!l3BBB!EBBgEBg6ga 


l^mzm  ■■n:t>i  ■  BBEy?™  ■t^-^  gHgUj 

B  MBiM  gy>v«iBMRii  M»y>M  ■;r^«  w^q 


B  ikf»  g^roni  m>vum  m%zmmT5mM®g}M 

B  BftiB  M  m»k.m  mmm  bti  BEEHB 

immmEssmMEmMEmmSMMM 


Ph..«  R*^ln»  Pol..  c^SL.  c^L  ZZZ  u*S2L< 

Can  V  oi  tag.  {mF)  Erwgy  (J)  Currant  uH) 


Total  Modul.  Modul.  Nljmbf.r  Total  Total 

Inductor  IQBT  Diod.  Inductor  Capacitor 

Volum.  Volum.  Volum.  °*  Volum.  Volum. 

(m*)  (m*)  (m*)  Modules  (m>)  {m>j 


BBE3BJEE5IM  BEEBB»E3iJBEBBiEM 


Sarvlc.  Total 

Volum.  Volum. 
(m’)  (m*) 


B  BBE3B  B^EB  W&MM  ■.^.^■■rT<,r/ca  a> B  TTrirm  b 

B  H'i>i'ip»«n  K'X»y»iB^>JLi^j  M38&M  BEgj  »ir'M  ■■y'H'Mgnr»;?<'-i  | 
B  1 1 J  i  ■  ga^B  MM  m g g B  gKnm  ■  t*M  U<i gjj  ■ 

■  EaElMEEggB  BEB  BBE5I1  E5E3  BEEEH  B 

B  Kgjjjl  gr>v.;  ■  jviJi^i  mzzwm  m^-.bm  m  w*m  w^TKm  yr?rm  a 

B  «■  BBK3B BSjEB »v<< ■  T  '  B l7T>ikM  B 

■  Kjggi  BEHiBB  DrrrrrB  ezi-a 3  E  * — BEjM  BKP1B  TTrsim  b 

B  ■<»:<■  Egg  i  ^  ikhi  b>  kM  a 

B  >/■>!»■  wta  B«x»y»ni  j^=^j  wlv-il  ■  kfti  ■*?■!  nm  nr>5  ki  b 

■  »«:«!■■  MgEM  PP’CTi  ■kpx»  ict<r»  m  Tyrrm  rrttic^  a 

B  |^CT  BE3I  »■!  :'»i  ffTfiM  i  Egm  » it  ii  ti  b>.-«  griirm  nr.Tfv'i  a 
B  MEE5M  ■ritPJB  ijEsgi  m- g  f  krti  t  girnirn  nriiitM  a 

B  bb^tm  ggjggi  f&mM  ■  i<?y  M  g'T»Tr»Ti  rn.i  ^  a 

B  p<^x^  M^TI  Egga  m  t.-rm  M-l  K'7'l  Q 1 

B  BEEUB  t=> BEM  biim  Triirm  g^r^kM  j 

I  >/■>][♦■  bum  groitia  Egm  »  a  a  g|  j^M  ■«x*sr»ia  ^  k.'^a  ■ 


B  mgSB  MF^M  ■>^a  MM 

BUIfff'V'lMgEM>KJB^WE3i 
■  BftjBBMiSjJB^IMBigaBBEiBBBBEai 


■  » ■  BE^*M  1»T-!  M 


3fEv1 


aBEBEaB 


BM^jjE^B3E3BBjgsBB3iiBB13g 

B  Bftaa  B[*K’M  BB*KJ  B*EM  BtfjLlM 


QEjBBSEBIEgBB 


BBEB  BEES3  B  BEE5B  gE3BI  BB^M  BEBSH 

1^ \wEzm 
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